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= concentration  variance  of  data  with  high  frequency 

fluctuations  removed  by  averaging  over  time  intervals  of 
duration  Ta,  kg^/m6 


= "effective"  concentration  variance,  reduced  from  the  true 
variance  c'2  by  a receptor  system  time  constant  tc,  see 

( A26 ) ; kg2/m6 

= mean  concentration  at  a fixed  height  z on  plume  centerplane 
y = 0,  kg/m3 

= maximum  ground  level  concentration  on  plume  centerline,  for 
fixed  plume  height  h and  windspeed  U,  see  (B4),  kg/m3 

= conditional  "inside  plume"  mean  excluding  periods  of 
zero  concentration;  kg/m3 


= conditional  "inside  plume"  fluctuation  variance  excluding 
periods  of  zero  concentration;  kg2/m6 
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= ensemble  variance  of  dose_  d of  individual  realizations  about 
the  ensemble  mean  dose,  D,  see  (E10);  kg2-s2/m6 

= constants  in  power  law  plume  spread  equations  (B2)(B3) 


Ec  = spectral  density  of  concentration  variance:  Ec(f)df  is  the 

~2 

amount  of  variance  c between  frequencies  f and  f+df; 
kg2-s/m6 
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; thickness  of  the  neutrally  stable  atmospheric  boundary  layer 
that  is  influenced  by  surface  roughness;  m 

h 

height  of  plume  centerline  above  ground;  m 

hv 

height  of  the  effective  variance  source  above  ground,  see 
( B20) ; m 

i=/c,2/c  = 

fluctuation  intensity  of  total  concentration  variance 

V^V 

conditional  fluctuation  intensity  excluding  periods  of 
zero  concentration. 

^ 00  * ^ poo 

total  and  conditional  fluctuation  intensities  measured  on 
the  plume  centerline,  y=0,  z=h  for  a plume  with  negligible 
surface  interaction,  az<<h 

1 D 

dosage  fluctuation  intensity,  see  E( 19 ) 

kl 

2-nf/U,  wavenumber,  a normalized  fluctuation  frequency;  m“* 

Pc 

probability  density  function  (pdf):  pc(c)dc  is  the 
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source  strength;  the  pollutant  release  rate;  kg/s 
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hazard  rate:  probability  per  unit  time  of  exposure  of 

exceeding  a specified  threshold  c*  for  the  first  time 

Rc 

normalized  autocorelation  coefficient,  see  (A7) 

S 

probability  of  not  exceeding  the  threshold  concentration  c* 
during  an  exposure  time  Te 
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= time  delay  in  autocorrelation,  see  (A7);  s 

= data  averaging  time  chosen  to  smooth  out  high  frequency 

o 

fluctuations  in  the  true  variance  c'  to  reduce  it  to 
c'g,  see  (A30);  s 

= exposure  time  over  which  a receptor  is  subjected  to 
concentrations;  s 

= sampling  time:  the  continuous  time  period  over  which  data 

is  collected,  and  which  sets  a lower  limit  on  the 
frequency  f of  fluctuations  that  can  contribute  to  the 

observed  variance  c'^,  and  mean  c-j.;  s 
= local  mean  wind  velocity;  m/s 

= mean  advection  (or  convection)  speed  of  a plume  past  a y-z 
plane;  m/s 

= along-wind  turbulent  velocity  variance;  m^/s^ 

= fraction  of  concentration  readings  (or  fraction  of  time) 
that  concentration  is  greater  than  the  specified  value  c* 
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= cross-wind  turbulent  velocity  variance;  m^/s^ 
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the  effective  variance  source,  see  (B20) 

intermittency  factor:  the  fraction  of  non-zero 

concentration  readings  in  a sample  at  a fixed  point  in 
space. 
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spread 

normalized  distance  at  which  maximum  ground  level 
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source,  see  (B36) 
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change  in  concentration;  s 

integral  Eulerian  length  scale  of  concentration 
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= integral  concentration  length  scale  at  a height  far  above 

2 

influence  of  surface  turbulence  and  wind  shear  on  c , 
when  oz«  h;  m 

VVAw 

= along-wind,  crosswind  and  vertical  Eulerian  integral 
scales  of  turbulent  velocity  fluctuations;  m 

Q 

= cumulative  probability  density  function;  q(c*)  is  the 
fraction  of  concentration  readings  less  than  a specified 
value  c*  see  (F2) 
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= integral  Eulerian  time  scale  of  concentration 
fluctuations;  see  (A8);  s 
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= integral  time  scale  of  concentration  fluctuations  with 
high  frequencies  removed  by  averaging  over  time  intervals 
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= apparent  integral  time  scale  of  concentration  fluctuations 
that  have  had  their  high  frequency  components  attenuated 
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Subscripts 


a = for  data  that  has  been  subjected  to  averaging  over  time 

blocks  of  length  Ta,  which  removes  contributions  from 
high  frequency  fluctuations 

e = for  a dose  exposure  time  Te  (this  is  equivalent  to  an 

averaging  time  process  with  Ta  e Te) 

eff  = "effective"  value,  with  high  frequency  fluctuations 

attenuated  by  the  receptor  system  response  time  constant 
Ts 

m = at  the  point  of  maximum  ground  level  concentration  for  a 

given  source  height  h and  windspeed  U 

o = measured  on  the  center-plane  of  the  plume,  i.e.  at  y = 0 

p = "inside  plume",  conditionally  sampled  value  with  periods 

of  zero  concentration  removed 

T = for  data  taken  with  a sampling  time  duration  Ts,  which 

excludes  contributions  from  low  frequency  fluctuations 
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EXECUTIVE  SUMMARY 


Objectives 

The  primary  objective  of  the  studies  begun  in  1979  has  been  to  develop 
a reliable  and  easy  to  use  prediction  method  for  estimating  ground  level 
mean  and  peak  concentrations  of  toxic  gas  from  the  rupture  of  sour  gas 
pipelines  and  well  blowouts.  The  reports  "Expansion  and  Plume  Rise  of  Gas 
Jets  from  High  Pressure  Gasline  Ruptures",  (April  1981)  and  the  report 
"Predicting  Risk  of  Exposure  to  Peak  Concentrations  in  Fluctuating  Plumes", 
(December  1982)  together  provide  the  theoretical  basis  for  this  model,  and 
verify  it  using  both  wind  tunnel  and  atmospheric  data.  The  present  study 
was  initiated  to  carry  out  more  fundamental  research  into  the  statistical 
theory  in  order  to  improve  the  model  and  allow  it  to  be  applied  to  a wider 
range  of  atmospheric  and  source  conditions.  Three  objectives  were  set: 


1.  To  develop  a form  for  the  fluctuating  concentration  model  that 
would  be  suitable  for  predicting  the  probability  of  peak 
concentrations  under  continuous  stack  gas  plumes.  Extending  the 
model  to  the  longer  exposure  times  that  are  associated  with 
continuous  plumes  required  methods  for  dealing  with  the  increase 
in  plume  meandering  with  exposure  time,  the  intermittent  periods 
of  zero  concentration  that  occur  due  to  meandering  and  incomplete 
mixing,  the  response  time  constants  of  biological  receptors,  and 
the  effect  of  atmospheric  stability  on  concentration 

fl uctuations. 

2.  To  use  the  existing  concentration  fluctuation  model  as  a method 
for  predicting  the  mean  of  all  non-zero  values  of  concentration 
at  a receptor.  This  conditional  mean  concentration  removes 
intermittent  periods  when  plume  meandering  reduces  the 
concentration  to  zero  by  moving  the  plume  completely  off  a 
receptor.  Because  this  conditional  mean  concentration  is  a more 
accurate  indicator  of  the  actual  dilution  within  the  meandering 
plume,  and  also  is  less  sensitive  to  varying  exposure  time,  its 
use  as  a superior  indicator  of  ambient  air  quality  compliance  was 
studied. 

3.  To  develop  new  statistical  techniques  which  deal  efficiently  with 
the  time  varying  risk  of  exposure  for  a fixed  receptor  under  a 
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transient  release  from  a pipeline  rupture.  Statistical  theories 
were  developed  to  predict  the  probability  of  being  exposed  to  a 
toxic  gas  plume  for  a fixed  time  without  ever  crossing  some  set 
concentration  threshold  c*.  Atmospheric  and  wind  tunnel  data  on 
concentration  fluctuations  were  used  to  develop  an  improved 
probability  distribution  function  to  describe  fluctuations  in 
intermittent  plumes.  Simple,  closed-form  equations  for  the 
threshold  crossing  probability  were  developed  to  reduce  the  time- 
consuming  numerical  integrations  previously  used  to  evaluate  risk 
of  exposure. 

The  emphasis  in  the  report  is  to  develop  statistical  and  diffusion  models 
which  form  the  basic  elements  of  a risk  analysis  calculation  procedure. 
Because  of  this  theoretical  orientation,  the  70  pages  in  the  body  of  the 
report  are  supplemented  by  93  pages  of  material  in  the  appendices  A - G 
show  the  various  theoretical  models  in  sufficient  detail  to  make  their 
assumptions  and  limitations  clear.  Where  possible,  experimental  data  on 
concentration  fluctuations  are  used  to  verify  these  models. 

Model  Refinements  to  Predict  Mean  Concentration 

The  existing  semi -empi ri cal  model  for  concentration  fluctuations  uses 
Gaussian  profiles  for  the  mean  concentration  and  the  fluctuation  variance. 

A steady  emission  rate  from  the  source  produces  a constant  mass  flux  of 
pollutant  at  each  downwind  plume  cross-section.  In  contrast,  the  model 
assumes  that  all  concentration  fluctuations  are  produced  close  to  the 
source  and  that  the  flux  of  concentration  variance  decays  rapidly  and 
continuously  in  the  downwind  direction.  As  a result,  the  effective  source 
strength  of  concentration  fluctuation  variance  is  a function  of  downwind 
distance,  while  the  mean  source  strength  remains  constant. 

The  fluctuating  concentration  theory  developed  by  Wilson  (1982)  was 
originally  optimized  to  produce  accurate  predictions  for  the  fluctuation 
intensity  and  intermi ttency  factor  within  a plume.  One  important  feature 
of  the  original  model  was  the  use  of  a pair  of  effective  variance  "sources" 
located  on  either  side  of  the  plume  centerline.  However,  when  the  model 
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was  applied  to  predict  conditional  mean  concentrations  it  was  found  that 
the  use  of  a pair  of  variance  sources  gave  unrealistically  shaped  profiles 
and  a poor  agreement  with  measured  wind  tunnel  data.  The  model  was 
modified  to  incorporate  only  a single  variance  source  located  at  the  same 
point  as  the  mean  emission.  The  use  of  a single  source  greatly  improved 
the  accuracy  of  prediction  for  the  conditional  mean  concentration  and  the 
intermittency  factor,  which  is  the  fraction  of  time  that  non-zero 
concentrations  occur.  With  this  small  change  in  source  geometry  the  model 
was  both  simpler  to  use  and  gave  better  estimates  for  the  conditional  mean 
and  fluctuation  statistics. 

One  major  problem  in  monitoring  atmospheric  pollution  levels  is  to 
interpret  the  highly  intermittent  concentration  records  that  occur  when  the 
wind  fails  to  carry  the  plume  directly  over  the  monitoring  station.  A 
possible  solution  is  to  set  air  quality  standards  on  the  basis  of  the 
conditional  mean  concentration  determined  by  discarding  the  time  intervals 
with  zero  concentration.  Because  probability  distributions  used  to 
estimate  the  fraction  of  time  concentration  exceeds  some  threshold  level 
require  the  conditional  mean  concentration  and  intermittency,  it  is 
sensible  to  express  air  quality  standards  in  terms  of  these  conditional 
stati sties. 

An  important  result  of  the  present  study  was  the  observation  that  this 
conditional  mean  concentration  will  take  on  a constant  value  rather  than 
decaying  to  zero  as  the  cross  wind  distance  off  the  centerline  increases. 
This  constant  value  for  the  conditional  mean  in  the  fringes  of  the  plume  is 
typically  1/5  to  1/2  of  the  centerline  conditional  concentration,  and 
depends  on  the  amount  of  meandering  that  occurs  at  the  plume  centerline. 
This  constant  value  for  the  conditional  mean  is  physically  realistic 
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because,  in  the  outer  edges  of  the  plume  the  dilution  process  can  be 
visualized  as  a series  of  random  steps  with  a contaminated  eddy  entraining 
uncontaminated  air.  The  dilution  of  these  eddies  should  be  a function  of 
travel  time,  and  not  how  far  the  puffs  are  off  the  centerline. 

From  the  practical  standpoint  of  monitoring  ambient  concentration,  the 
conditional  mean,  with  zero  concentration  periods  discarded,  appears  to  be 
a logical  choice  because  it  remains  large  in  the  fringes  of  the  plume  and 
so  gives  a better  indication  of  concentrations  on  the  plume  centerline  for 
a monitoring  station  located  off  the  centerline.  The  major  disadvantage  of 
using  conditional  mean  values  is  the  need  for  long  sampling  periods.  For 
non-zero  values  that  occur  only  10%  of  the  time,  a 10  hour  data  record  will 
be  required  to  extract  a one  hour  conditional  mean.  The  unsteady  nature  of 
atmospheric  wind  and  stability  conditions  may  make  it  very  difficult  to 
interpret  such  long  data  records. 

The  Effect  of  Sampling  Time  on  Concentration 

The  sampling  time,  also  called  the  exposure  time,  is  the  continuous 
time  interval  over  which  the  receptor  samples  the  atmosphere.  In  most 
engineering  and  laboratory  situations,  the  sampling  time  is  long  enough  to 
include  all  the  significant  turbulent  fluctuations,  and  steady  values  for 
the  mean  and  fluctuation  statistics  are  obtained.  In  the  atmosphere, 
sampling  times  are  typically  a few  minutes  to  a few  hours  in  length,  and  do 
not  sample  all  the  low  frequencies  in  the  turbulent  spectrum.  As 
atmospheric  sampling  time  increases,  both  the  mean  and  fluctuating 
statistics  will  vary,  even  though  the  mean  windspeed  and  atmospheric 
stability  remain  constant. 

Changing  this  sample  time  significantly  affects  the  observed  plume 
spread,  mean  concentrations  and  fluctuation  statistics.  The  influence  of 
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sample  time  is  due  to  the  large  scale  turbulence  and  unsteadiness  of 
atmospheric  conditions  which  cause  plume  meandering,  mainly  in  the  cross 
wind  direction.  As  sampling  time  increases,  more  of  these  low  frequency 
fluctuations  are  included  in  the  sample  and  the  mean  concentration  is 
reduced  as  crosswind  spread  increases.  Plume  fluctuation  statistics,  in 
particular  the  intermittency,  are  also  strongly  affected  by  sampling  time. 
As  the  sampling  time  increases  meandering  induced  by  large  scale  turbulence 
and  unsteadiness  exposes  a receptor  on  the  plume  centerline  to  the  highly 
intermittent  fluctuations  on  the  fringes  of  the  instantaneous  plume. 

The  theoretical  model  developed  by  Wilson  (1982)  to  predict 
concentration  fluctuations  was  optimized  using  a data  base  of  wind  tunnel 
observations.  The  use  of  wind  tunnel  data  to  adjust  the  model's  parameters 
limit  its  use  to  neutral  atmospheric  stability,  and  to  short  sample  times  a 
few  minutes  in  duration.  To  overcome  the  model's  limitation  to  short 
sample  times,  imposed  by  the  use  of  wind  tunnel  data,  a meandering  plume 
model  similar  to  that  suggested  by  Gifford  (1959)  was  used.  Long  sampling 
times  were  simulated  by  allowing  the  crosswind  profiles  of  the  short 
sampling  time  statistics  to  meander  back  and  forth  with  a Gaussian 
probability  distribution  over  a receptor  on  the  plume  centerline.  Using 
this  meandering  plume  model  it  was  found  that  the  centerline  becomes  more 
intermittent  and  the  mean  concentration  decreases  with  increasing  sample 
time,  as  expected.  In  contrast,  the  conditional  mean  concentration  with 
zero  readings  removed  is  virtually  independent  of  sampling  time,  which 
demonstrates  that  plume  meandering  mainly  affects  intermittency  and  not  the 
pollutant  concentration.  The  weak  dependence  of  conditional  mean 
concentration  on  sampling  time  lends  further  support  to  its  use  as  a 
measure  of  ambient  air  quality. 
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Because  the  concentration  fluctuation  theory  had  its  adjustable 
parameters  set  using  wind  tunnel  data,  there  is  some  uncertainty  in 
applying  it  to  situations  where  the  atmospheric  stability  is  not  neutral. 

To  help  overcome  the  theory's  limitation  to  neutral  atmospheric  stability, 
downwind  distance  was  expressed  in  terms  of  a normalized  plume  spread. 
Hopefully,  the  variation  of  plume  spread  with  changing  atmospheric 
conditions  would  automatically  correct  the  theory  for  stability  effects. 

The  present  study  carried  out  a detailed  examination  of  this  normalized 
spread  parameter  for  varying  atmospheric  stabilities,  and  concluded  that, 
at  least  for  mean  concentration,  this  normal i zation  greatly  reduces  the 
variability  caused  by  changing  atmospheric  stability. 

Influence  of  Data  Averaging  Time  and  Receptor  Response  Time 

Two  other  time  factors  complicate  the  interpretati on  of  atmospheric 
concentration  fluctuations.  The  first  of  these  is  the  receptor  response 
time,  which  is  a measure  of  the  ability  of  a receptor  to  respond  to  sudden 
changes  in  concentration.  Both  biological  receptors  and  measuring 
instruments  have  inertia  which  causes  their  response  to  lag  behind  sudden 
changes  in  input  conditions.  The  effect  of  this  response  time  constant  is 
to  cause  the  system  to  lose  its  sensitivity  to  high  frequency  concentration 
fluctuations,  (in  contrast  to  sampling  time  which  changes  the  contribution 
of  low  frequency  fluctuations). 

To  account  for  this  insensitivity  to  high  frequency  fluctuations 
methods  were  developed  in  the  present  study  to  correct  the  actual 
concentration  fluctuation  variance  to  an  effective  variance  with  high 
frequency  fluctuations  smoothed  out  by  the  averaging  time. 
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One  example  of  receptor  time  response  is  apparent  in  the  difference  in 
fluctuation  levels  observed  indoors  compared  to  their  outdoor 
counterparts.  Here,  the  mixing  process  within  a building  produces  a large 
time  constant  whose  slow  response  can  reduce  the  indoor  fluctuation 
intensity  to  10%  or  less  of  its  outdoor  level. 

To  account  for  this  sensitivity  to  high  frequency  fluctuations  methods 
were  developed  in  the  present  study  to  correct  the  actual  concentration 
fluctuation  variance  by  reducing  it  to  an  effective  variance  with  the  high 
frequencies  attenuated.  It  was  found,  as  expected,  that  the  receptor  time 
constant  prevented  short  intermittent  periods  of  zero  concentration  from 
being  felt  by  the  receptor  and  caused  the  effective  fluctuations  to  be  less 
intermittent  when  filtered  through  the  receptor's  time  response. 

The  second  time  factor  which  influences  concentration  fluctuations  is 
the  data  averaging  time.  In  a sense,  this  is  an  artificial  adjustment, 
because  it  is  simply  the  time  over  which  we  choose  to  average  a concen- 
tration record.  For  example,  many  ambient  concentration  records  are 
averaged  over  successive  one  hour  intervals.  This  process  of  hourly 
averaging  removes  the  high  frequency  components  from  the  concentration 
fluctuation  data.  To  compare  theoretical  predictions  from  a dispersion 
model  to  these  hourly  averaged  values  it  is  necessary  to  account  for  the 
averaging  process.  Methods  are  developed  for  correcting  the  variance  by 
assuming  a frequency  spectrum  function  and  then  attenuating  the  high 
frequency  components  for  a specified  averaging  time. 

Sampling  time  and  averaging  time  are  two  quite  different  parameters 
that  are  often  confused  with  one  another.  The  sampling  time  increases 
turbulent  concentration  fluctuations  by  exposing  a receptor  to  more  of  the 
low  frequency  components  in  the  meandering  plume.  Averaging  time  has 
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exactly  the  opposite  effect,  reducing  the  variance  as  high  frequency 
components  are  removed  as  the  data  is  smoothed  by  averaging  over  successive 
time  intervals.  Sampling  time  (also  called  exposure  time)  is  influenced  by 
the  natural  process  of  atmospheric  dispersion  while  averaging  time  is  set 
to  smooth  out  fluctuations  by  a conscious  choice  for  an  averaging  interval 
to  reduce  data  sets  to  manageable  size. 

Probability  Distributions  and  Threshold  Crossing  Probability 

For  pollutants  whose  effects  are  cummulative,  a receptor  responds  to 
the  total  dose  integrated  over  the  time  of  exposure.  At  the  other  extreme, 
the  response  time  to  a highly  toxic  gas  such  as  hydrogen  sulfide  is  rapid 
enough  that  the  damage  may  be  best  described  by  a threshold  concentration 
level  below  which  there  is  no  response  and  above  which  short  exposures 
cause  immediate  damage. 

Two  distinctly  different  probabilities  are  required  to  estimate  the 
risk  from  these  hazards.  The  first  is  a probability  distribution  function 
which  uses  the  conditional  mean,  conditional  variance,  and  i ntermi ttency  in 
the  plume  to  predict  the  fraction  of  time  that  instantaneous  concentrations 
or  an  integrated  dosage  will  exceed  a specified  limit.  In  addition,  for 
receptors  and  pollutants  that  produce  a rapid  response  it  is  also  essential 
to  know  the  probability  of  being  exposed  for  a specified  time  without  ever 
once  crossing  a critical  concentration  threshold.  Analytical  relationships 
are  developed  between  these  two  probabilities,  and  it  is  shown  that  simple 
closed  form  estimates  of  the  threshold  crossing  probability  can  be  made 
directly  from  the  fraction  of  time  exceedance  probability  distribution. 

Several  probability  distribution  functions  are  reviewed  as  candidates 
to  describe  the  variability  in  concentration  fluctuations.  The  Gamma 
distribution  for  fraction  of  time  exceedance  is  shown  to  be  superior  to  the 
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more  common  exponential  and  log-normal  distributions  used  by  previous 
investigators.  The  Gamma  distribution  is  tested  against  several 
independent  data  sets  from  the  atmosphere  and  wind  tunnel,  and  shown  to 
give  good  estimates  of  the  frequency  of  occurrence  of  concentrations  both 
higher  and  lower  than  the  mean. 

Frequency  of  occurrence  probability  distributions  are  skewed,  with  the 
most  probable  concentration  occurring  at  a value  somewhat  less  than  the 
conditional  mean  value.  Because  of  this  skewed  frequency  distribution,  the 
instantaneous  concentration  will  exceed  the  conditional  mean  only  30  to  40% 
of  the  time  during  which  non-zero  concentrations  occur.  This  result 
emphasizes  the  need  for  an  accurate  prediction  of  the  frequency  of 
occurrence  of  concentrations  less  than  the  conditional  mean,  and  supports 
the  conclusion  that  the  gamma  distribution  is  a better  choice  than  the 
log-normal  or  exponential  both  of  which  are  less  accurate  at  low 
concentrations. 

A probability  model  of  the  dilution  process  of  pollutant  released  from 
a point  source  is  developed.  The  dilution  is  modeled  as  a Wiener  process, 
which  is  a continuous  approximation  to  a random  walk  sequence  of 
dilutions.  The  results  of  the  model  show  that  in  an  initially  Gaussian 
distribution  function  is  rapidly  skewed  with  time  of  travel,  until 
meandering  produces  an  exponential  frequency  of  occurrence  distribution. 

The  Gamma  distribution,  which  reduces  to  an  exponential  distribution  when 
the  conditional  fluctuation  intensity  approaches  unity,  is  shown  to  be  a 
reasonable  approximation  to  this  downwind  distance  dependent  probability 
distribution  function. 

Threshold  crossing  probability  represents  a different  approach  to 
estimating  risk  of  exposure.  Because  there  is  no  data  available  from  which 
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to  predict  this  first  time  of  crossing  probability,  a simple  model  was 
developed  to  allow  their  calculation  directly  from  the  known  probability 
distribution  for  fraction  of  time  that  concentration  exceeds  this 
threshold.  By  assuming  that  concentrations  in  successive  time  intervals 
are  statistically  independent,  it  is  shown  that  the  probability  of 
exceeding  the  threshold  for  some  fixed  exposure  time  is  simply  the  negative 
exponential  of  the  product  of  the  length  of  exposure  time  and  the  fraction 
of  time  that  the  probability  exceeds  the  threshold.  This  simple  closed 
form  solution  allows  a rapid  calculation  of  first  threshold  crossing 
probabilities  from  the  fraction  of  time  exceedance  distributions. 
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CHAPTER  1 


INTRODUCTION 


This  report  presents  the  results  of  efforts  to  refine  and  extend  the 
concentration  fluctuation  model  of  Wilson  (1982)  to  predict  the  risk  of 
exposure  to  peak  concentrations  in  plumes.  The  model  is  based  on  wind 
tunnel  data,  which  we  will  show  simulates  exposure  times  of  a few 
minutes,  to  estimate  the  hazard  from  short  duration  transient  releases  of 
toxic  gas  from  pipeline  ruptures. 

Our  primary  objective  is  to  extend  this  model  to  the  much  longer 
sampling  times  of  one  hour  to  several  days  that  are  useful  in  evaluating 
dosage  and  peak  exposure  probabilities  for  continuous  plumes  typical  of 
industrial  sources.  To  carry  out  this  extrapolation  it  is  necessary  to 
develop  a new  theoretical  model  for  the  effect  of  long  term  plume  meandering 
on  mean  and  fluctuation  statistics.  The  concepts  of  sampling  time,  data 
averaging  time,  and  receptor  response  time  must  be  included  so  that  raw 
fluctuation  statistics  in  the  plume  can  be  adjusted  to  determine  the  response 
of  animals  and  vegetation. 

To  make  the  existing  model  useful  for  long  term  exposures,  a whole  new 
set  of  analytical  tools  were  required.  This  report  reflects  the  development 
of  those  tools  by  the  detailed  appendices  which  present  their  derivation. 

Each  of  these  appendices  deals  only  with  a specific  topic,  and  in  the  report 
we  will  bring  them  together  to  show  how  they  interact,  and  how  they  can  be  used 
to  develop  realistic  predictions. 

As  the  model  becomes  more  complex,  one  philosophical  question  is 
whether  concentration  fluctuation  predictions  are  really  necessary  for 
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continuous  plumes.  Barry  (1977)  points  out  the  dramatic  increase  in  the 
catastrophic  extreme  value  episodes  that  will  be  caused  by  small  changes 
in  the  long  term  mean  value  of  concentration.  It  is  only  by  understanding 
the  statistics  that  lead  to  these  extremes  that  we  can  effectively  assess 
the  cost-benefit  ratio  for  pollution  control  measures  to  lower  this 
mean  value.  The  need  for  concentration  fluctuation  probabilities  is  also 
evident  when  we  consider  the  mean  values  only  predict  the  concentration 
that  will  be  exceeded  50%  of  the  time.  However,  it  is  the  extreme  values 
that  occur  only  a few  percent  of  the  time  that  may  be  responsible  for  most 
problems  related  to  health  and  vegetation  damage.  We  should  determine 
not  only  a range  of  acceptable  concentration  levels,  but  also  specify 
the  frequency  with  which  these  levels  may  be  exceeded. 
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CHAPTER  2 


REFINEMENT  OF  THE  CONCENTRATION  FLUCTUATION  MODEL 


In  the  preceding  report  Wilson  (1982)  developed  a model  to  predict 
concentration  fluctuation  statistics  across  a plume  profile.  The  adjustable 
constants  in  this  semi -empirical  model  were  set  by  matching  wind  tunnel 
measurements  of  Fackrell  and  Robins  (1982)  who  used  tracer  gases  in  a 
neutrally  stable  wind  tunnel  boundary  layer  to  simulate  atmospheric  rel eases 
from  ground  level  and  elevated  sources.  The  major  features  of  this 
modified  Gaussian  dispersion  model  are: 

• The  use  of  variance  "source"  with  a strength  proportional  to  the 
mean  concentration  source  and  to  the  ratio  of  source 
size  to  turbulent  scale.  Unlike  the  mean  source  which 
remains  constant,  the  variance  source  first  increases 
and  then  dissipates  with  downwind  distance. 

t The  use  of  an  image  sink  term  for  variance  to  simulate  the 
dissipation  of  concentration  fluctuations  in  the  highly 
sheared  flow  near  the  ground  surface. 

t The  use  of  a pair  of  variance  sources  displaced  on  either 
side  of  the  plume  axis  and  located  above  the  actual  point 
of  release  to  account  for  distributed  production  of  concentra- 
tion fluctuations  in  the  central  region  of  the  plume. 

The  details  are  discussed  in  two  papers  by  Wilson,  Fackrell  and  Robins 
(1982).  It  was  found  later  by  Wilson  (1982)  that  the  theory  could  be  extended 
to  predict  conditional  plume  statistics  such  as  intermittency , conditional 
intensity  and  conditional  mean  without  changing  any  of  the  adjustable  constants 
used  in  the  original  model.  This  suggests  that  the  model  is  not  simply  an 
exercise  in  curve  fitting  with  adjustable  constants,  but  has  a physical  basis. 

One  important  characteristic  of  the  model  is  that  concentration 
fluctuations  do  not  decrease  to  zero  at  the  ground  surface  as  might  be 
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expected  from  physical  reasoning.  Instead,  the  wind  tunnel  data  on 
which  the  model  is  based  indicates  that  the  fluctuation  variance  at  ground 
level  is  about  half  the  magnitude  of  the  variance  that  would  be  expected 
with  no  dissipating  surface  present.  This  result  has  been  confirmed  by 
independent  measurements  of  Bara  (1985)  in  a water  channel  simulation  of 
a ground  level  release.  Bara's  measurements  extended  close  enough  to  the 
surface  to  be  within  the  roughness  elements  themselves,  at  a height  of 
about  twice  the  log-law  roughness  length  Zq  , and  demonstrated  the  existence 
of  large  concentration  fluctuations  in  this  region.  Why  do  these  fluctua- 
tions persist  at  the  ground  surface  where  high  rates  of  dissipation  should 
have  removed  all  concentration  variance?  The  reason  appears  to  lie  in  the 
crosswind  meandering  of  the  plume,  which  exposes  a receptor  to  concentration 
fluctuations  caused  by  the  crosswind  variation  of  mean  concentration.  As 
this  crosswind  c profile  meanders  back  and  forth  over  a receptor  it  generates 

the  low  frequency  fluctuations  observed  near  the  ground  surface,  and  causes  a 
~r2 

variance  c to  appear. 

Crosswind  Spacing  of  the  Variance  Sources 

The  model  uses  a pair  of  variancesources  displaced  by  some  fraction 
cp  of  the  crosswind  plume  spread  on  either  side  of  the  plume  axis  (for 
plumes  normalized  by  their  half  width  6 = 1.77  a the  displacement  is  a 

y y 

fraction  3 = 0.85<f>).  Figure  1 shows  that  a displacement  3 = 0.6  (or<})=  0.706) 
gives  a better  prediction  of  the  crosswind  profile  than  a single  variance 
source  with  3 = 0 (or  <j>  = 0). 

Closer  examination  of  the  overall  performance  of  the  model  in  predicting 
conditional  plume  statistics  indicates  that  a single  variance  source,  first  sug 
gested  by  Csanady  (1973)  is  a better  choice.  This  is  apparent  in  Figure  2 
which  shows  an 
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Figure  1 Crosswind  variance  profiles  from  an  elevated  source  (normalized 
by  extrapolated  centerline  variance  c^)  from  Wilson,  Fackrell 
and  Robins  (1982). 
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Figure  2 Effect  of  source  pair  spacing  3 on  crosswind  intermittency 
profiles  from  a ground  level  source  in  a wind  tunnel.  Data 
of  Robins  and  Fackrell  (1979). 
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improved  agreement  with  crosswind  intermittency  when  a single  variance 

source  on  the  plume  centreline  is  substitued  for  the  source  pair  with 

spacing  3 = 0.6.  The  three  most  important  parameters  predicted  by  the 

model  are  the  conditional  "plume"  mean  c , plume  intensity  i and  inter- 

P P 

mittency  y . These  are  the  required  input  variables  needed  for  the  probabil- 
ity distribution  used  to  predict  the  risk  of  exceeding  some  specified  peak 
level.  The  choice  of  variance  source  spacing  has  no  effect  on  i , which  is 
assumed  to  be  constant  across  the  plume. 

Figure  3 shows  that  the  crosswind  variation  of  the  conditional  plume 
mean  c^  is  better  estimated  using  the  single  variance  source  with  $ = 0. 

In  Figure  4 this  single  variance  source  continues  to  give  a good  fit  to  the 
data  at  a location  further  downwind,  close  to  the  position  of  ground  level 
maximum  concentration.  When  comparing  the  model  with  measurements  in  the 
figures,  it  should  be  kept  in  mind  that  with  the  exception  of  mean  concentration, 
none  of  the  predicted  profiles  are  forced  to  agree  with  experimental  data  by 
normalizing  with  centreline  values.  Differences  between  predicted  curves  and 
measurements  are  due  to  both  data  scatter  and  the  deficiencies  in  predicting 
variance  source  strengths  in  the  plume.  Considering  this,  the  agreement 
between  experiment  and  theory  is  remarkably  good. 

Because  most  receptors  are  on  the  ground,  it  is  less  important  to 
predict  accurately  vertical  plume  profiles.  On  the  other  hand,  predicting 
vertical  profiles  provides  the  most  stringent  test  of  the  model's  ability  to 
simulate  dissipation  effects  on  concentration  fluctuations  near  ground  level. 
Vertical  profiles  of  Fackrell  and  Robins'  wind  tunnel  dataare  shown  in  Figures 
5,  6 and  7 at  three  downwind  distances  from  an  elevated  source.  For  the  source 
height  of  h = 0.192  H the  maximum  mean  ground  level  concentration  occurs  at 


Crosswind  Distance 


Figure  3 Crosswind  profiles  of  mean  and  fluctuation  statistics  at 

x/H  = 0.83  in  the  plume  from  an  elevated  source  d = 0.30  cm 
in  a neutrally  stable  wind  tunnel  boundary  layer,  H = 120  cm, 
z /H  = 4 x 10~4.  Data  of  Fackrell  and  Robins  (1982a). 
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Figure  4 Crosswind  profiles  of  mean  and  fluctuation  statistics  at  x/H  = 

2.88  in  a plume  from  an  elevated  source  d = 0.85  cm  at  h = 23  cm 
in  a neutrally  stable  wind  tunnel  boundary  layer  H = 120  cm, 
zQ  = 4 x 10-4.  Data  of  Fackrell  and  Robins  (1982a). 
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Figure  6 Vertical  profiles  of  mean  and  fluctuation  statistics  at  x/H  = 2.88 
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Figure  7 Vertical  profiles  of  mean  and  fluctuation  statistics  at  x/H 
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about  x/H  - 3.0,  with  measured  profiles  shown  at  x/H  = 0.83,  2.88  and 
4.76.  These  vertical  profiles  are  independent  of  the  crosswind  source  pair 
spacing  0,  but  do  show  that  a vertical  variance  source  displacement 
<j>z  = 0.706  is  still  better  for  predicting  variance  and  i ntermi ttency  in  the 
upper  region  of  the  plume  above  source  height.  Because  some  vertical  dis- 
placement of  the  variance  source  is  essential  in  predicting  the  vertical 
profiles  for  a ground  level  source,  we  will  continue  to  use  a value  of 
<j)z  = 0.706  (or  3Z  = 0.6). 

In  summary,  a review  of  the  data  indicates  that  a single  variance 

source  with  0=0,  but  with  vertical  displacement  0z  = 0.706  is  the  best 

~r2 

choice.  Although  a single  variance  source  gives  a poor  fit  to  c in  the 
central  core  of  the  plume,  it  gives  significantly  better  agreement  for  variance, 

r2  _ 
c , intermittency  y and  plume  mean  concentration  c^  in  the  fringes  of  the  plume 

where  y > 2 a . This  does  not  imply  that  the  true  profile  of  variance  is 

Gaussian  in  the  outer  regions,  because  our  assumption  of  constant  plume 

intensity  i^  in  the  outer  regions  may  not  hold,  and  the  use  of  a single  variance 

source  may  simply  provide  a compensating  error. 

Finally,  the  most  important  benefit  of  using  a single  variance  source 
is  that  it  prevents  unrealistic  increases  in  the  probability  of  exceeding  a 
threshold  level  in  the  outer  crosswind  regions.  These  oscillations  in  exceed- 
ence probability  predicted  by  Wilson  (1982),  completely  disappear  when  a 
single  variance  source  is  used.  From  this,  we  may  conclude  that  simple  models 
are  often  best,  and  that  correction  factors  which  improve  their  performance 
in  limited  regions  often  do  more  harm  than  good  to  the  overall  accuracy. 
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Accounting  for  Atmospheric  Stability 

The  concentration  fluctuation  model  expresses  the  downwind  variation  of 
source  strength  in  terms  of  a dimensionless  plume  spread  parameter  X defined 
in  terms  of  the  vertical  and  crosswind  plume  spreads  a and  az 


where  H is  the  thickness  of  the  atmospheric  boundary  layer,  and  may  be  taken 
as  H - 650  m.  This  boundary  layer  thickness  is  simply  a means  of  including 
the  scale  of  atmospheric  turbulence,  and  Wilson,  Fackrell  and  Robins  (1982) 
recommend  for  neutrally  stable  wind  tunnel  boundary  layers 


in  terms  of  the  neutral  log  law  surface  roughness  length  ZQ. 

It  is  assumed  in  the  model  that  X accounts  for  all  effects  of  turbulence 
intensity  and  scale  caused  by  changes  in  stability  and  surface  roughness.  To 
examine  this  assumption,  the  equations  for  mean  ground  level  concentration 
along  the  plume  axis  are  transformed  in  Appendix  B to  a coordinate  system 
based  on  the  spread  parameter  X . The  influence  of  stability  was  examined  by 
considering  the  extremes  of  "very  unstable"  (Class  A)  and  "very  stable" 

(Class  F)  with  plume  spread  parameters  recommended  by  Angle  (1978).  Normalized 
ground  level  concentration  profiles  shown  in  Figure  8 confirm  that  the  effects 
of  stability  are  negligible  when  downwind  distance  is  expressed  in  terms  of  the 
plume  spread  parameter  X rather  than  downwind  distance  x.  In  plume  spread 
coordinates  the  concentration  profiles  have  a maximum  variation  of  about  5% 


(1) 


H * 2500  Z 
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Figure  8 Ground  level  profiles  of  mearr  concentration  normalized  in 
downwind  distance  x and  plume  spread  X coordinates  using 
plume  spreads  a and  az  from  Angle  (1978). 


TABLE  1 


Variability  of  Location  of  Maximum  Ground  Level 
Concentration  with  Atmospheric  Stability  for  Two  Source  Heights 


Stability 

Class 

Source  Height  h = 10  m 

Source  Height  h = 65  m 

xm 

A 

m 

X 

m 

A 

m 

X 

3 

O 

A n 
m,D 

o 

E 

X 

X n 
m,D 

A 

0.37 

1.16 

0.23 

0.94 

D 

1.00 

1.00 

1.00 

1.00 

F 

3.58 

1.20 

5.76 

1.48 
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from  the  neutral  stability  (Class  D)  conditions  on  which  our  model  is  based. 

Another  indication  of  how  well  the  plume  spread  parameter  A performs 

as  a universal  distance  scale  is  the  location  of  maximum  mean  concentration 

A^.  Table  1 compares  the  variability  of  the  distance  to  ground  level  maximum 

in  both  A and  x coordinates.  It  is  clear  from  this  that  A is  much  closer  to 

m 

a universal  constant,  varying  by  only  about  + 20%  while  the  distance  to 
maximum  concentration  x varies  by  more  than  + 400%  over  the  same  range  of 
stabilities.  From  this  we  expect  that  the  concentration  fluctuation  model 
might  give  reasonable  predictions  for  all  atmospheric  stability  classes  when 
downwind  distances  are  expressed  in  the  form  of  plume  spread  parameter  A . 

Limitations  Imposed  by  Wind  Tunnel  Simulation 

Wind  tunnel  simulation  provides  the  basis  for  setting  adjustable  constants 
in  the  concentration  fluctuation  model  developed  by  Wilson  (1982).  The  working 
equations  for  this  model  are  presented  in  Appendix  B for  a receptor  located 
at  ground  level.  While  a wind  tunnel  boundary  layer  is  essential  to  provide 
the  reproducible  conditions  needed  to  measure  fluctuation  statistics  accurately, 
this  small  scale  simulation  does  not  adequately  represent  atmospheric  eddies 
that  cause  the  plume  meandering  over  periods  of  a few  minutes  to  a few  hours. 
Also,  the  wind  tunnel  is  incapable  of  simulating  the  long  term  meandering 
of  wind  direction  which  occurs  over  sample  time  periods  of  a few  hours  to 
several  days.  The  mean  concentration  on  the  centreline  decreases  continuously 
as  meandering  causes  the  short  term  "instantaneous"  plume  to  move  back  and 
forth  over  the  receptor.  This  nonstationary  decrease  in  mean  concentration 
is  shown  schematically  in  Figure  9,  where  the  limitation  of  wind  tunnel 
simulation  is  clearly  evident. 


Cj  Mean  Concentration  for  SampleTime 


Figure  9 Sampling  time  limitations  in  wind  tunnel  simulation 
of  atmospheric  dispersion 
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The  wind  tunnel  data  of  Fackrell  and  Robins  are  for  a boundary  layer 
thickness  of  H = 1.2  m to  simulate  the  atmospheric  boundary  layer  in  about 
550:1  scale.  At  a scale  downwind  distance  of  about  1.5  km  the  wind  tunnel 
model  produced  a full  scale  equivalent  crosswind  spread  a = 1 28  m.  The 
rough  surface  of  the  wind  tunnel  boundary  layer  generated  a mean  wind  profile 
with  an  equivalent  full  scale  log  law  roughness  height  of  Zq  = 26  cm. 

Using  the  power  law  crosswind  spread  equations  recommended  by  Smith  from 
Pasquill  (1974)  p.  375,  the  predicted  crosswind  spread  for  an  unspecified  sample 
time  is  about  a = 120  m.  Pasquill  suggests  that  these  crosswind  spread 
values  in  the  atmosphere  are  typical  of  3 minute  sample  times,  and  the  close 
agreement  between  the  measured  and  predicted  values  suggest  that  the  wind  tunnel 
simulation  is  limited  to  turbulence  scales  that  are  typical  of  three  minute 
samples  in  the  atmosphere.  These  short  exposure  times  are  exactly  what  is 
needed  to  simulate  the  passage  of  a transient  plume  of  toxic  gas  from  a pipe- 
line rupture.  On  the  other  hand,  long  term  exposures  of  several  hours  to 
several  days  that  are  required  for  determining  damage  to  plants  and  animals 
from  continuous  releases  cannot  be  predicted  from  the  existing  model.  In 
Chapter  3 a procedure  will  be  developed  for  extrapolating  the  existing  three 
minute  sampling  time  predictions  to  estimate  concentration  statistics  for 
longer  exposure  times. 

Wind  tunnel  simulation  can  create  the  illusion  of  predicting  concentra- 
tions for  sampling  times  as  long  as  one  hour.  If  only  mean  concentration  is 
simulated,  the  loss  of  large  scale  eddies  which  cause  plume  meandering  can  be 
compensated  for  by  adding  an  excess  of  turbulence  to  the  wind  tunnel  boundary 
layer.  The  mean  concentration  depends  only  on  the  vertical  and  crosswind  plume 
spreads  oz  and  a , and  the  vertical  spread  remains  relatively  constant  for 
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sampling  times  longer  than  about  three  minutes  in  the  atmosphere.  The 


crosswind  spread  o j continues  to  increase  with  sampling  Ts  as  the  effective 
scale  of  the  velocity  fluctuations  becomes  larger.  A simple  power  law  is 


a 

yT 


v y. 


'2 


(B  A T)1_ay  ay 

o vT  J x 


(3) 


where  V v /U  is  the  intensity  and  Bq  = Lv.v'/AvjU  is  the  Lagrangian  scaling 
constant,  and  Avy  is  the  integral  scale  of  crosswind  turbulence 
for  a sample  time  duration  Tg.  From  this  it  is  easy  to  see  that  an  ingeneous 
wind  tunnel  modeller  can  produce  the  required  crosswind  spread  for  any 
sample  time  simply  by  generating  more  upstream  turbulence  v to  compensate 


for  a value  of  A^j  that  is  too  small.  Typically  a - 0.88  so  that 
0 1 2 

~ A vy  . This  weak  dependence  requires  only  a 10%  change  in  tu 
intensity  to  account  for  a factor  of  two  error  in  simulating  scale. 


However,  if  a wind  tunnel  is  used  to  model  concentration  fluctuation 
statistics,  it  can  never  properly  simulate  the  high  intermi ttency  caused  by 
meandering  over  long  sample  times.  In  the  atmosphere,  crosswind  meandering 
exposes  a receptor  to  the  highly  intermittent  fluctuations  in  the  outer  fringes 
in  the  crosswind  plume  profile.  Adjustments  to  wind  tunnel  turbulence  levels 
will  only  produce  more  spread,  and  not  the  intermittency  needed  to  simulate 
sample  periods  of  several  hours. 

Comparison  With  Full  Scale  Atmospheric  Data 

The  most  important  aspect  of  the  concentration  fluctuation  model  is  its 
ability  to  predict  the  crosswind  variation  of  plume  intensity  and  intermittency. 
The  field  measurements  of  Ramsdell  and  Hinds  (1971)  used  a radioactive  tracer 
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gas  to  study  crosswind  profiles  from  a ground  level  source.  Thirty-nine 
samplers  were  located  at  200  m and  800  m downwind,  and  averaged  over 
38.4  seconds  for  sampling  times  of  10  to  20  minutes.  Normalized  fluctuation 
statistics  were  found  to  be  similar  at  both  of  the  downwind  locations,  and 
the  data  was  combined  to  produce  a single  normalized  data  set. 

Their  crosswind  profiles  of  total  variance  in  Figure  10  show  the 
same  trend  observed  in  the  wind  tunnel  measurements  in  Figure  1.  In  both 
cases,  the  use  of  a pair  of  variance  sources  gives  a better  prediction  in  the 
core  of  the  plume  for  y < 2 o^.  In  the  fringes  beyond  two  standard  deviations 
a single  variance  source  (source  pair  spacing  <j>  = 0 ) gives  a better  fit.  The 
conditional  plume  statistics  of  intermittency  y,  plume  intensity  i and 
total  plume  intensity  i in  Figure  11  show  good  agreement  with  theoretical 
predictions  normalized  by  experimental  values  on  the  plume  centreline.  Here, 
the  use  of  a single  variance  source  gives  significantly  better  estimates 
of  total  intensity  and  intermittency  for  the  crosswind  profiles.  It  is 
gratifying  to  find  that  the  field  data  in  Figure  11  showed  an  approximately 
constant  value  for  conditional  plume  intensity  i^.  This  assumption  of  constant 
intensity  across  the  plume  is  the  most  important  element  used  in  constructing 
the  plume  fluctuation  model. 

Although  the  shape  of  crosswind  profiles  for  plume  statistics  is 
reasonably  well  simulated,  the  model  underestimates  the  magnitude  of  both 
ip  and  i.  A comparison  of  centreline  values  is  shown  in  Table  2.  Two 
probable  reasons  for  the  model's  failure  are: 

t The  predictions  are  limited  to  a short  sampling  time  of  about 
three  minutes,  while  the  field  data  was  typical  of  fifteen 
minute  sampling  times.  These  longer  sampling  times  allowed 
more  plume  meandering  and  accounted  for  the  lower  value  of 
intermittency  factor  y in  the  measurements. 


Total  Fluctuation  Variance 


23 


Figure  11 


Field  measurements  of  crosswind  profiles  of  total  and  plume 
intensities  and  intermi ttency . Data  of  Ramsdell  and 
Hinds  (1071). 
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TABLE  2 

PLUME  CENTERLINE  CONCENTRATION  FLUCTUATIONS  FROM 
RAMSDELL  AND  HINDS  (1971)  GROUND  LEVEL  SOURCE  FIELD  MEASUREMENTS 


FIELD  MEASUREMENTS 

MODEL  PREDICTION* 

VARIABLE 

Combined  Data  for 
x = 200  m and  800  m 

Sample  Time 
Ts  = 10  to  20  min 

Sample  Time 
Ts  = 3 min 

Intermittency 

y = 0.64 
To 

y = 1.00 
0 

PI  ume 
Intensity 

ipl  = 0.95 

i = 0.38 
P 

Total 

Intensity 

—J. 

—1 

O 

II 

• 

o 

i = 0.38 
0 

r 

*Using  average  A = 0.063  of  the  two  locations,  x = 200  m (A  = 0.031) 
and  x = 800  m ( A = 0.095)  for  d = 0.02  m,  h = 1.0  m az  = 35  m. 
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• The  model  predicts  no  effect  of  source  size  for  a ground 
level  release.  The  wind  tunnel  simulation  on  which  this 
was  based  considered  full  scale  equivalent  sources  of 
1.5  m to  20  m.  In  the  small  scale  turbulence  near  ground 
level  the  very  small  source  size  of  d ~ 0.02  m in  the  field 
measurements  may  have  significantly  affected  plume 
meandering. 

Neither  of  these  are  serious  deficiencies,  but  some  adjustment  to  source 
size  effects  in  the  model  are  necessary,  at  least  for  ground  level  sources. 
The  effect  of  sampling  time  will  be  accounted  for  by  a crosswind  meandering 
correction  presented  in  the  next  chapter. 
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CHAPTER  3 

SAMPLING,  AVERAGING  AND  RESPONSE  TIME  EFFECTS 

The  theory  (Appendix  B)  for  mean  and  fluctuating  concentration 
assumes  that  the  wind  blows  directly  from  the  source  to  the  receptor,  and 
that  fluctuation  statistics  are  measured  over  a long  enough  sampling  time 
to  include  all  relevant  large  scale  meandering  effects.  These  fluctuations 
are  then  assumed  to  act  on  a receptor  with  an  instantaneous  response  time 
which  follows  the  fluctuation  peaks.  In  practice,  neither  of  these  ideal- 
ized limits  is  ever  realized.  Barry  (1977)  suggests  that  a sample  time  of 
about  four  days  in  duration  would  be  required  to  include  all  fluctuations 
that  contribute  to  plume  dispersion.  Current  estimates  of  plume  spread 
parameters  a and  az  apply  only  for  short  sample  times  of  about  three  minutes. 
While  these  short  sample  times  are  adequate  for  predicting  concentration 
fluctuations  in  a brief  release  of  toxic  gas,  they  are  inappropriate  for 
long  term  estimates  of  dosage  and  peak  concentration  levels  over  several 
hours . 

Accounting  for  the  response  time  of  a receptor  is  equally  important 
in  estimating  the  effect  of  fluctuating  concentration.  The  response  time 
of  humans  and  animals  to  highly  toxic  gases  lies  in  the  range  of  a few 
seconds  to  a few  minutes.  At  the  other  extreme,  the  response  of  slow 
growing  vegetation  such  as  trees  to  time  varying  concentration  will  have  a 
time  response  measured  in  hours  or  in  days.  The  damage  to  such  a slowly 
responding  system  is  often  expressed  in  terms  of  the  total  integrated  dose 
received  over  the  exposure  time.  When  the  fluctuating  concentration  is 
averaged  over  these  exposure  time  intervals  high  frequency  fluctuations  are 
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removed  and  dosage  level  is  much  less  variable  than  the  instantaneous 
concentration.  The  parameters  required  to  account  for  these  time  factors 
are: 

T - The  sampling  time  duration  over  which  a receptor  is 
exposed  to  a plume. 

T - The  averaging  time  over  which  integrated  dosage  is 
a computed.  Ta  will  always  be  less  than  or  equal  to 

the  sampling  time  Ts. 

t$  - The  system  response  time  constant  of  the  receptor. 

3 - The  integral  time  scale  of  the  concentration  fluctuation 

c spectrum. 

Integral  Time  Scale  <7C  of  Concentration  Fluctuation 

In  order  to  estimate  the  effect  of  sampling  averaging  and  response 

~T2 

times  we  must  specify  how  the  spectrum  of  c varies  with  frequency.  In 
Appendix  A the  Markov  spectrum  is  used  as  the  simplest  reasonable  approxi- 
mation to  the  concentration  spectrum.  This  spectrum  is  consistent  with  an 
exponential  "inertialess"  autocorrel ation  function,  and  its  scale-^7,  is 

the  only  variable  required  to  specify  the  relative  contribution  of  each 

_2 

frequency  to  the  concentration  variance  c . 

The  wind  tunnel  data  of  Fackrell  and  Robins  is  used  in  Appendix  C to 
estimate  the  variation  of  this  integral  scale  with  position  in  the  plume 
and  downwind  distance.  Their  data  clearly  shows  a dramatic  increase  in  the 
size  of  concentration  scale  near  the  ground  where  the  highly  sheared 
atmospheric  turbulence  dissipates  small  (high  frequency)  fluctuations 
leaving  only  large  scale  turbulence  to  contribute  to  the  variance.  The 
empirical  model  is  used  to  relate  the  increase  in  integral  scale<Jc  to 
the  decrease  in  variance  near  the  ground  caused  by  the  same  dissipation 
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effects.  The  data  indicates  that  in  the  absence  of  ground  dissipation 
effects,  the  scale  of  concentration  fluctuations  remains  constant  across 
a plume  cross-section,  and  gradually  increases  with  downwind  distance  as 
dissipation  removes  small  scale  concentration  eddies.  Because  it  is  the 
size  of  eddies,  not  their  time  of  passage,  that  remains  constant  with 
changing  windspeed,  the  empirical  relations  in  Appendix  C are  presented  in 
terms  of  the  integral  length  scale  = <{7cU,  rather  than  the  time  scale. 

The  eddy  convection  velocity  is  assumed  to  the  local  mean  velocity  at  the 
height  z at  which  the  spectrum  is  measured. 

The  empirical  correlation  for  concentration  fluctuation  length  and 
time  scales  given  in  Appendix  C is  only  valid  for  the  sampling  time  T$  of 
about  three  minutes,  for  which  the  wind  tunnel  data  is  appropriate.  When 
a receptor  is  exposed  for  longer  sampling  times  the  slow  crosswind  meandering 
of  the  plume  will  contribute  more  large  scale  low  frequency  turbulence  to  the 
variance  and  cause  the  time  scale<-7c  to  increase.  The  correction  factors 
necessary  to  adjust  the  time  scale  for  sampling  times  longer  than  three 
minutes  are  developed  in  Appendix  D. 

Adjustment  for  Sampling  Time  - T$ 

For  a receptor  exposed  to  a continuous  plume  we  are  usually  interested 
in  the  response  over  a relatively  long  exposure  time  ranging  from  several 
hours  to  several  days.  Unfortunately,  there  is  no  simple  method  available 
for  extrapolating  the  mean  and  fluctuation  statistics  predicted  by  the  model 
for  an  ensemble  of  three  minute  sample  times  to  these  longer  exposures. 

To  allow  this  extrapolation,  a meandering  plume  model  is  developed  in 
Appendix  D to  simulate  the  effect  of  long  sampling  times.  The  position 
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of  the  plume  caused  by  meandering  is  assumed  to  have  a normal  probability 

distribution  with  a crosswind  standard  deviation  a which  increases  with 

m 

sample  time,  the  same  as  Gifford's  (1959)  meandering  plume  model. 


This  crosswind  meandering  process  is  shown  schematically  in  Figure 
12,  which  shows  by  the  symmetric  meander  probability  about  y = 0 that 
the  long  term  average  wind  direction  is  assumed  to  be  constant  during  the 
sample  period.  Because  we  are  interested  in  an  ensemble  average  of  many 
realizations  of  sample  periods  of  length  T$ , this  constant  wind  direction 
assumption  is  not  a severe  limitation  for  sample  times  up  to  several  hours. 
For  longer  sample  durations  the  effects  of  topography,  drainage  flows  and 
synoptic  trends  should  be  included  by  using  wind-rose  weighting  of 
directions.  By  using  the  suggested  empirical  relation  between  sample  time 
Ts  and  meander  width  given  in  (D44),  the  short  term  statistics  for 
sampling  time  T^^  are  related  to  longer  sampling  times  as  follows: 

From  (D42)  the  mean  concentration  cqj 


c 

o 


sref 


(4) 


where  p is  an  empirically  determined  exponent.  From  (D32)  and  (D44)  the 

total  fluctuation  intensity  i‘  is 
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Figure  12  Crosswind  meandering  model  for  sampling  time  adjustments  to 
plume  centerline  statistics. 
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For  long  sampling  times  where  Y /T  e.p  > 10  this  may  be  approximated  by 


i2T  + 1 = (0.707  + i 2 
oT  v o 


rJ 


(6) 


sref^ 


Meandering  increases  the  integral  scale^cT  of  concentration  fluctuations 
by  the  ratio  of  intensities.  From  (D41  ) 


coT 


.2 
Vr 
. 2 


(7) 


The  crosswind  avera 


ged  intermittency  y j must  be  found  by  numerical 


integration  of  (D19)  in  terms  of  am/a  , and  t^ien  re^ated  to  sample  time 
with  (D44).  The  conditional  "plume"  intensity  i ^ for  sampling  time  Ts  is 
calculated  from  the  values  of  y-j-  and  intensity  iQ^  and  the  relation 

1.2 

yoT  “ 


i t + 1 

Pi 


i2  + 1 
oT 


(8) 


The  intermittency  y in  the  3 minute  plume,  and  the  crosswind  average 
yQy  are  shown  in  Figure  13.  When  normalized  by  yQ,  the  3 minute  centerline 
value,  these  intermittencies  are  functions  of  the  centreline  intensity  iQ. 

For  sampling  times  long  enough  to  make  a m > 2a^,  the  plume  centerline 
intermittency  yQy  falls  to  about  half  its  three  minute  value  yQ.  Using  the 
recommended  value  of  p = 0.2  for  the  exponent  Figure  14  shows  the  predicted 
change  in  plume  centerline  statistics.  We  see  the  conditional  mean  value 
CpQ  levels  off  at  a constant  value,  while  the  mean  concentration  cq^.  and 
intermittency  y^  decrease  continuously  with  exposure  time. 
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Figure  13  Gaussian  model  prediction  of  intermi ttency  y across  a plume,  and  for  sampling  time  meander 
averaging  on  centerline  y0j,  with  single  variance  source  <J>  = 0. 
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Figure  14  Meandering  plume  model  prediction  of  sampling  time  effects  on  plume  centerline  statistics 
for  a single  variance  source  (<J>  = 0)  and  p = 0.2  using  equations  (D13),  (D19)  and  (D44). 
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This  sampling  time  adjustment  might  account  for  the  observed 
differences  between  Ramsdell  and  Hinds  field  measurements  and  the  model 
predictions  for  three  minute  plumes.  Table  3 shows  the  predicted  effect 
of  sampling  time  on  intermittency  and  fluctuation  intensities  for  the 
model  and  ground  level  field  measurements.  The  values  of  p = 0.2 
and  0.3  are  reasonable  estimates  of  plume  meander  exponent.  While  the 
correction  significantly  improves  agreement  between  the  model  and  measure- 
ments, the  theory  still  grossly  underpredicts  conditional  and  total 
fluctuation  intensity.  In  turn,  this  underestimates  the  intermittency. 

If  i = 0.95  and  y =1.0  are  used  for  the  three  minute  plume,  a 15  minute 
p 'o 

sample  time  has  iQ^  = 1.32  and  y^  = 0.71  for  p = 0.3.  These  are  in  much 
better  agreement  with  the  measurements,  and  suggest  that  the  empirical 
constants  in  the  three  minute  plume  model  be  adjusted  to  produce  larger 
ground  level  intensities  and  a more  intermittent  plume. 


Effect  of  Averaging  Time  - Tq 

The  averaging  time  T is  the  time  interval  chosen  by  an  investigator, 
a 

over  which  the  instantaneous  concentration  is  averaged  to  produce  a single 

value  for  each  interval.  Because  this  averaging  process  removes  high 

frequency  fl uctuations,  the  variance  of  the  averaged  values  will  be  less  than 

the  variance  of  the  instantaneous  fluctuations.  The  attenuation  factor  for 

the  frequency  spectrum  is  discussed  by  Pasquill  (1974)  pp.  12-16,  and  is 

applied  to  the  Markov  spectrum  in  Appendix  A.  For  the  most  common  case 

where  the  averaging  time  T is  much  longer  than  the  intergral  time  scale 

a 

<7C  of  concentration  fluctuations,  equation  (A31 ) is  a good  estimate  of 

~TZ 

reduced  variance  c for  averaging  time  T . 
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TABLE  3 

PLUME  CENTERLINE  STATISTICS  FOR  RAMSDELL  AND  HINDS  (1971) 
GROUND  LEVEL  SOURCE  FIELD  DATA  COMPARED  WITH  MODEL  PREDICTIONS 
CORRECTED  FOR  SAMPLE  TIME 


VARIABLE 

FIELD  MEASUREMENTS 
Ts  = 10  to  20  min 

MODEL  PREDICTION* 

Ts  = 3 min 

Ts  = 15  min 

C\J 

o 

II 

CL 

p = 0.3 

Y0T 

Intermittency 

0.64 

1.00 

0.92 

0.84 

1 PT 

PI  ume 
Intensity 

0.95 

0.38 

0.47 

0.52 

V 

Total 

Intensity 

1.40 

0.38 

0.58 

0.71 

* At  average  A = 0.063 


36 


7?  T?  2^ct 
ca  ~ CT  — 


(9) 


7 

for  T > 5<7-r,  where  cT  and  ^7  T are  the  variance  and  scale  for  a sampling 
a cT  I cl 

time  Ts  equal  to  the  averaging  time.  The  time  averaged  mean  cq  is  by 


definition  the  mean  c-j.  for  sample  time  T$ 


T . Dividing  (9)  by  the 

3 

square  of  these  means  yields  the  time  averaged  intensity  i'a  = cfl  /c^ 


. 2 


2 2^cT 


(10) 


Because  increasing  sampling  time  causes  the  seal  e <c7cq-j-  to  increase  due  to 
plume  meandering,  it  is  more  convenient  to  express  (10)  in  terms  of  the 
time  scale  co  of  the  three  minute  plume  using  (7) 


rioV 


L 0 J 


(11) 


2 2 

for  T >5  i <^7  / 1 and  the  ratio  of  intensities  computed  using  (6)  with 

3 0 1 CO  0 

Ts  = T^.  Typically  Acq  - 200  m at  z = 10  m in  the  atmospheric  boundary 

layer  (see  Appendix  C)  so  for  U = 3 m/s  at  10  m,  = Aco/U  is  about 

1 minute.  Usually  we  have  i T/i  < 3,  so  (11)  is  valid  for  T > 45  min. 

o i o a 

In  most  practical  situations  averaging  times  of  an  hour  or  more  are  used, 
and  (11)  is  a good  approximation. 


Influence  of  Averaging  Time  on  Intermi ttency 


To  produce  a zero  concentration  over  an  averaging  interval  the 
instantaneous  concentration  must  be  zero  throughout  the  enti re  averaging 
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period  T . The  chance  of  observing  zero  concentration  will  decrease  with 

increasing  averaging  time,  and  the  intermi ttency  factor  y will  approach 

a 

1.0  for  long  averaging  times.  This  increase  in  y arises  directly  from 

a 

its  definition 


y 


a 


1 + i 
1 + i 


2 

£_ 

2 

a 


(12) 


and  the  decreasing  value  of  i in  (11).  For  total  fluctuation  intensity  of 

a 

iQ  > 3,  (6)  predicts  that 


.2 

V 


C T 


sref 


for  large  i 


(13) 


Also,  for  ip>  1.0  the  conditional  plume  fluctuation  intensity  i^  is  independent 

of  sampling  time  (see  Table  6).  So,  i should  also  be  independent  of  averag- 

pa 

ing  time,  because  T and  T have  complementary  spectral  attenuation  functions 
as 

(see  Pa s qu ill  (1974)  p.  16).  Using  (13)  in  (11)  and  the  result  in  (12) 

2 


1 + i 


^a  = 


1 + 


2 yV, 


_pa_ 


co 


sref 


T ->1-2P 
sref 


(14) 


Comparison  With  Field  Measurements  of  Averaging  Time 

The  full  scale  data  of  Barry  (1977)  in  Figure  15  allows  us  to  estimate 

i and  the  exponent  p in  (14).  Barry  found  that  the  conditional  plume 

intensity  was  approximately  constant  at  i - 1 .00  for  averaging  times  of 

pa 

1 hour  to  100  hour.  This  confirms  the  prediction  of  our  sampling  time 
model  that  i^  remains  constant  with  T$  for  large  fluctuation  intensity  iQ 
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on  the  centerline  of  the  three  minute  plume.  Over  the  range  from 
T = 1 hour  to  100  hours  the  data  in  Figure  15  follows  the  empirical 

d 

relation 


^a  = 


2.0 


1 + 


13.3 
T 0.4 


(15) 


where  is  expressed  in  hours.  Using  the  measured  value  of  ipa  - 1.0  in 

(14)  and  estimating  - 0-02  hr  and  Tsref  - 0.05  hr  we  find  that 

iQ  - 7.4  and  p = 0.3  by  direct  comparison  of  (14)  and  (15).  This  is 
significantly  larger  than  the  value  p = 0.20  recommended  by  Angle  (1978),  and  cl( 
to  the  inertial  subrange  value  of  p = 0.33  suggested  by  Pasquill  (1974) 
p.  65.  It  is  interesting  to  note  that  Barry's  simple  power  law  fit  to  the 
intermittency  in  Figure  15  predicts  that  the  value  of  unity  will  occur  for 
an  averaging  time  of  11  days.  Barry  suggests  that  this  asymptotic  limit 
should  occur  within  about  4 days,  which  is  the  average  time  for  the  wind 
vector  to  swing  through  360°.  From  our  point  of  view,  what  is  most 
surprising  is  that  the  simple  power  law  model  developed  for  the  effective 
sampling  time  appears  to  be  valid  for  sampling  and  averaging  times  of  at 
least  100  hours.  This  conclusion  is  further  supported  by  the  data  of 
McGuire  and  Noll  (1971)  who  found  that  a power  law  with  an  exponent 
of  p = 0.2  applied  for  averaging  times  from  a few  minutes  to  one  year. 
However,  their  data  was  for  an  urban  receptor  surrounded  by  multiple  point 
sources,  and  it  does  not  seem  likely  that  the  same  exponent  should  apply  over 
this  wide  range  of  Ts  for  a single  source.  While  such  long  term  extrapolation 
must  be  subject  to  local  topography  and  climate,  it  appears  to  be  reasonable  to 
use  the  power  law  predictions  of  (4)  through  (11)  for  sampling  and  averaging 
times  up  to  4 days. 


0.1 55  T 


Field  measurements  of  effect  of  data  averaging  time  on  intermi ttency  for  a monitor 
x = 1500  m from  a stack  hs  = 50  m at  Chalk  River  Ontario.  Data  from  Barry  (1977). 


- 40 


Dosage  Fluctuations 

The  prime  purpose  of  using  an  averaging  time  T is  to  express  the 

a 

concentration  as  an  integrated  dose  experienced  by  receptor  over  a specified 
exposure  time.  The  implication  of  using  dosage  as  a measure  of  damage  is 
that  the  receptor  has  a finite  response  time  that  integrates  the  concentra- 
tion history.  The  major  difference  between  a dosage  concept  and  a receptor 
time  constant  is  that  dosage  provides  a memory  of  al 1 previous  events,  while 
the  receptor  time  constant  has  a receding  memory.  If  for  example  the 
receptor  is  exposed  to  a fluctuating  concentration  for  several  hours, 
after  which  the  concentration  remains  at  zero,  a response  function  based 
on  receptor  time  constant  x^  will  indicate  that  the  effective  concentration 
is  zero  after  about  three  time  constants.  A response  factor  based  on  dosage 
will  indicate  that  the  dose  received  remains  constant  after  the  concentration 
has  gone  to  zero. 

It  is  probably  more  accurate  to  express  the  long  term  response  of  a 
receptor  by  using  an  integrated  dose  based  on  a concentration  fluctuation 
history  that  has  been  attenuated  by  a short  term  time  constant  response. 

The  integrated  dose  would  then  represent  the  irreversibl e effects  that  had 
occurred.  A further  refinement  would  be  to  use  a nonlinear  function  of 
concentration  to  represent  the  increasing  sensitivity  and  damage  caused  by 
high  concentration.  However,  given  our  current  state  of  ignorance  on  the 
correlation  between  a fluctuating  concentration  field  and  the  long  term 
effects  on  a receptor,  this  functional  form  is  difficult  to  specify. 

The  smoothing  of  concentration  fluctuations  caused  by  the  long  exposure 

time  over  which  the  dose  is  integrated  is  exactly  the  same  as  the  effect  of 

~T2 

averaging  time.  This  relationship  between  fluctuating  dose  D and  the 
concentration  variance  has  been  derived  in  Appendix  E.  The  result  (E22) 
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is  identical  to  (11)  for  an  averaging  time  T equal  to  the  exposure  time 

a 

T . As  the  exposure  time  increases,  the  fraction  of  the  ensemble  that  has 

zero  dosage  will  decrease  as  it  becomes  more  unlikely  that  concentration 

will  remain  at  zero  value  over  the  entire  exposure  time.  This  dosage 

intermittency  is  identical  to  the  time  averaged  intermi ttency  y given  in 

a 

(14)  when  the  averaging  time  is  equal  to  the  exposure  time  Tg. 

The  conditional  dosage  fluctuation  intensity  is  assumed  to  be  independent 
of  exposure  time  and  equal  to  the  plume  intensity  i ^ for  T^  = Tg.  The 
conditional  mean  dosage  is 


where  the  mean  dosage  is  defined  by  (E8)  for  an  exposure  time  Tg 


T 

e 


with  c-j.  the  mean  concentration  for  a sample  time  T$  = Tg,  from  (4).  The 

average  intermittency  y may  be  computed  from  (14)  with  T = T . Using 

a a e 

the  value  of  i j9  and  y^  the  probability  of  exceeding  some  dosage  peak 
★ 

D may  be  estimated  from  the  probability  distribution. 


Fluctuation  Damping  by  Receptor  Response  Time 

The  response  of  physical  systems  to  fluctuating  inputs  can  often  be 
simulated  by  assuming  that  the  system  has  a time  constant  which  causes 

it  to  lag  behind  sudden  changes  in  input  conditions.  The  effect  of  this 
time  constant  is  to  cause  the  system  to  see  an  effective  concentration 
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fluctuation  variance  c ^ that  is  less  than  the  variance  of  the 

eff 

instantaneous  fluctuations.  For  biological  systems  such  as  trees  and 
other  slow  growing  vegetation,  the  time  constant  may  be  long  enough  to 
effectively  damp  most  concentration  fluctuations  and  cause  the  system  to 
respond  to  the  average  concentration.  The  reduction  in  total  variance 
for  a Markov  spectrum  of  concentration  fluctuations  is  discussed  in 
Appendix  A.  For  this  spectrum  the  effect  of  time  constant  on  variance 
produces  the  remarkably  simple  result 


eff 


(16) 


1 + 


3c  T 


1 2 

where  c-j.  and  are  variance  and  scale  for  sampling  time  T . The 
concentration  c^^  is  unaffected  by  the  time  constant  so  that 


mean 


c 


eff 


(17) 


Dividing  (16)  by  the  square  of  (17)  and  using  (7)  to  allow  the  integral 
seal e^/co  for  the  3 minute  plume  to  be  used.  The  centerline  value  is 
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where  i and  ^7  are  intensity  and  scale  in  the  3 minute  plume.  For 
o ^co 

slowly  responding  systems,  where  »<7co  the  unity  term  in  the  denominator 
may  be  neglected  and  (18)  becomes 
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As  an  illustration  consider  some  vegetation,  which  might  have  a time 

constant  of  - 3 hour.  Typical ly  ^/cq  - 1 minute  near  the  ground,  and 

for  sampling  times  of  several  days  iQ^  - 3 iQ.  Then  from  (19) 

i - 0.45  i which  is  only  about  a factor  of  two  decrease, 
eft  o 


Correcting  for  Instrument  Time  Response 

Another  use  for  the  time  constant  ts  is  to  correct  measured  data 
from  slowly  responding  instruments  to  determine  the  true  concentration 
variance.  Figure  16  shows  the  response  of  an  S0£  monitor  to  concentra- 
tion fluctuations  in  a plume.  The  time  constant  of  the  instrument  is  about 
ts  =1.0  min  so  the  slow  rise  and  fall  of  concentration  at  the  beginning 
and  end  of  episodes  must  be  due  to  plume  meander  and  not  instrument 

response.  Equation  (16)  may  be  applied  to  the  measured  data  to  determine 
~rZ 

the  true  variance  c over  sampling  time  Because  the  instrument  time 

constant  will  also  increase  the  apparent  scale  of  turbulence,  (16)  may  be 
put  in  a more  useful  form  by  using  (A42)  to  determine  the  true  integral 
scale  from  the  effective  found  from  the  measured  data.  Combining 

these  equations  gives  the  simple  correction  factor 


ci 


'2 


ci 


^ci 


(20) 


1 2 1 2 
where  c^  and  xj  are  the  measured  variance  and  integral  scale  and  c 

is  the  true  variance.  The  true  scale  is  related  to  the  measured  scale  by 


(21) 
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Figure  16  Response  of  an  S02  field  monitor  to  concentration  fluctuations. 
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Time  Constants  for  Indoor  Air  Pollution 

A third  useful  application  of  the  system  time  constant  response  is  in 
estimating  the  effect  of  fluctuating  outdoor  concentration  on  indoor 
pollution  levels.  Because  buildings  are  generally  tightly  sealed  and  have 
forced  air  circulation  within  them  they  usually  are  approximated  by  a 
well  mixed  chamber.  If  we  express  the  air  infiltration  into  a well  mixed 
volume  by  the  number  of  air  changes  per  hour  (ACH),  it  is  easy  to  show 
that  the  system  response  time  constant  is 

Ts  = w hours  (22) 

As  an  illustrative  example,  consider  a house  that  is  exposed  to  fluctuating 
concentrations  with  3 minute  plume  intensity  iQ  = i =1.0  so  that  intermi ttency 
factor  yQ  = 1.0,  and  fluctuations  with  an  integral  scale  = 1 minute. 

For  a 12  hour  exposure,  using  an  exponent  p = 0.2  we  find  from  (4),  (5),  (7), 
(D19)  and  (D44)  that  yoj  = 0.40,  iQj=  2.04,  i^  = 1.03.  Then,  fora  typical 
house  with  0.5  ACH  (air  changes  per  hour),  =120  min.  Using  (18),  we  find 
that  the  indoor  fluctuation  intensity  i^^  is  only  18%  of  the  outdoor 
over  a typical  (ensemble  averaged)  12  hour  exposure.  It  is  clearly  better 
to  stay  indoors,  avoiding  the  high  peak  concentrations  that  will  be  felt 


outsi de. 
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CHAPTER  4 


CONDITIONAL  MEAN  CONCENTRATION 


The  major  problem  in  monitoring  atmospheric  pollution  levels  is 

the  intepretation  of  the  highly  intermittent  concentration  fluctuations 

when  the  mean  wind  direction  fails  to  carry  the  plume  directly  over  the 

monitoring  station.  One  possibility  is  to  set  air  quality  standards  on 

the  basis  of  conditional  "plume"  mean  concentration  cp  which  is  determined  by 

discarding  the  time  periods  with  zero  concentration.  Because  the  probability 

★ 

distributions  used  to  estimate  the  fraction  of  time  over  some  threshold  c 
use  conditional  statistics  c^  and  y it  is  sensible  to  express  air  quality 
standards  in  terms  of  this  conditional  value  and  the  intermittency.  From 
the  practical  standpoint  of  monitoring  concentration  readings,  the  condi- 
tional mean  concentration  is  a logical  choice  because  it  remains  large 
even  in  the  intermittent  fringes  off  the  plume  centreline.  For  a single 
variance  source  (<j>  = 0)  the  model  for  short  three  minute  sampling  times 
(D8)  and  (D14)  to  (D17)  predicts  a crosswind  variation  of 
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where  y and  c are  the  plume  centreline  values  for  short  sample  times  of 
oo  r 

Tsref  - 3 minutes.  By  definition,  the  conditional  mean  is 
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so  from  (23)  and  (24) 

c 
c 

pc 

At  large  distances  off  the  plume  centreline  both  the  mean  c and  inter- 
mittency  rapidly  approach  zero.  However,  for  y/a  > 3.0  the  conditional 
mean  concentration  has  an  asymptotic  value  of 


exp 
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Figure  17  shows  these  profiles,  which  have  a constant  c^  in  the  outer 
fringes  of  the  plume.  This  crosswind  variation  is  physically  realistic 
because  we  can  visualize  two  distinct  types  of  dilution  by  entrainment: 


• In  the  outer  fringes  of  the  plume,  dilution  consists 
of  a series  of  random  steps  with  a contaminated  eddy 
entraining  uncontaminated  air.  The  dilution  of  these 
eddies  should  be  a function  of  travel  time,  and  not  how 
far  the  puffs  are  off  the  centerline. 

• Contaminated  eddies  in  the  central  core  of  the  plume  near 
the  centerline  will  experience  a similar  dilution  process, 
but  will  mix  with  other  contaminated  eddies,  and  so  main- 
tain a higher  level  of  concentration  for  the  same  travel 
time. 


The  essential  difference  of  puffs  of  plume  material  entraining  contaminated 
air  on  the  plume  centerline  and  pure  air  in  the  fringes  is  a plausible 
explanation  for  the  shape  of  a crosswind  profile  shown  in  Figure  17.  The 
plume  fluctuation  intensity  i , on  which  the  asymptotic  value  of  the 
conditional  mean  concentration  depends,  depends  on  the  plume  centerline 
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Figure  17  Gaussian  model  prediction  of  conditional  (plume)  means  CD0  across  a plume,  and  for  sampling 
time  meander  averaging  on  centerline  Cp0y,  with  single  variance  source  <f>  = 0 . 
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intermittency  yQ.  Because  plume  fluctuation  intensity  i remains  constant 
across  the  plume  at  a fixed  downwind  position,  (D16)  shows  that  the  inter- 
mittency factor,  yQ  decreases  for  increasing  fluctuation  intensity  iQ.  When 
the  centerline  fluctuations  become  more  intermittent  (yQ  becomes  small), 
the  core  of  the  plume  has  more  eddies  that  are  entraining  pure  air  rather 
than  mixing  with  contaminated  material.  This  causes  the  diluted  eddies  in 
the  plume  core  to  have  concentrations  much  the  same  as  those  in  the  outer 
fringes,  which  are  also  mixing  with  uncontaminated  air.  Thus,  for  the 
high  intermittency  (small  y ) consistent  with  large  centerline  intensity  i 
we  expect  less  variation  in  the  conditional  mean  concentration  across  the 
plume.  This  is  exactly  what  we  see  in  (27),  plotted  in  Figure  17. 

Conditional  Mean  for  Long  Sample  Times 

The  plume  centerline  concentration  for  long  sample  times  will  be 
reduced  by  crosswind  meandering  as  the  sample  time  increases  from  our 
reference  value  of  three  minutes  to  several  hours.  The  meandering  plume 
model  in  Appendix  D expresses  the  meander  width  in  terms  of  the  standard 
deviation  a m of  the  meander  position  probability  distribution.  The 
conditional  mean  concentration  c^  on  the  plume  centerline  for  long  sample 
times  may  be  calculated  from  the  intermittency  y obtained  from  (D19) 
by  numerical  integration,  and  the  nonstationary  mean  c^  from  (D13).  These 
values  are  shown  in  Figure  17,  where  we  see  that  the  centerline  conditional 
mean  c^^  takes  on  an  asymptotic  value  for  long  sampling  times.  For  a 
sampling  time  exponent  of  p = 0.2,  Figure  14  shows  the  rapid  approach  of 
the  plume  centerline  conditional  mean  to  this  constant.  Table  4 gives  the 
asymptotic  values  for  two  different  cases:  the  far  crosswind  fringes  of 
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TABLE  4 

ASYMPTOTIC  VALUES  OF  CONDITIONAL  MEAN  CONCENTRATIONS 
NORMALIZED  BY  THREE  MINUTE  CENTERLINE  VALUE  Cpo 


Centerline  Intensity 

For  Reference 
Sample  Time 

Tsref  = 3 min 
i 

0 

Conditional  ' 

"Plume"  Mean 

cp/cpo 
Sample  Time 
Ts  = 3 min 

y/ay  + 00 

cPoT^cpo 
Sample  Time 

T$  -*■  00 

y/Oy  = 0 

0.01 

0.0001 

0.32 

0.25 

0.059 

0.55 

0.50 

0.20 

0.68 

1.00 

0.50 

0.83  ‘ 

1.50 

0.69 

0.90 

2.00 

0.80 

0.94 
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the  three  minute  sampling  time;  and  the  plume  centerline  for  large 
sampling  times.  In  most  practical  situations  the  centerline  intensity 
i is  greater  than  0.5  so  the  plume  centerline  conditional  mean  does 
not  vary  significantly  with  sample  time.  The  effect  of  increased  sampling 
time  is  mainly  to  cause  the  intermi ttency  factor  to  decrease,  without 
changing  any  of  the  other  conditional  plume  statistics. 


Crosswind  profiles  for  long  sampling  times  can  be  computed  from 
(D45)  to  (D49 ) . These  use  the  same  crosswind  distribution  functions  as 
the  three  minute  plume,  but  with  the  time  averaged  centerline  statistics 
as  normalizing  factors,  and  an  increased  crosswind  spread  a ^ from  the 
power  law  in  (D2) 


O t 

J£l  = 


sref 


(28) 


The  exponential  form  of  the  crosswind  distribution  in  (23)  through  (26) 
makes  them  very  sensitive  to  changes  in  time  averaged  crosswind  spread. 
For  example,  if  the  mean  wind  direction  places  a monitoring  station  at 
y = 3.0  Oy  off  the  axis  of  the  3 minute  plume,  for  p = 0.2  it  will  only 
be  a distance  of  y = 1.09  off  axis  after  an  8 hour  sample  time.  We 
see  that  for  long  sample  times  crosswind  meandering  will  make  off-axis 
monitors  give  a good  estimate  for  the  conditional  mean  value.  Taking  a 
typical  intensity  of  iQ * 1.0  (26)  gives  our  8 hour  off  axis  reading  as 
0-78  c^y.  Thus,  standards  written  for  8 hours  and  based  on  conditional 
values  are  easier  to  enforce,  because  compliance  is  easier  to  monitor. 
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Required  Sampling  Times  for  Conditional  Means 

For  long  sampling  times  the  conditional  mean  remains  relatively 
constant,  which  suggests  that  it  is  a better  measure  of  what  is  happening 
on  the  plume  centerline  than  is  the  mean  concentration  c which  decreases 
rapidly  with  increasing  sampling  time  as  the  plume  becomes  more  intermittent. 
Unfortunately,  we  cannot  escape  the  effects  of  i ntermi ttency  simply  by 
using  conditional  averages.  If  for  example,  we  require  one  hour  of  non- 
zero data  to  obtain  a good  estimate  for  the  conditional  mean  c^,  we  need 
only  wait  an  hour  when  the  intermi ttency  y^  = 1.0.  However,  in  a highly 
intermittent  plume  where  y^  = 0..2,  we  must  wait  five  hours  to  obtain  the 
same  data  record  of  nonzero  concentrations.  The  mean  wind  direction,  speed 
and  stability  often  vary  considerably  over  periods  of  hours,  requiring  several 
measurements  to  form  an  ensemble  average.  For  this  reason,  it  may  not  be 
any  easier  to  obtain  reliable  data  sets  for  conditional  means  than  for 
conventional  mean  values.  What  does  make  the  conditional  mean  value 
useful  is  that  it  is  a better  representation  of  plume  centerline  conditions 
even  when  the  monitoring  station  is  at  an  off-axis  location. 

Variation  of  Conditional  Mean  with  Distance  from  a Stack 

The  model  for  alongwind  variation  of  plume  statistics  from  Appendix  C 

was  combined  with  the  plume  spread  parameters  recommended  by  Angle  (1978)  to 

predict  total  and  conditional  statistics  for  two  elevated  plumes.  The 

results  are  shown  in  Figure  18  for  neutral  (Class  D)  atmospheric  stability 

and  a 3 minute  sample  time.  When  profiles  were  normalized  with  the  ground 

level  maximum  mean  concentration  c and  the  distance  x^  to  this  maximum, 

o m 

profiles  for  stability  Classes  C,  D and  E were  virtually  identical.  The 
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Figure  18  Normalized  plume  (conditional)  and  total  mean  concentrations 
for  3 min.  samples  for  neutral  stability  Class  D dispersion. 
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shape  of  the  curves  will  depend  on  sampling  time,  because  for  large  T$ 

the  conditional  mean  c T becomes  constant  while  the  overall  mean  cT 

pT  T 

continues  to  decrease  as  the  plume  becomes  more  intermittent. 

It  is  clear  from  Figure  18  that  the  conditional  concentration  is 

very  sensitive  to  the  ratio  of  source  diameter  to  source  height.  When  the 

source  is  small  compared  to  its  effective  plume  height,  ground  level 

concentrations  are  highly  intermittent  and  conditional  mean  values  are 

large.  However,  becuase  ground  level  receptors  lie  many  standard 

deviations  below  the  plume  axis  at  positions  close  to  the  stack,  the  model 

is  probably  not  reliable  there.  In  any  case.  Figure  18  shows  that  some  care 

must  be  taken  in  estimating  the  effect  of  the  source  size  to  have  any 

hope  of  accurately  predicting  conditional  statistics.  The  asymptotic 

conditional  plume  concentration  at  locations  far  off  the  axis  is  also 

shown  in  Figure  18.  In  the  highly  intermittent  plume  from  x = 0 to 

x - 0.5  x the  conditional  concentration  c is  relatively  constant  across 

m p 

the  plume,  making  it  easier  to  estimate  these  conditional  concentrations 
when  the  mean  wind  direction  is  not  directly  from  the  source  to  the 
monitoring  station. 

Barry  (1977)  p.366  presents  field  measurements  of  conditional  plume 
statistics  downwind  of  the  Sudbury,  Ontario  smelter  stack.  The  ground  level 
conditional  mean  concentrations  remained  almost  constant  over  the  range  of 
downwind  distance  of  x - 5 to  40  km,  with  maximum  mean  ground  level 
concentration  at  xm  = 1 5 km.  This  is  close  to  the  shape  predicted  in 
Figure  18  for  five  to  one  ratio  of  stack  height  to  source  diameter,  with 
short  sampling  time. 
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CHAPTER  5 


PROBABILITY  DISTRIBUTIONS  AND  THRESHOLD  CROSSING  PROBABILITIES 


We  now  have  a plume  dispersion  model,  verified  using  field  data, 
that  allows  us  to  predict  the  spatial  variation  of  plume  mean  concentration, 
its  variance  and  intermittency.  These  statistics  can  be  adjusted  for  any 
sampling  time  period  from  a few  minutes  to  several  days,  and  corrected  for 
the  smoothing  effects  of  system  time  response  and  averaging  time  for 
integrated  dosage.  The  final  element  required  for  our  fluctuation  model 
is  a probability  distribution  function  which  uses  these  plume  statistics 
to  predict  the  fraction  of  time  that  instantaneous  concentration  or 
integrated  dosage  will  exceed  a specified  limit.  For  exposure  to  highly 
toxic  gases  it  is  also  essential  to  know  the  probability  of  being  exposed 
for  a specified  time  without  even  once  exceeding  a critical 
concentration  limit.  We  will  show  in  this  chapter  that  these  two 
probabilities,  one  of  which  specifies  the  fraction  of  concentrations  that 
will  exceed  some  threshold,  and  the  other  the  probability  of  not 
crossing  the  threshold  are  related  in  some  situations. 

Log-Normal  and  Exponential  Probability  Distributions 

In  the  preceding  report  Wilson  (1982)  found  that  measured  probability 
distributions  of  concentration  fluctuations  of  turbulent  mixing  suggest 
that  there  are  two  physically  different  types  of  dilution  process: 

• In  the  core  of  the  plume  eddies  are  diluted  by  mixing 
with  other  contaminated  eddies.  Here  the  plume  has 
few  intermittent  periods,  and  the  intermittency  factor 
y + 1.0. 
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• In  the  intermittent  fringes  of  the  plume  where  the 
intermittency  factor  y -*  0,  contaminated  eddies  are 
isolated  from  each  other  and  are  diluted  mostly  by 
entrainment  of  uncontaminated  air. 

These  two  processes  should  give  rise  to  different  probability  distributions, 
with  the  core  of  the  flow  following  a log-normal  pdf  and  the  highly  inter- 
mittent fringes  an  exponential  pdf.  Barry  (1977)  supported  this  approach 
by  noting  that 

"We  cannot  expect  one  single  distribution  to  apply  to 
all  cases;  a multiplicity  of  distributions  may  be 
pertinent  depending  on  the  averaging  time  or  averaging 
space." 

In  addition  to  the  different  types  of  dilution  process,  Barry  found  that  the 

averaging  time  used  to  define  an  individual  sample  point  will  also  affect 

the  shape  of  the  probability  distribution.  For  a dilution  process  which 

produces  a random  series  of  uncorrelated  concentration  points  he  showed  that 

an  exponential  pdf  should  result.  However,  when  n sequential  points  were 

averaged  to  simulate  the  effect  of  grouping  data  in  averaging  times,  a Gamma 

distribution  resulted.  While  this  indicates  that  some  care  must  be  taken 

not  to  extrapolate  observed  distributions  far  beyond  the  range  of  averaging 

times  for  which  they  were  observed,  Barry's  field  measurements  downwind 

of  a stack  show  that  for  averaging  times  from  T = 1 to  1 00  hours  there 

a 

is  virtually  no  effect  of  averaging  time  on  the  shape  of  the  distribution 
function.  Very  long  averaging  times  of  several  hours  cause  the  short  term 
pdf  to  overestimate  the  probability  of  exceeding  extreme  values. 

Wilson  (1982)  suggested  a linear  combination  of  exponential  and  log- 
normal distribution  functions  be  used  to  describe  an  intermittent  plume. 
While  this  approach  was  able  to  simulate  the  bimodal  pdf's  found  in  a highly 
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intermittent  jet,  it  is  computationally  cumbersome  and  requires  the 

specification  of  yet  another  adjustable  constant  to  determine  the  relative 

contribution  of  the  two  functions.  To  predict  extreme  values  larger  than 

about  twice  the  conditional  mean  concentration  c it  was  found  that  a 

P 

simple  log-normal  distribution  gave  sufficient  accuracy.  For  concentrations 
less  than  the  conditional  mean  c^  the  log-normal  consistently  under- 
predicted the  probability  of  observing  low  concentrations. 

The  Gamma  Probability  Distribution 

The  Gamma  distribution  is  an  attractive  compromise  between  the  widely 
used  log-normal  pdf  for  plumes  with  y -*  1 .0,  and  the  exponential  pdf 
observed  during  highly  intermittent  conditions  when  y -*  0.  Appendix  G 
presents  a detailed  discussion  of  the  Gamma  distribution  and  shows  it  is 
consistent  with  a physically  realistic  model  for  a dilution  process. 

The  probability  density  has  the  functional  form 


and  r (k)  is  the  gamma  function  that  gives  the  distribution  its  name 


(31) 


where  the  variance  appears  through  k 


(32) 
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The  total  pdf  including  the  spike  at  c = 0 due  to  intermittency  is 

p(c)  = (1  - y)  6 (c)  + yPp(c ) (34) 

where  6(c)  is  the  delta  function.  The  cumulative  pdf  Q,  (c)  is  the  fraction 
of  concentrations  less  than  c and  is  the  integral  of  (34) 


ft(c)  = (1  - y ) + Y^p(c) 


where 


Vc> = 


Pp(c)  dc 


(35) 


(36) 


and  the  complementary-cumulative  fraction  exceeding  c is 


V(c ) =1  - fi(c) 


= Y 0 - ^p(c)) 


(37) 


For  the  Gamma  distribution  (31),  the  integral  in  (36)  must  be  evaluated 
numerically  using  a series  approximation  to  the  incomplete  Gamma  integral, 
see  Appendix  G. 

When  the  plume  intensity  ip  = 1.0,  k = 1.0  the  Gamma  distribution 
reduces  to  the  exponential  form  with  r (1 ) =1.0 


It  is  this  feature  that  allows  the  Gamma  distribution  to  represent  both 
highly  intermittent  and  continuous  concentration  fluctuations.  Bencala 
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and  Seinfeld  (1976)  in  an  analysis  of  concentration  distributions  in 
cities,  noted  that  the  Gamma  and  log-normal  distributions  predicted 
similar  probabilities  for  extreme  values.  This  is  evident  in  Table  5 
which  compares  the  normal  , log-normal  and  Gamma  complementary-cumulative 
distributions  for  exceedance  probabilities. 

Measured  Probability  Distributions 

Wilson  (1982)  reviewed  the  observations  of  several  investigators 
who, with  the  exception  of  Barry  (1977),  fitted  their  data  to  log-normal 
distributions.  However,  most  of  this  data  was  for  multiple  sources  in 
urban  areas.  For  these  area  sources  much  of  the  log-normal  behavior  can 
be  explained  by  a log-normal  distribution  of  wind  speed,  which  is  the 
major  factor  in  determining  urban  pollution  levels. 

Recent  measurements  of  Berger,  Mel  ice  and  Demuth  (1  982)  tested  the 
log-normal  and  Gamma  distributions  on  an  extensive  data  set  for  suspended 
particulate  levels  in  Gent,  Belgium.  Typical  data  sets  of  more  than  seven 
hundred  daily  averages  were  used  to  demonstrate  conclusively  that  the  Gamma 
distribution  was  superior  to  the  log-normal.  They  found  that  the  log-normal 
consistently  underestimated  the  fraction  at  low  concentrations  and  over- 
estimated the  fraction  of  extreme  values  at  high  concentration.  For  their 
eleven  different  measuring  stations  the  conditional  fluctuation  intensity 
ip  ranged  from  0.5  to  0.8  with  an  average  of  0.64.  The  multiple  dis- 
tributed fossil  fuel  sources  which  contributed  to  the  suspended  particulate 
level,  and  the  long  averaging  time  of  24  hours,  allows  us  to  make  the 

assumption  that  the  i ntermi ttency  factor  y was  unity. 

a 
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TABLE  5 

PEAK  EXCEEDANCE  PROBABILITIES  FOR  THREE  FREQUENCY  DISTRIBUTIONS 


Fluctuation 

Intensity 

i 

P 

Concentration  Peak  c* 
Above  Mean 

g*  - Cp 
^;2 

Fracti on 

of  Time  Exceeded 
V(C*) 

Normal 

Log-Normal 

Gamma 

0 

0.50 

0.406 

0.433 

1 

0.158 

0.137 

0.151 

0.5 

2 

0.0228 

0.0442 

0.0424 

3 

0.00136 

0.0148 

0.0103 

4 

0.00003 

0.00521 

0.00229 

0 

0.339 

0.368 

1 

0.106 

0.135 

1.0 

2 

ll 

0.0413 

0.0498 

3 

0.0187 

0.0183 

4 

0.00940- 

0.0674 

0 

0.294 

0.308 

1 

II 

0.0827 

0.117 

1.5 

2 

0.0344 

0.0501 

3 

0.0173 

0.0225 

4 

0.00977 

0.0104 
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The  tendency  of  the  log-normal  to  underestimate  the  fraction  at 
low  concentrations  is  apparent  in  the  wind  tunnel  data  in  Figure  19,  for 
values  less  than  half  the  mean  concentration.  The  Gamma  pdf  gives  good 
agreement  for  both  high  and  low  values,  although  the  peak  appears  to  be 
somewhat  too  large.  Comparing  the  areas  under  the  curves  suggests  that 
this  difference  may  be  due  to  improper  scaling  of  the  probability  measure- 
ments, rather  than  any  deficiency  in  either  distribution.  Both  the 
Gamma  and  the  log-normal  distributions  are  good  approximation  for  the 
extremes  larger  than  twice  the  mean  Cp. 

Hinds  (1969)  measurements  of  short  range  diffusion  from  ground  level 
sources  are  shown  in  Figure  20,  where  the  Gamma  distribution  gives  a 
consistently  better  estimate  than  the  log-normal,  particularly  for  low 
concentrations.  For  long  averaging  times  Barry  (1977)  found  that  the 
conditional  intensity  i was  always  close  to  unity,  and  probability 
distributions  were  exponential.  Because  exponential  and  Gamma 
distributions  are  the  sameif  i = 1.0,  Barry's  data  in  Figure  21  shows 
that  the  Gamma  distribution  is  again  superior  to  the  log-normal  in  predict- 
ing fluctuations  for  time  averaged  data.  Barry's  observation  that  plume 
fluctuation  intensity  i - 1.0  for  long  sampling  times  is  properly 
simulated  by  the  plume  meandering  model  in  Appendix  D.  For  long  sample 
times  the  back  and  forth  meandering  causes  the  crosswind  profile  of  mean 
concentration  c to  generate  fluctuations  on  the  plume  centerline.  When  the 
"internal"  variance  in  the  three  minute  plume  is  relatively  small,  the 
variance  iQ  induced  by  meandering  tends  to  significantly  increase  i^Q  for 
long  sample  times,  as  shown  in  Table  6.  We  see  that  if  "internal" 
fluctuations  with  periods  less  than  three  minutes  are  removed  either  by 


C • p(c),  Normalized  Probability  Density 


62 


■=■  Normalized  Instantaneous  Concentration 


Log-normal  and  Gamma  probability  densities  compared  with 
ground  level  measurements  in  a wind  tunnel  boundary  layer 
h/H  = 0.44,  x/H  = 3.2,  X = 0.175. 


Figure  19 


ft  (c)  Cumulative  Probability 


63 


Relative  Concentration 


Figure  20  Field  measurements  of  exceedance  probability  for  diffusion 
from  a ground  level  source  in  flat  terrain  and  in  the  wake 
of  a building  immediately  downwind  of  the  source,  with  sample 
time  6 to  15  min.  and  averaging  time  of  5 sec.  Data  of  Hinds 
(1969)  in  Barry  (1977). 
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Figure  21  Field  measurements  of  averaging  time  effects  on  exceedance 

probability  at  x = 1500  m from  a stack  h = 50  m.  Data  from 
Barry  (1977).  s 
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TABLE  6 

CHANGES  IN  CENTERLINE  CONDITIONAL  PLUME  FLUCTUATION 
INTENSITIES  WITH  SAMPLING  TIME 


"Internal"  Intensity 

Sample  Time 
T - 3 min 

Yo  : i-00 

Long  Sample  Time 
In tens i ty 

Conditional 

V 

Total 

1oT 

1.50 

1.5 

3.8 

1.00 

1.0 

2.8 

0.75 

0.8 

2.3 

0.50 

0.7 

1.7 

0.10 

0.8 

1.6 

0.05 

0.9 

1.6 

0.01 

1.1 

1.6 
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mixing,  receptor  response,  or  averaging  time,  that  plume  fluctuation 
intensity  tends  to  values  near  unity  due  to  the  meandering-induced 
fluctuations  from  the  mean  profile.  Sawford  (1983),  using  Monte-Carlo 
methods  for  particle  pair  statistics  in  a dispersing  cloud  found  that 
the  total  fluctuation  intensity  ranged  from  about  0.5  to  2.0  depending 
on  the  initial  size  of  the  cloud.  This  further  supports  the  assumption 
of  plume  intensities  near  unity. 

The  more  accurate  prediction  given  by  the  Gamma  distribution  is 
particularly  important  when  the  threshold  for  damage  of  a receptor  is  less 
than  the  mean  concentration,  as  might  be  the  case  for  sensitive  vegetation. 
Even  in  in  a continuous  plume  with  no  intermi ttency,  Tabl e 5 shows  that 
the  concentration  exceeds  the  conditional  mean  c ^ only  30%  to  40%  of 
the  time  for  log-normal  and  Gamma  distributions.  This  emphasizes  the 
need  for  accurate  predictions  of  these  low  concentrations,  and  supports 
our  choice  of  the  Gamma  distribution. 

Threshold  Crossing  Probability 

There  are  many  situations  where  knowing  the  fraction  of  time 
concentration  exceeds  a specified  threshold  c*  is  not  sufficient  to 
determine  risk.  What  is  required  in  some  cases  is  a knowledge  of  the 
probability  of  exceeding  a threshold  level  at  least  once  during  an  exposure 
time  T . The  usefulness  of  this  threshold  crossing  probability  is  obvious 
when  we  consider  exposure  to  a toxic  gas  like  hydrogen  sulfide.  In  this 
case,  the  concept  of  not  once  exceeding  a threshold  is  critical. 

In  Appendix  F we  show  that  this  probability  may  be  estimated  from  the 
cumulative  exceedence  distribution  V in  (37)  using  (F24)  for  the  crossing 
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probability  S 


S(T  ) = exp  f-  V T 1 
e'  K L avg  e J 


(39) 


where  the  average  exceedence  probability  is  for  some  threshold  concentra- 
tion c* 


V (c*)  = ~ 
avg'  Te 


Vf' 


PT 


Vl dt 


(40) 


The  time  dependence  of  the  exceedence  probability  V ( t ) arises  when  the 

concentration  statistics  c y and  i vary  with  time  during  a transient 

P P 

release  of  a toxic  gas.  An  alternative  way  of  looking  at  the  probability 
of  threshold  crossing  is  to  view  in  (4)  as  the  probability  "dose" 

seen  during  passage  of  a transient  plume. 


If  the  plumeis  continuous,  we  will  eventually  exceed  the  threshold 
c*  if  we  wait  long  enough.  However,  because  the  plume  becomes  more  inter- 
mittent as  sample  time  T$  increases,  V decreases  in  direct  proportion  in 
(37).  Combining  (37)  and  (40), 


V 


avg 


dt 


P 


dt 


(41) 


Application  Of  Threshold  Crossing  Probability  to  Air  Quality  Monitoring 

The  concept  of  crossing  probability  can  also  be  directly  applied  to 

planning  field  monitoring  programs.  The  concentration  of  interest  in  this 

case  is  the  sequence  of  time  averaged  concentrations  c with  an  averaging  T 

a a 
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equal  to  the  time  interval  over  which  the  air  quality  standards  are  set. 

For  this  sequence  of  time  averaged  values  we  want  to  know  the  length  of 

time  we  must  wait  for  a violation  of  the  standard  c*  to  occur.  Because 

plume  statistics  change  as  sample  time  increases,  we  must  carry  out  an 

integrated  average  of  (41)  to  find  VQVg.  This  value,  determined  by 

numerical  integration  of  (41)  is  used  in  (39)  to  determine  the  probability 

of  completing  the  exposure  time  Te  from  the  beginning  of  monitoring  to  the 

present  moment  without  having  a violation.  It  is  clear  from  both  the 

probability  distribution,  and  from  common  sense,  that  when  ground  level 

concentrations  c are  a factor  of  two  smaller  than  the  limit  c*, 
a 

it  may  be  necessary  to  sample  for  a long  period  before  seeing  the  violation 
on  a fixed  monitor. 

One  troublesome  approximation  in  Appendix  F is  the  assumption  (F6) 
that  the  hazard  rate  is  equal  to  the  fractional  exceedence  probability. 
Because  this  approximation  does  not  explicitly  state  the  time  units,  we 
must  make  a logical  choice  for  the  units  used  for  time  measurements  in  (39). 
One  reasonable  choice  is  to  take  the  averaging  time  over  which  successive 
steps  are  taken  as  one  time  unit,  so  that  exposure  time  for  one  hour  standards 
is  also  expressed  in  hours.  For  example,  if  the  adjusted  statistics  for 
a one  hour  exposure  indicate  that  the  air  quality  standard  will  be  exceeded 
10%  of  the  time,  then  V&  = 0.1,  and  (39)  indicates  that  we  must  wait  a 
total  of  46  hours  to  be  99%  certain  of  seeing  this  violation.  Flowever,  this 
choice  of  time  units  is  rather  arbitrary  and  illustrates  the  difficulty  in 
inferring  the  time  series  behavior  from  the  fractional  exceedence  pdf  without 
making  many  restrictive  and  often  questionable  assumptions. 
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EFFECT  OF  RECEPTOR  TIME  CONSTANT  AND  DATA  AVERAGING  TIME 
INTERVAL  ON  EFFECTIVE  FLUCTUATION  VARIANCE 

All  measuring  systems  have  a limit  to  their  frequency  response,  usually 
expressed  as  the  -3dB  attenuation  point  in  their  amplitude  gain,  G.  For  a 
response  (i.e.  "output")  amplitude  AQut  to  a system  input  amplitude  A-n  the 
"gain"  is 

G = jsui  (Al) 

i n 

Using  the  definition  of  attenuation  decibels 

dB  = 10  log  G (A2) 

and  we  see  that  the  -3dB  point  on  the  frequency  response  curve  is  the  location 

where  G = 0.5,  and  the  amplitude  has  decayed  by  a factor  of  two.  For  many 

electronic  systems  this  definition  is  often  confused  with  the  "power  bandwidth" 

of  the  device.  Here,  the  amplitude  A is  voltage,  and  because  power  output  into 

2 

a constant  resistive  load  is  proportional  to  the  square  of  the  voltage  A , 
the  -3dB  attenuation  on  the  power  occurs  when 


0.5 


or 


A 4. 
out  _ 

" A 

in 


0.707 


The  listed  frequency  response  of  electronic  devices  often  fails  to  distinguish 
between  amplitude  and  power  bandwidths  and  some  care  is  required  in  interpreting 
specified  -3dB  frequency  limits. 
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System  Frequency  Response 

Consider  the  frequency  response  of  a simple  resistance-capacitance 
system  which  is  described  by  a single  time  constant  ts>  with  no 
dissipative  damping.  The  response  of  this  measuring  system  to  a step  change 
in  input  from  zero  to  Vin  is  an  exponential  rise  with  time,t 


Vout  = Vin(1-eXp(Ts)) 


(A3) 


For  a sinusoidal  input  with  amplitude  A.  and  frequency  "f"  cycles  per 
second,  the  gain  G is 


out 

A. 

in 


1/2 


1 + (t$  2rrf ) ‘ 


(A4) 


The  -3dB  attenuation  frequency  occurs  where  G = 0.5  at  a frequency  fg  ^ which 
from  (A4)  is 


f 

0.5  2tt 


(A5) 


The  Markovian  Spectrum 

If  a variable,  c1  , with  no  inertia  is  subjected  to  random  fluctuations, 
the  auto  correlation  of  these  fluctuations  for  varying  time  delay,  t',  will 
be  exponential.  The  essential  characteristic  of  such  inertialess  systems 
is  that  the  time  rate  of  change  of  c'  has  no  autocorrelation.  That  is. 


ac'  (t)  3c'  ( t+t ' ) 

at  * at 


is  zero  except  when  t'  = 0.  In  simple  terms  this  says  that  knowing  c' 
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is  increasing  at  time  t tells  us  nothing  about  whether  c'  will  continue  to 
increase,  stay  constant,  or  decrease  in  the  next  instant.  Hence  the  term 
"inertialess".  Such  first  order  Markov  processes  have  an  exponential  auto  time 
correlation  given  by 


Rc(t')  = exp 


where 


Rc(t) 


c ' ( t)  c'(t+t') 
c‘2 


(A6) 


( A7 ) 


The  constant  is  the  integral  time  scale  of  concentration  defined  by 


Rc(f) 


dt 


( A8) 


Frozen  Turbulence 

If  the  turbulent  eddies  are  not  too  small,  they  lose  their  coherence 
slowly  as  they  are  advected  past  a point.  In  this  case  we  can  assume  that  the 
turbulence  is  "frozen"  into  the  flow  and  is  carried  past  as  point  at  the 
local  mean  velocity,  U.  Because  it  is  the  size  of  eddies  not  their  time  of 
passage  that  remains  constant  when  the  speed  varies,  it  is  more  useful  to 
think  in  terms  of  the  concentration  integral  length  scale  A^.  For  turbulence 
frozen  into  the  flow  which  passes  at  speed  U, 

Ac  = U^7c  (A9) 


In  all  the  analysis  here,  a non-intermi ttent  turbulence  is  assumed.  That  is 
c > 0 and  y = 1.0  at  all  times.  For  intermittent  fluctuations  adjustments 
are  required  to  deal  with  the  periods  of  i ntermi ttency . 
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Freguency  Spectrum  of  Concentration  Fluctuations 

To  see  the  effect  of  frequency  dependent  gain  G on  the  input  concentration 
fluctuations,  we  first  use  a Fourier  transform  to  express  the  exponential 
correlation  as  a frequency  spectrum  Ec(f),  for  frequency,  f 

T r 00 

Ec(f)  = 4C  Rc(t ' ) cos  ( 2irf t ' ) dt1  (A10) 

o 


which  for  the  exponential  correlation  is 


E 

c 


(f) 


o 


exp 


cos  (2nft 1 ) dt ' 


(All) 


Integrating  this  twice  by  parts  gives  the  Markov  spectrum 
4 c'2^c 

E (f ) = (A12) 

C 1 + (^2-nfr 


This  can  be  expressed  in  terms  of  the  integral  length  scale  of  concentration 
fluctuations  Ac  by  using  (A9) 


U Ec(f)  = 


1 + (Ac 


2rf  T 

U 1 


(A13) 


For  turbulence  "frozen"  into  the  mean  flow  the  fluctuation  frequency  f 
increases  in  direct  proportion  to  the  flow  speed  U.  The  wave  number 
defined  as 


ki  = ^ T — H fir  (A14> 

I U wavelength 


remains  constant  with  changes  in  flow  speed 
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Observed  Spectrum  and  Variance  with  Gain  Attenuation 

The  spectrum  is  defined  in  terms  of  the  square  of  the  fluctuation 
amplitude  c'  about  the  mean.  The  spectrum  will  be  attenuated  by  the  square 
of  the  system  response  gain  G so  that  the  observed  "effective"  spectrum  is 

EC;eff(f)  = G2(f)  • Ec(f)  (A15) 


Combining  (A4)  and  A ( 1 2 ) 
Ec,eff  (f) 


1 + (t  2nff  - 


4 C'2^c 


1 + (^7.  2TTf ) 


( AT  6 ) 


This  system  response  attenuation  is  shown  schematically  in  Figure  A1 . 
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To  determine  the  "effective"  variance  c'  ^ observed  by  the  system, 
the  integral  of  the  spectral  density  is  taken 


Ec,eff 

o 


(f)  df 


(A17) 


Using  (A4)  and  (A16),  noting  that  the  integral  time  scale  ^7cis  a function 
of  location  and  not  frequency  f. 


,2 

eff 


= 4^7 


df 

0 (1  + (27TTs)2f2)  (1  + ( 27T  ^7c  ) 2f 2 ) 


To  evaluate  the  integral,  Bara  (1983)  defines 
1 


a = 


2tt  t 


and 


b = 


(A18) 


(A19) 


(A20) 
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Then  factor  these  parameters  out  of  the  integrand  to  obtain 


- 4aV^c 


oo 


0 


df 

(a2  + f2)(b2  + f2) 


(A21) 


Then  separate  the  integrand  into  two  parts  using  partial  fractions 


(a2  + f2)(b2  + f2)  (a2  + f2)  + (b2  + f2) 

? ? 2 ? 

1 = Ab  + Af  + Ba  + Bf 

1 = (Ab2  + Ba2)  + (A  + B)f2 
The  only  solution  of  which  must  be 


(A22) 


A + B = 0 
Ab2  + Ba2  = 1 


which  leads  to 


.2  ,2 

b - a 


B = 


*2  h2 

a - b 


so  that  (A22)  becomes 


(a2  + f2)(b2  + f2)  (b2  - a2)(a2  + f2)  (a2  - b2)(b2  + f2 


(A23) 


The  integral  in  (A21 ) then  becomes 


ceff 


(b2  - a2 


df 


4a2b2<7 


S—  + 


(a2  + f2)  (a2  - b2) 


df 


( A24) 
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The  integrals  can  how  be  evaluated  from  standard  integral  tables  as 


df 


0 (a2  + f2)  a 


- tan 


[ f 1 
a 


with  a similar  equation  for  the  second  integral 
Using  these  in  (A24), 


f=c 


eff  = 4a2b2^7 


,2 


a(b2  - a2) 


tan 


(A 


f=0 


noting  tan  \o)  =0  and  tan  "*(<»)  = u/2 


2tt  a2b2  J 


a(b2  - a^)  b(a2  - b2) 


2tt  a2b2^7 


ab(a  + b) 


From  (A20)  ^<JC  = 1 /b  so  that  the  equation  reduces  to 


,2 


eff  1 


^ 0 + 


Using  the  definitions  in  (A19)  and  (A20)  this  reduces  to  the  simple  result 


(A26) 


The  spectral  attenuation  is  shown  in  Figure  A1 . 
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This  can  be  interpreted  in  two  ways.  When  applied  to  the  response  of  a 

2 

sensing  instrument  it  shows  that  the  actual  variance  c'  is  just  the  measured 

variance  multi  pled  by  (1  + TS/^7C)-  On  the  other  hand,  for  a biological 

receptor  with  a response  time  constant  of  t$ , the  fluctuation  level  seen  by 

2 

the  receptor  is  just  the  actual  variance  c‘  di vided  by  this  same  factor  of 

O + Ts/^7c) 

Effect  of  Data  Averaging  Time  Interval  on  Observed  Variance 

As  we  have  shown  in  the  previous  analysis,  the  response  time  constant, 
of  the  receptor  (a  measuring  instrument,  or  biological  receptor)  causes 
the  observed  concentration  fluctuation  variance  to  be  less  than  the  actual 
variance.  This  time  constant  can  seldom  be  controlled,  and  must  be  considered 
as  part  of  the  system. 

The  variance  is  also  reduced  by  another  quite  different  parameter;  the 

averaging  time  selected  for  sequential  data  samples.  Because  of  the  hour 

to  hour  variability  of  wind  speed,  direction  and  stability,  full  scale  data 

is  collected  and  analyzed  over  short  intervals  of  duration  T . Because  the 

a 

variance  computed  over  an  interval  T will  exclude  low  frequency  fluctuations 

a 

with  frequencies  less  than  1/T  , some  of  the  full  scale  variance  will  be 

a 

missed.  The  problem  we  will  deal  with  here  is  how  to  estimate  the  actual 

atmospheric  fluctuation  variance  c from  the  measured  variance  c , an 

a 

ensemble  average  of  variance  each  of  which  is  measured  over  a time  interval 
T . Here  we  will  assume  that  the  instantaneous  data  used  to  compute  the  sample 

a 

averages  is  taken  with  an  ideal  sensor  that  has  a zero  response  time  constant, 
ts  = 0,  so  that  there  is  no  attenuation  due  to  sensor  response. 
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The  effects  of  sample  averaging  time  T are  similar  to  those  of 

a 

response  time  constant  x . In  both  cases: 

1. )  The  peak  concentration  excursionsare  reduced 

2. )  The  observed  periods  of  zero  concentration 

(i.e.  intermittency)  occur  for  a smaller  fraction 
of  time,  so  observed  y is  larger. 

3. )  The  observed  variance  is  reduced. 

These  effects  are  illustrated  in  Figure  A2.  In  the  following  analysis, 

only  the  influence  of  averaging  time  on  variance  will  be  estimated.  Pasquill 

(1974,  pp.  12-16)  discusses  the  attenuation  of  variance  for  averaging  time 

T by  analyzing  the  fluctuation  spectrum.  The  attenuation  factor  at 
a 

frequency  f for  the  variance  is  given  by  the  square  of  the  gain  G,  where 
sin(Tr  f T ) 

G(f)  = “Tn (A27) 

3 

This  may  be  compared  to  (A4)  which  is  the  gain  for  a receptor  with  a response 
time  constant  of  x . The  attenuated  concentration  fluctuation  spectrum 
E_  a(f)  is  given  by: 

Ec>a(f)  = G2(f)  • Ec(f)  ( A28) 


The  variance  for  samples  of  duration  Ta  is  simply  the  integral  of  the  spectrum 


Ec  Jf)  df 

c ,a 


(A29) 


o 
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This  variance  is  the  smoothed  average  of  a large  ensemble  of  samples  each  of 

duration  T . For  the  Markov  spectrum  in  (A12),  Shoji  and  Tsukatani  (1973) 

2 

have  evaluated  the  integral  of  G • Ec  using  (A27)  to  obtain 


(A30) 


For  very  long  averaging  times  relative  to  the  fluctuation  time  scale,  the 
exponential  approaches  zero  and  the  ratio  of  the  observed  to  actual  variance 
reduces  to 


for  large 


(A31) 


For  short  times  the  exponential  may  be  expanded  as  a Taylor 


T = 0 to  obtain  to  fourth  order 

a 


exp 


LV 


Inserting  this  in  (A30)  yields 


for  smal 


3c 


series  about 


(A32) 


(A33) 


which  is  correct  to  fourth  order  in  . This  series  approximation  in  (A33) 
will  be  within  ]%  of  the  exact  ratio  predicted  by  (A30)  for  ( Ta/<7C )<  1.0. 
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The  ratio  of  observed  to  actual  variance  can  also  be  computed  directly 
from  the  time  correlation  Rc ( t 1 ) defined  in  (A7).  This  approach  gives  a 
better  physical  insight  into  the  effect  of  time  averaging  interval  than  does 
the  spectral  analysis  method.  Fackrell  (1978)  presents  a concise  derivation 
which  we  will  reproduce  here.  The  concentration  fluctuation  c'  about  the 
mean  c when  averaged  over  the  time  interval  T^  is 


o 

The  auto  correlation  of  this  is 


(A34) 


;2=u 


c (t)  dt 


a , ~ 
c (t)  dt 


TaT 


c (t ) c ' (t)  dt  dt 


(A35) 


o o 


For  a stationary  random  variable  we  can  define  t 1 = t - t without  loss  of 
generality,  and  then  take  the  ensemble  average  to  write 


T _ 1 

• V 


Va 


c ' ( t ) c' ( t+t ' ) dt  dt' 


(A36) 


o o 


Using  the  definition  of  the  auto  correlation  R ( t ' ) in  (A7)  this  becomes. 


T T 


R (t‘)  dt  dt' 
c 


(A37) 
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This  can  be  maniputed  using  Kampd  de  Fdriet's  transformation  (see  Hinze 
(1975),  p.52)  to  a single  integral 


( 1 - R_ ( t 1 ) dt ‘ (A38) 

a L 
o 

For  the  exponential  correlation  in  (A6)  that  produces  the  Markov  spectrum, 

(A38)  can  be  integrated  by  parts  to  produce  (A30). 

Effect  of  Response  Time  Constant  and  Averaging  Time  on  Measured  Fluctuation 
Scales 

Both  the  system  response  time  constant  x and  the  averaging  time  T 

s a 

attenuate  the  high  frequency  end  of  the  fluctuation  spectrum,  leaving  the 
low  frequency  fluctuations  unaffected.  This  high  frequency  attenuation 
will  cause  the  measured  concentration  integral  scale  <Jceff  or  <J/ca  to  be 
larger  than  the  true  value^.  These  integral  time  scales  are  related  to 
the  fluctuation  variance  through  their  spectra  by  the  relation  from  Appendix 
C,  equation  (C3) 

Ec(0)  = 4 c '2<7c  ( A39 ) 

Because  the  averaging  time  T and  the  system  response  time  constant  affect 

a 

only  the  high  frequencies,  this  zero  frequency  intercept  Ec(0)  will  be  the 

1 2 1 2 

same  for  spectra  of  c „ and  c with  time  constant  x or  averaging  time  T . 

e t t a s a 

Equating  the  spectral  densities  Ec(0)  yields  the  simple  relationship  from  (A40) 


( A40) 
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( A41 ) 

Using  (A26)  and  (A40)  produces  the  remarkable  result  that  the  effective 
fluctuation  time  scale  ^ceff  is  simply  the  sum  of  the  true  time  scale  and 
the  system  time  constant 


^ceff  ^7c  + Ts 


(A42) 


for  the  Markov  spectrum.  The  effect  of  averaging  time  T may  be  found  by 

a 

combining  (A30)  and  (A41).  However,  in  many  situations  the  averaging  time  is 
much  larger  than  the  fluctuation  scale^,  and  the  approximate  form  in  (A31) 
may  be  used.  This  leads  to  the  simple  result  that  the  measured  integral  scale 
of  time  averaged  data  is 


because  for  long  averaging  times  the  measured  integral  time  scaled/  has 

ca 

nothing  whatever  to  do  with  the  actual  integral  scale  <7c,  because  the 
fluctuations  which  contribute  to  the  scale  ^ have  been  smoothed  out  by 
the  averaging  process. 


Normalized  Spectral  Density 
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Figure  A1  Attenuation  of  the  High  Frequency  End  of  the  Concentration 
Fluctuation  Spectrum  by  Receptor  Time  Response. 
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Figure  A2  Influence  of  Sensor  Response  on  Observed  Interim' ttency  in 
Water  Channel  Simulation  of  the  Plume  from  a Ground  Level 
Source. 


APPENDIX  B 


UNIVERSAL  GROUND  LEVEL  PROFILES  OF  MEAN 
AND  FLUCTUATING  CONCENTRATION 

The  objective  of  this  note  will  be  put  the  along-wind  variations 
of  mean  concentration  and  fluctuation  intensities  in  terms  of  the  single 
dimensionless  plume  spread  parameter  A,  where 

/ a a 

(Bl) 

and  H is  the  boundary  layer  thickness. 

Mean  Concentration 

For  a Gaussian  plume  with  plume  spreads  in  the  form  of  power 

laws. 


a 

c = D x y (B2) 

y y 

a 

a = D x z (B3) 

z z 


Pasquill  (1974,  p.273)  and  Weil  and  Jepson  (1977)  find  the  maximum  ground 
level  concentration  CQm  for  a given  windspreed  U and  constant  source  height 
h to  be 
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om 
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it  exp(|-)  Uhe 


(Dz/e) 

D D 

z y 
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(B4) 


at  the  downwind  location  x_ 


m 


x = 
m 


h 

D~7t 


1/a. 


(B5) 


where 


e = 1 + -^ 
a 

z 


(B6) 


The  alongwind  variation  of  ground  level  mean  concentration  on  the  plume 
axis  is  given  by  the  normalized  equation 


c ea 

om  /_x_\  z 
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exp 


r x "2a 
(— ) z- 1 
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(B7) 


given  by  Weil  and  Jepson  (1977).  These  equations  can  be  transformed  to 
plume  co-ordinates  by  combining  (Bl),  (B2)  and  (B3)  to  write 


0 D D a + a 

X2  = X y Z 


(B8) 


Then,  evaluating  (B8)  at  the  position  of  maximum  mean  concentration. 


to  find  A 
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9 D D a + a 

x 2 = -x_z  . y z 
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o x_ 

,2  m 


( B9) 


Combining  this  with  (B5) 
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from  which 


Xm  = 


(D  D ) 
x y z' 

H 


0.5  r 


L ^ 


Taking  the  ratio  of  (B8)  to  (B9)  yields 

2 


a + a 

y z 


i Am 
^ m 


(BIO) 


(BID 


Using  this  to  transform  (B7)  gives  the  normalized  mean  concentration  as 


(B12) 


By  solving  (B9)  for  and  inserting  the  result  in  (B4),  the  maximum  ground 

level  concentration  may  be  expressed  as  a function  of  A m 
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om 


(1 

TT  exp(f)  UAm2  H2 


(B13) 


The  appearance  of  the  boundary  layer  thickness,  H,  is  a result  of  the  de- 
finition of  A,  which  uses  H as  a normalizing  factor.  Wilson,  Fackrell  and 
Robins  (1982)  suggest  that  in  neutral  stability  H - 650  m.  In  terms  of  the 
neutrally  stable  log-law  surface  roughness  ZQ  they  recommend,  as  a rough 
estimate 

H = 2500  ZQ  (B14) 


While  the  value  of  H will  affect  the  shape  of  the  alongwind  concentration 


-64- 


prof  i 1 e because  it  changes  the  scaling  between  A and  x in  (B9),  it  has  no 

2 2 

effect  on  the  maximum  concentration  c because  A H = a a , which 

om  m ym  zm 

has  no  direct  dependence  on  H. 

Fluctuation  Intensity 

The  square  of  the  fluctuation  intensity  is  given  by 
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(B15) 


At  ground  level,  where  z = 0 the  vertical  and  crosswind  distribution  functions 
reduce  to 


Fm(0)  = exp 


Gm(^-)  = exp 


Fy(0)  = exp 
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Gv  = 7 exp  f-  Ko1,  - 
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(B16) 


(B17) 


(BIS) 


(B19) 


The  effective  height,  hv,  of  the  variance  source  in  (B18)  is 


h = 
v 


f.2  , A 2 2 

h + <PZ  oz 


0.5 


(B20) 
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Using  (B20)  in  (B18),  and  then  taking  the  ratio  of  (B18)  to  the  square 
of  (B16) 
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The  ratio  of  crosswind  distribution  functions  from  (B19)  and  (B17)  is 
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Using  these  in  (B15) 
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The  ratio  h/az  can  be  expressed  in  terms  of  plume  size  X by -writing  (B5)  as 


h = /I  D x 
z m 


(B24) 


Using  a from  (B3)  the  ratio  is 


= /T 


(B25) 
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Then  using  the  relation  between  x and  A in  (Bll) 


/T 


(B27) 


Using  this  (B23)  becomes 
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Then,  along  the  plume  centerline  where  y - 0,i  = i the  axis  ground  level 

intensity,  from  (B28) 
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i02  = (J)  exp 


Am  e <b._  -] 
£ f_E\  _ 

2 A 2 


ex 


p(-VKV)  (B29) 


Conditional  Fluctuation  Intensity  and  Intermi ttency 

Far  above  the  influence  of  ground  surface  effects,  the  conditional 
fluctuation  intensity  ip^  is  constant  across  a plume  cross  section  at 
a fixed  x (or  A)  location.  It  is  given  by 


i 2 = (l)2  exp  (-  4-) 

poo  'Q'oo  r d 


(B30) 


where  the  variance  source  strength  is  given  by 
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0.337  X 


X2  + ( 0 . 1 1 9 ) 2 


exp(£) 


(B31 ) 
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The  effect  of  the  receptor  height  z above  ground  level  on  the  conditionally 
averaged  fluction  intensity  i is  given  by 


(B32) 


At  z = 0 this  reduces  to  the  simple  result  that 


(B33) 


The  plume  intermi ttency  may  be  computed  from  the  general  relationship 


.'hi 

i1 2  + 1 


(B34) 


Variance  Source  Strength 

The  source  strength  is  given  by 


3.  = 
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0.337  X 

2 2 
A + Ag 


exp(^-) 


(B35) 


with  the  virtual  origin  for  fluctuations  at  Xq, 


Ao  = 


0.119 


1 + 0.033 


with  the  reference  height  given  by 


(B36) 
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This  can  be  expressed  in  terms  of  the  plume  size  parameter  A by  using 
(B27) 


exp 


1 (A) 

e 'V 
m 


(B38) 


From  this  we  see  that  the  only  parameter  other  than  A and  A needed  to  charact 

m 

ize  the  normalized  alongwind  mean  and  fluctuating  concentration  on  the  plume 

axis^is  the  source  height  to  diameter  ratio  h/d  in  (B38).  The  value  of  A is 
’ m 

computed  from  source  height  h and  the  plume  spread  constants  using  (BIO). 


APPENDIX  C 


LENGTH  AND  TIME  SCALES  OF  CONCENTRATION  FLUCTUATIONS 


The  time  scale  of  concentration  fluctuations  is  defined  by 


^ = 


R (t ' )dt 
c 


(Cl 


Because  intermiitency  is  like  a square  wave  distribution  superposed  on  plume 
fluctuations,  the  random  frequency  of  the  i ntermi ttency  will  be  superposed  on 
the  plume  fluctuations,  and.  contribute  to  both  spectrum  and  scale.  For  this 

To  TJ 

reason  we  use  c rather  than  the  conditional  average  cp  to  characterize  scale. 
The  spectrum  and  correlation  are  related  by  the  Fourier  integral  transform 


Ec(f)  = 4^ 


Rc  ( t ' ) cos  (2irft 1 ) dt ' 


(C2) 


The  time  scale  may  be  obtained  from  the  zero  frequency  intercept  of  the  spectrum, 
obtained  experimentally  by  extrapolating  to  f = 0,  and  in  (C2) simply  by  setting 
f = 0 to  obtain;  noting  cos(0)  = 1.0. 


Ec(0)  = 4^ 


4c 


R (t ' ) dt 1 
cv 


( C 3 ) 


Assuming  frozen  turbulence,  convected  by  at  the  local  mean  velocity  U,  we  can 
! define  a length  scale  A 


K = ^7cu 


(C4) 
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Unlike  the  Eulerian  time  scaler?  which  becomes  shorter  as  the  convection 

°c 

velocity  U increases,  the  length  scale  Ac  measures  the  physical  size  of 
concentration  eddies,  and  remains  constant  at  a given  location  in  the  plume. 
Combining  (C3)  and  (C4) 


A 


c 


Ec(0).  U 


(C5) 


Variation  of  Length  Scale  Across  a Plume 

Because  only  a few  measurements  of  length  scale  are  available  from  the 
wind  tunnel  data  of  Fackrell  and  Robins  (1982),  a semi -empirical  model  must  be 
developed  to  interpolate  between, and  extrapolate  from  existing  data.  Before 
proposing  a model,  it  is  worthwhile  reviewing  Fackrell  and  Robins  observations 
of  concentration  scale  for  a ground  level  and  elevated  source.  Their  data  has 
been  transformed  from  a normalized  time  scale  <7  Uu/H  to  length  scale  A0/H  by 

^ C n C 

using  (C4)  and  the  mean  velocity  profile  in  their  boundary  layer.  The  data  on 

the  plume  centerline  y = 0 is  presented  in  Figures  Cl  and  C2  for  the  ground 

level  and  elevated  sources.  The  measurements  have  more  scatter  than  the  mean  c 
~r2 

and  variance  c measurements,  emphasizing  the  difficulty  of  measuring 
turbulence  scale.  This  scatter  is  due  to  the  long  sampling  periods  required 
to  define  the  low  frequency  end  of  the  spectrum  needed  for  extrapolation  to 
get  E (0)  in  (C5).  In  spite  of  this,  three  features  are  clear: 

1)  The  small  scale  turbulence  and  mean  wind  shear  near  the 

ground  surface  dissipates  and  smears  out  the  small  scale 

concentration  eddies  in  the  plume.  This  leaves  only 

large  scale  fluctuations  and  causes  the  dramatic  increase 

in  A near  the  surface, 
c 


-C3  - 


2)  Far  above  the  surface  dissipation  effects,  where  either 
z foz  or  h/a2  is  large,  the  length  scale  remains  relatively 
constant  across  a plume  cross  section  at  a fixed  downwind 
distance.  Measurements  at  varying  crosswind  y/o^  locations 
confirm  this. 


3)  The  scale  AQ  away  from  the  influence  of  surface  effects 
increases  gradually  with  downwind  distance  as  the  plume 
spread  becomes  larger. 

The  effect  of  surface  dissipation  on  the  variance  is  accounted  for  by  a variance 

~r2 

sink  term  in  the  Gaussian  model  for  c . Because  these  surface  effects  decrease 
~r2 

c while  increasing  the  scale  Ac»  we  will  assume  in  the  model  that  the  effects 
are  directly  related  through  their  effect  on  the  spectrum. 

~~r2 

As  a physical  model  of  the  relationship  between  Ac  and  c we  assume 
that  the  decreasing  turbulence  scale  and  the  increasing  mean  wind  shear 
dU/dz  near  the  surface  act  only  on  the  smaller  concentration  eddies,  leaving 
the  large  eddies  at  the  low  frequency  end  of  the  spectrum  unaffected.  Because 
these  surface  related  dissipation  effects  act  only  on  the  high  frequency  end 
of  the  spectrum,  the  zero  frequency  intercept  Ec ( 0 ) should  be  the  same  at 

p 

all  heights  above  ground  when  normalized  by  c , the  square  of  the  mean 
concentration.  Then  (C5)  can  be  written  as 


Ec(°)  U 

-2 
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4c  A( 

-2 
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(C6) 


which  is  constant.  Equating  values  far  above  the  ground  where  surface 
dissipation  effects  are  negligible  this  leads  to 


-C4  - 


.2  . 2 
l A.  = i 
£ 00  £00 


(C7 ) 


where  the  subscript  refers  to  the  same  location  but  with  the  ground 
plane  and  surface  dissipation  effects  removed,  and  "i"  is  the  fluctuation 
in ten si ty 

i = ~~  (C8) 


Wilson  (1982  gives  the  ground  surface  dissipation  effect  on  the 
fluctuation  intensity  as 

h.  z , Bo. 5 


1 - a exp 
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1 + exp 
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h z 
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where  hv  is  the  effective  variance  source  height,  which  is  related  to 


the  actual  source  height  by 

[ 


hv  = I h2  + ^\Z 


(CIO) 


The  values  of  <J)  = 0.706  and  a = 0.6  were  derived  by  Wilson,  Robins,  and 
Fackrell  (1982)  from  experiment.  Using  (C9)  in  (C7)  gives  the  scale  ratio 
as 
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1 + exp 

9 h z 1 
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1 - a exp 


- 2 
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It  is  interesting  to  note  that  regardless  of  source  height,  the  ground 
level  concentration  fluctuation  scale  at  z = 0 always  takes  on  the 
constant  values  of 


(C12) 


For  a = 0.6  this  ratio  is  exactly  10.0,  a very  large  increase  indeed. 

Alongwind  Variation  of  Concentration  Scale 

The  variation  of  concentration  scale  with  downwind  distance  is  a 

complex  process  in  which  the  turbulence  scales  of  the  velocity  fluctuations 
~2  T ~2 

u , v and  w interact  with  source  size  effects  and  decay  of  concentration 
fluctuations  to  produce  the  scale  Ac-  Sirivat  and  Warhaft  (1982)  measured 
the  decay  of  concentration  variance  of  tracer  gas  injected  uniformly  from  a 
plane  source  in  grid  turbulence.  They  found  that  the  variance  decay  rate  was 
a function  of  the  ratio  of  concentration  length  scale  to  velocity  length  scale, 
showing  the  interaction.  The  semi -empi ri cal  model  proposed  here  is 


A 

-jp  = B2(X  + X0)  (Cl  3) 

The  effect  of  velocity  fluctuation  scale  is  contained  implicitly  in  A 
! because 


/a  a 


( Cl  4) 
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! and  by  dimensional  reasoning 

; fluctuation  length  scales  A 
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the  plume  spreads  and  oz  contain  the  velocity 

and  A and  turbulence  intensities 
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where  Bq  is  the  Lagrangian-Eulerian  constant.  Source  size  effects  are 

contained  in  the  virtual  origin  term  A , which  is  assumed  to  be  the  same  functior 

y o 

• p 

used  for  the  along-wind  changes  in  c . From  Wilson  (1982),  the  virtual  origin 


Ao  s 


X = 
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_ '/°yo°zo 


0.119 
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(Cl  7) 
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The  ratio  of  reference  height,  zq,  to  source  diameter,  d is 


( C 1 9 ) 


From  the  definition  of  A and  Aq  we  see  that  the  concentration  length  scale  Ac 
is  not  a direct  function  of  the  boundary  layer  thickness  H.  Equation  (Cl 3) 
can  be  written  using  ( C 1 4)  and  ( Cl 7 ) as 


-C7- 


V B2  + /ayoazo)  (C2°) 

with  no  dependence  on  H,  and  no  explicit  dependence  on  turbulent  velocity 

scale  A , A or  A . 
u v w 

Comparison  With  Measurements 

Using  the  data  of  Fackrell  and  Robins  (1982)  from  Figures  Cl  and  C2, 

corrections  were  applied  from  (Cl 2)  to  estimate  the  scale  Acoo  for  each  vertical 

profile  of  Ac-  These  values  are  shown  in  Figure  C3  along  with  the  linear 

correlation  equation.  The  data  suggests  that  B-j  = 2.0  and  82  = 0.34.  This 

figure  shows  that  considerably  more  uncertainty  exists  in  predicting  concentration 

~~r2 

length  scale  than  in  estimating  the  variance  c . The  final  form  for 
estimating  scale  is 


A 

-TT  = °-34(A  + V (C21) 

Using  this  equation  and  ( Cl 2 ) vertical  profiles  of  scale  Ac  were  calculated,  and  th 
predictions  are  compared  with  the  measurements  in  Figures  C4  and  C5.  Consider- 
ing the  scatter  in  the  data,  the  agreement  is  satisfactory. 

Using  the  semi -empi ri cal  theory  of  (Cl 2 and  (C21),  concentration 
fluctuation  length  scales  near  the  ground  can  be  estimated.  However,  an  un- 
certainty of  about  a factor  of  two  may  be  expected  due  to  data  scatter  and  the 
use  of  variance  source  location  and  dissipation  constants  AQ  from  (C18), 

a = 0.6  and  <p  = 0.706.  These  empirical  constants  were  derrived  for  the 

~r2 

variance  field  to  predict  c , and  refined  values  may  be  required  to  optimize 
~rZ 

their  use  for  both  c and  A estimation. 

c 

Finally,  it  should  be  noted  that  the  behaviour  of  the  concentration  and 
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ana  velocity  length  scales  near  the  ground  is  quite  different.  The  high 
intensity  small  scale  eddies  that  are  produced  by  the  mean  flow  wind  shear 
dU/dz  near  the  ground  decrease  the  turbulent  velocity  scales  Au»  Av  and 
as  shown  in  Figure  C6.  At  the  same  time,  these  mean  wind  shear  and 
small  scale  velocity  eddies  are  dissipating  concentration  fluctuations  and 
increasi ng  the  concentration  scale  A . Comparing  Figure  C6  to  Figure  Cl  and 
C2  clearly  shows  the  different  effect  of  the  flow  near  the  ground  on  velocity 
and  concentration  scales. 

So  far  we  have  considered  only  the  spatial  variation  of  length  scale 
Ac  for  some  unspecified  sampling  time  interval.  Because  these  estimates 
use  the  wind  tunnel  data  of  Fackrell  and  Robins, the  sampling  time  is  probably 
about  3 to  30  minutes  in  the  full  scale  atmosphere.  We  expect 
that  as  the  sampling  time  interval  T s increases  that  the  crosswind  meandering 
of  the  plume  will  generate  low  frequency  fluctuations  that  will  increase  the 
concentration  time  scale^..  These  effects  are  discussed  in  Appendix  D. 
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Figure  C3  Along-wind  variation  of  concentration  length  scale  A above 
ground  surface  dissipation  region.  c°° 


Figure  C4  Ground  level  source:  Comparison  of  concentration  length  scales 

with  predictions  for  h/H  = 0,  d = 15  mm,  H = 1 .2  m. 
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Figure  C5  Elevated  source:  Comparison  of  concentration  length  scales 

witn  predictions  for  h/H  = 0.192,  d = 8.5  mm,  H = 1.2  m. 
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Figure  C6  Turbulence  velocity  length  scales  and  mean  velocity  profile 
for  wind  tunnel  boundary  layer  of  Fackrell  and  Robins  (1982) 


APPENDIX  D 


EFFECT  OF  WIND  DIRECTION  MEANDERING 
ON  MEAN  AND  FLUCTUATING  CONCENTRATION 

As  the  sample  time  over  which  concentrations  are  observed  goes  up, 
the  mean  concentration  goes  down  due  to  increased  meandering  caused  by  large 
scale  atmospheric  turbulence.  The  decrease  in  mean  concentration  continues 
even  at  long  sample  times  when  the  plume  has  sampled  all  the  large  scale 
atmospheric  turbulence.  This  long  term  decrease  in  concentration  is 
caused  by  shifts  in  mean  wind  direction  which  force  the  plume  centerline 
to  meander  back  and  forth  over  a fixed  receptor  location. 

Wind  tunnel  data  of  plume  dispersion  and  concentration  fluctuation,  on 

which  the  theory  in  this  report  is  based, is  incapable  of  simulating  the 

plume  centre  line  meandering  caused  by  the  two  effects  of  large  scale 

atmospheric  turbulence  and  changing  mean  wind  direction.  The  wind  tunnel 

data  of  Fackrell  and  Robins  is  limited  by  the  available  turbulence  scales 

to  simulation  of  full  scale  sample  times  3 to  30  minutes.  This 

short  full  scale  sample  time  is  appropriate  for  estimating  the  risk  of 

exposure  to  peak  concentrations  in  short  term  accidental  releases,  such  as 

a pipeline  rupture.  However,  the  long  term  exposure  of  plants  and  animals 

to  pollution  from  a continuous  industrial  source  cannot  be  accurately  estimated 

using  the  theories  developed  from  these  wind  tunnel  simulations.  In  this 

appendix  we  will  develop  a method  for  extrapolating  predictions  from  short 

sample  periods  to  longer  exposure  times  typical  of  continuous  sources.  This 

sampling  time  T$  for  wind  meandering  will  be  used  to  account  for  the  non- 

— 'T. 

stationary  decrease  in  mean  c^  and  variance  c^.  . This  effect  must  be  combined  wit 
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the  variance  smoothing  caused  by  sample  averaging  time  T for  the  stationary 

a 

statistics  which  have  a constant  mean  c for  all  T . 

a 

Accounting  for  Long  Sampling  Time  by  Increased  Crosswind  Spreading 

One  widely  used  method  of  accounting  for  sampling  time  is  to  assume 
that  for  samples  longer  than  three  minutes  the  vertical  spread  remains 
constant,  while  the  crosswind  spread  increases  as  some  empirical  power 
function  of  sample  time.  For  any  fixed  receptor  location  under  an  elevated 
plume  with  an  effective  source  height  h,  ground  level  mean  concentration 
on  the  plume  centre  line  is  given  by 


c 


o 


ttU  o c 
c y z 


exp 


(Dl) 


For  averaging  time  T^,  the  crosswind  spread  is  usually  assumed  to  be 
proportional  to  averaging  time  raised  to  some  power,  p. 


°y.T  = 

cr  r 

y ref 


T 

s 

Ts  ,ref 


P 


032) 


while  the  vertical  spread  oz  remains  constant.  This  model  produces  the 
empirically  observed  result  that  the  mean  concentration  sample  over  T 
behaves  like 


(D3) 


Unfortunately,  this  approach  leads  to  serious  difficulties  when  it  is  applied 
to  concentration  fluctuation  statistics.  The  difficulty  becomes  apparent 
when  the  plume  intermittency  is  calculated.  The  intermi ttency  is  directly 
related  to  the  fluctuation  intensities  by 
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Y = 


i2  + 1 


(D4) 


which  is  valid  for  all  sampling  times.  On  the  plume  axis  both  the 
conditional  (plume)  intensity  i ^ and  the  total  intensity  i of  concentration 
fluctuations  depends  only  on  the  ratio  of  variance  source  strength  to  mean 
source  strengths  and  is  completely  independent  of  a ^ or  a Using  the 
increased  crosswind  spread  model  to  account  for  sampling  time, 
there  should  be  no  effect  of  sampling  time  on  the  observed  plume  intensities 
or  on  the  intermittency  y . This  result  is  clearly  incorrect  because  we 
know  by  observing  long  time  records  of  full  scale  concentration  that  the 
plume  becomes  increasingly  intermittent  with  time,  as  changes  in  mean  wind 
direction  cause  the  plume  to  meander  back  and  forth  across  the  receptor. 
Clearly,  the  use  of  increased  crosswind  spread  is  an  incorrect  model  for 
the  effects  of  sampling  time. 

The  reason  for  the  failure  of  an  increased  crosswind  spread  model  is 
that  as  sampling  time  increases  a plume  will  meander  back  and  forth  allow- 
ing a receptor  to  sample  the  outer  edges  of  the  short  term  plume  fluctua- 
tions. Because  the  concentration  statistics  on  the  plume  axis  carry  no 
information  about  the  structure  in  the  outer  edges  of  the  plume  it  is  not 
possible  to  account  for  increased  sampling  time  simply  by  adjusting  cross- 
wind  spread  and  applying  this  to  values  computed  on  the  axis.  The  outer 
edge  of  the  plume  has  different  turbulence  structure  than  the  plume  axis, 
so  turbulence  properties  such  as  intermittency  y will  change  as  the  receptor 
sees  the  fringes  of  the  plume.  In  the  following  section  we  will  present  a 
more  realistic  model  using  crosswind  averaging  to  allow  the  long  term  time 
average  statistics  to  "sample"  across  the  entire  width  of  the  plume. 


-D4- 


Crosswind  Integration  to  Simulate  Long  Sampling  Times 

The  effect  of  increasing  the  sampling  time  T s over  which  the  mean 

concentration  is  measured  is  equivalent  to  allowing  the  short  sample  time 

plume  to  meander  back  and  forth  over  the  receptor.  This  exposes  the 

receptor  to  the  tails  of  the  crosswind  distributions  of  mean  and  fluctuating 

concentration,  and  changes  both  the  mean  and  fluctuating  statistics  of 

the  plume.  Starting  from  the  short  term  statistics  c , i , i and  y 
r a o o p,  o 

measured  on  the  plume  centerline  for  a short  sample  time  Tsre.p  the  crosswind 
profiles  of  these  parameters  must  be  averaged  to  produce  iQ^  i ^ and  y^ 

for  a longer  sample  time  . The  following  assumptions  will  be  made  to 
model  the  meandering : 

• Correlations  of  effect  of  sampling  time  on  full  scale  data 
use  the  maximum  mean  concentration  cQjObserved  for  a fixed 
sampling  period  T^.  To  model  these  maximum  values  it  will 
be  assumed  that  the  plume  blows  directly  toward  a receptor 
located  on  the  plume  centerline  at  y = 0. 

• The  effect  of  meandering  is  to  expose  the  receptor  to  the 
entire  crosswind  profile  of  the  short  term  plume.  The 
probability  that  the  receptor  is  exposed  to  the  part  of  the 
short  term  plume  located  at  + y off  the  centerline  is  pm  dy, 
where  the  probability  density  p^  is  assumed  to  be  normally 
distri buted . 


where  a is  the  crosswind  meander  standard  deviation, 
m 
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• The  short  term  model,  derived  from  wind  tunnel  measurements, 
represents  a sampling  time  Tsref  that  is  sufficiently  long 
that  increasing  the  sample  time  will  cause  no  further  vertical 
meandering.  Thus,  there  will  be  no  contribution  from  the 
vertical  profile  to  long  term  meandering  effects. 

• The  additional  meander  width  a acts  only  on  the  crosswind 

m 

concentration  profiles.  The  additional  meander  width 
must  increase  with  increasing  sample  time,  and  cr^  -*  0 as 
Ts  Tsref»  where  ^sre^  is  the  sample  period  of  the  short 
term  pi  ume  model  . 

Using  our  assumption  that  the  probability  of  observing  any  part  of  the  short 
term  plume  is  pmdy,  the  crosswind  average  of  any  appropriate  statistic 
represented  by  A(y)  is 


A(y) 


pm(y)  dy 


(D6) 


where  is  the  statistic  for  a sampling  time  Ts , and  A is  the  same  statistic 
for  the  short  term  sample  Tsre^.  Our  crosswind  averaging  model  has  a great 
deal  in  common  with  Gifford's  (1  959),  (1  960)  classic  fluctuating  plume  model 
for  predicting  peak  to  mean  concentration  ratios.  There  are,  however, 
several  fundamental  differences  between  them: 

• Gifford  assumed  an  i nstantaneous  gaussian  plume  profile  with 
vertical  and  crosswind  spreads  az-  and  a ..  for  sampling  time 
T$  = 0.  This  instantaneous  profil e, which  has  no  concentration 
fluctuations,  is  then  allowed  to  meander  over  a long  but 
unspecified  sample  time  to  produce  a total  spread  that  is 
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2 2 

the  sum  of  squares  of  instantaneous  meandering  a^m 

components.  The  meandering  component  a is  zero  at 
T = 0. 


t The  crosswind  averaging  method  assumes  that  there  is  an 
existing  model  for  the  mean  and  fluctuation  statistics, 
and  that  this  existing  model  combines  both  the  instantaneous 
and  meandering  components  for  some  short  sampling  time  Tsref, 
In  the  crosswind  averaging  model  is  the  additional 


meandering  caused  by  sampling  for  times  longer  than  Tsre^. 

Because  crosswind  averaging  relies  on  an  experimentally  based  model  for  the 
short  term  plume,  it  should  provide  a more  accurate  extrapolation  to  longer 
sampling  times  than  the  meandering  plume  model  with  its  hypothetical 
"instantaneous"  profile  which  has  no  internal  turbulent  fluctuations.  Sawfor 
(1983)  in  a Monte-Carlo  simulation  of  puff  dispersion  and  the  data  of 
Fackrell  and  Robins  (1982a),  suggest  that  i and  ip  may  take  on  a finite 
value  far  downwind,  rather  than  going  to  zero  as  the  meandering  component 
becomes  negligible. 

Another  point  that  must  be  clarified  is  the  definition  of  an  "appropriate 
statistic  which  can  be  crosswind  averaged.  The  mean  concentration  c and  the 
intermittency  y (which  represents  the  fraction  of  time  the  plume  is  present) 
are  two  variables  which  can  be  crosswind  averaged.  However,  conditional 
statistics  such  as  cannot  be  obtained  by  direct  crosswind  averaging, 
because  the  fraction  of  time  they  are  present  varies  with  crosswind  position  ' 
Also,  ratios  such  as  the  fluctuation  intensities  i and  i cannot  be  averaged 
because  the  average  ratio  is  not  the  same  as  the  ratio  of  the  average  value 
from  which  it  is  formed.  To  compute  the  average  fluctuation  intensity  iQy 
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it  is  first  necessary  to  determine  the  crosswind  averaged  mean  cQj  and 
*2 

variance  cQj  and  use  these  in  the  definition 


, 1 2 

'oT 


oT  - 


t>T 
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Finally,  we  will  show  later  than  even  the  total  variance  c , for  which  we 
have  an  explicit  equation,  cannot  be  directly  crosswind  averaged  because  each 
local  variance  at  y is  taken  with  respect  to  its  local  mean  c. 


Crosswind  Averaging  of  Mean  Concentration 

The  crosswind  concentration  at  any  fixed  height  z is  given  by  the 
Gaussian  profile  for  the  short  sample  time  Tsre^. 
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where  cq  is  the  concentration  on  the  plume  centerline  at  y = 0.  Integrating 
this  across  the  plume  using  the  meandering  probability  pm  dy  from  (D5) 
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Combining  terms  and  noting  the  function  is  symmetric  about  y = 0 produces 
the  .single  sided  integral 
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Defining  a combined  "spread"  by 
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the  crosswind  average  reduces  to 
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This  integral  is  just  the  cumulative  probability  of  all  y locations  and  by 
definition  is  equal  to  unity.  The  crosswind  averaged  concentration  reduces 
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Crosswind  Averaging  of  Intermi ttency 

Intermi ttency  is  a particularly  important  parameter  because  much  of 
the  effect  of  increased  sampling  time  on  mean  concentration  is  to  reduce  the 
i ntermi ttency  factory  by  introducing  more  periods  of  zero  concentration  as 
meandering  causes  the  receptor  to  see  the  fringes  of  the  crosswind  profile. 
Most  of  the  reduction  of  cQj  with  longer  T$  is  caused  by  decreasi ng  YQy. 

The  conditional  (plume)  concentration  c^defined  by 


°PT  = yT 


(D14) 
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remains  relatively  constant  with  increasing  T^.  Using  the  definition  of 
y from  ( D4)  we  can  write 

2 


1 + i 


Y = Y 


0 1 + i2 


(D15 


where  y is  the  plume  centerline  concentration  for  the  short  term  sample 


time  T 


sref 


1 + i 


(D16) 


For  a single  variance  source  the  source  pair  spacing  is  <p  = 0,  and  the 
fluctuation  intensity  varies  in  the  crosswind  direction  at  any  fixed  height 
z as  a simple  exponential 


.2  . 2 

i = iQ  exp 
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In  this  appendix  only  the  case  of  a single  variance  source  (<|>  - 0)  will  be 
presented.  The  crosswind  averaged  i ntermi ttency  YoTis>  from  (D5)>  (D6)  and 
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Using  the  transformation  c = y/tfw  and  noting  that  the  integrand  is  symmetric 


about  y = 0 leads  to 
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for  variance  source  spacing  0=0 
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This  integral  may  be  easily  evaluated  by  numerical  methods.  Accuracy  to 
four  significant  figures  was  obtained  using  the  trapezoidal  rule  in  50 
equal  steps  from  c = 0 to  c = 5 an  alternative  approach,  using  an 

explicit  approximating  function  is  given  in  (D53). 


Crosswind  Averaging  of  Fluctuation  Variance 


ine  fluctuation  variance  normalized  by  the  plume  centerline  variance 

~r2 

c at  y = 0 is 
o J 


( D20) 


It  is  tempting  to  follow  the  same  procedure  used  for  mean  concentration  and 

integrate  this  over  the  meander  width  a . However,  this  would  not  give  the 

correct  value  for  the  variance  observed  over  the  sample  time  T . The  reasor 

the  variance  cannot  be  averaged  directly  is  because  the  mean  value  c,  about 

which  the  variance  is  measured,  also  varies  across  the  plume.  We  visualize 

the  averaging  process  as  a series  of  discrete  steps  each  of  duration 

Tsref  = 3 min.  to  make  up  a total  sample  time  Then,  each  of  these  three 

1 2 — 

minute  samples  will  have  a different  c and  c,  and  the  average  variance 
over  T^  must  include  the  variance  induced  by  the  changing  mean  c from  step 
to  step. 

To  find  the  crosswind  averaged  variance,  we  will  first  compute  the 

~~2 

averaged  mean  square  concentration  c . Because  this  total  mean  square 

concentration  is  measured  relative  to  zero  it  can  be  crosswind  averaged  in 

the  same  way  as  the  mean  concentration  c was.  For  any  sampling  time  the  mee 

— 1 ? 

square  is  related  to  the  mean  c and  variance  c by 


(D21) 
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In  particular,  for  a sample  time  T during  which  the  plume  meanders 


2 - 2 '2 

oT  " coT  + coT 


( D22) 
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From  this  the  averaged  fluctuation  variance  c^  can  be  determined  once 

2 — 2 — 
c and  c ^ are  known.  Because  we  already  have  c^from  (D13),  it  is  only 

2 

necessary  to  compute  the  mean  square  c j.  Integrating  (D21)  using  (D6) 

( D23) 


oT 


— 2 f '2 
p c dy  + p c dy 


Because  both  integrands  are  symmetric  about  y = 0 they  may  be  written  as 
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4 = 2 pm  c dy 


2 j rm 
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With  the  crosswind  meander  probability  density  from  (D5)  and  c from  (D8) 
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Defining  a*  = a / /2  this  can  be  written  as 
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This  is  exactly  the  same  form  as  (DIO)  with  a*  instead  of  c^.  Carrying  out 
the  transformation  of  (Dll)  and  (D12)  on  (D26)  leads  to  a result  exactly 
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To  evaluate  the  integral  I^,  use  (D20)  for  a single  variance  source  with 
<J>  = 0 to  write 


'2  '2 
c = cQ  exp 
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This  is  exactly  the  same  as  the  form  for  c/cq  so  produces  the  same  average 
as  (D13) 
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Then,  equate  (D22)  and  (D23)  to  show 
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ividing  by  cQj  yields  a total  crosswind  average  fluctuation  intensity 
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Squaring  (D13)  and  using  (D27)  and  (D29)  this  becomes 


( D32) 


for  variance  source  spacing  <p  = 0 
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The  conditional  "plume"  fluctuation  intensity  ipT  may  be  computed  from 
the  exact  relation 


( D33) 


2 

from  the  known  short  term  values  of  y0  and  iQ  » using  (D19)  for  y^and 
(D32)  for  ioT2. 


Influence  of  Sampling  Time  on  Integral  Scales 

As  the  sampling  time  period  increases,  low  frequency  concentration 
fluctuations  are  introduced  by  the  back  and  forth  meandering  of  the  plume 
centerline.  These  fluctuations  will  cause  an  increase  in  the  integral  time 
scale^7cT  observed  over  this  sampling  period.  To  estimate  the  effect  of 
Ts  on  the  observed let  us  look  first  at  the  concentration  spectrum 
E c ( f ) for  a sample  time  Ts>  Pasquill  (1974)  p.  14  shows  that  the  low 
frequency  components  in  the  spectrum  are  attenuated  to  produce  a spectrum 


where  Ecoo,  is  the  spectrum  for  T$  = °°.  The  spectral  density  at  zero  frequency 
is  related  directly  to  the  integral  scale,  see  Appendix  C equation  (C3)  by 


However,  at  f = 0 the  spectral  attenuation  in  (D34)  gives  the  disconcerting 
result  that  E^  (0)  = 0 implying  that ^7cj  z 0 as  well. 


EcT<f>  = EcJf>  1 - 
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A more  useful  approach  is  to  assume  that  both  the  short  sample  time 
spectrum  Ec  for  T$ref  and  the  spectrum  for  sampling  time  T may  be 
represented  as  Markov  spectra,  which  from  Appendix  A,  (A12)  are 
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on  the  centerline  of  the  short  sample  time  plume  and 
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for  the  crosswind  meander  averaged  plume  at  sample  time  T . Then,  the  high 
frequency  components  of  these  two  spectra  must  be  equal,  because  these  high 
frequencies  are  contributed  by  small  scale  turbulence  in  the  short  time 
average  plume  that  meanders  back  and  forth  to  generate  the  longer  sample 
time  statistics.  At  high  frequencies  2TT<^yf  » 1 and  the  spectra  reduces  to 

(D38) 


and 


'2 


Ed(f>  = 


(D39) 


*<7c/ 

These  two  spectra  must  be  the  same,  once  the  effect  of  the  non-stationary 
mean  c-j.  is  accounted  for.  We  will  assume  that  this  may  be  done  by  normal izi 
with  the  square  of  the  mean  so  that 
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at  high  frequencies.  Using  (D38)  and  (D39)  yields,  for  sampling  time  Ts 
and  reference  sample  time  Tsref  on  the  plume  centerline 


<J  CO  <7coT 


(D41  ) 


where  iQ  and^7cQ  are  the  intensity  and  scale  on  the  plume  centerline  for  a 
sample  time  Tsref  and  i0jand^coT  are  values  for  sample  time  during 
which  the  short  term  plume  is  crosswind  averaged  by  meandering. 


Relating  Meandering  Width  to  Sampling  Time 

All  the  analysis  developed  so  far  can  only  be  useful  if  the  meander  width 
ratio  can  de  directly  related  to  the  sampling  time  . We  will  begin 

by  assuming  a power  law  relationship  between  the  crosswind  averaged  mean 
concentration,  and  the  sampling  time  in  the  form, from  (D3) 


(D42 ) 


Pasqui 1 1 (1974)  p.  65  carries  out  an  analysis  of  the  effect  of  sampling 
time  on  the  crosswind  turbulent  velocity  spectrum  that  contributes  to  a j 
in  (D2).  For  short  sampling  times  where  the  attentuation  factor  in  (D34) 
cuts  the  spectrum  off  in  the  inertial  subrange  he  shows  that  p = 1/3.  For 
longer  sampling  times  p will  gradually  decrease,  unti  1 when  Ts+  °°,  p ->  0. 
For  a Markov  spectrum  a value  of  p = 0.5  is  the  upper  limit  for  short 
sampling  times. 


The  spectral  attenuation  caused  by  averaging  time  T is  just  the 

a 

complementary  value  of  the  sampling  time  attenuation-  compare  Appendix  A, 
(A27)  with  (D34).  For  this  reason,  vari ati on  of  maximum  mean  concentration 
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with  averaging  time  T will  have  the  same  exponent  p for  a range  of  averaging 

a 

and  sampling  times  where  they  are  acting  on  the  same  part  of  the  turbulence 
spectrum,  with  averaging  time  cutting  off  the  high  frequencies  and  sampling 
time  tnelow  frequencies.  Time  averaged  air  quality  data  analyzed  by 
McGuire  and  Noll  (1971)  gave  values  of  p from  0.175  to  0.369  with  a mean 
of  p = 0.197  for  urban  S0£  concentrations.  The  power  law  fit  their  data 
for  averaging  times  from  a few  minutes  to  a year.  Angle  (1978)  on  the 
basis  of  other  data  recommended  p = 0.20  for  averaging  times  from  3 min. 
to  several  hours. 

Most  of  the  data  for  the  effect  of  sampling  time  on  mean  concentration 
is  for  urban  data  with  many  overlapping  plumes.  The  effect  of  multiple 
sources  ranged  around  the  receptor  causes  mean  concentrations  to  decrease 
more  slowly  than  would  be  the  case  for  a single  source.  In  addition,  for 
sampling  times  longer  than  about  an  hour  local  terrain  and  drainage  flows 
will  strongly  influence  the  meandering  of  wind  direction,  even  for  large 
ensemble  averages.  These  factors  in  mind,  the  value  of  the  exponent  p is 
expected  to  lie  somewhere  between  0.1  and  0.5.  Angle's  estimate  of  p = 0.2 
may  be  conservative,  because  it  will  produce  higher  dosage  and  probabilities 
of  theshold  exposure  by  producing  a slower  decrease  in  Y0-j-and  c^ythan  would 
occur  for  a larger  exponent. 


The  relationship  between  meander  width  g^  and  the  sampling  time  is 
found  by  simply  equating  the  empi ri cal  power  law  in  (D42)  to  the  analytical 
prediction  of  (D13)  to  obtain 
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( D43 ) 
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which  reduces  to 
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Crosswind  Profiles  at  Long  Sample  Times 

So  far,  we  have  used  the  crosswind  profiles  of  plume  statistics  at  a 
short  sample  time  Tsre^  to  estimate  the  statistics  on  the  centerl ine  at 
longer  sample  times.  Once  these  centerline  statistics  cq^.,  i^,  i^  are 
known  it  is  a simple  matter  to  compute  off-axis  values  at  time  T . It  seems 
reasonable  to  assume  that  the  functional  forms  for  the  crosswind  profiles 
are  the  same  at  all  sampling  times,  with  the  only  change  being  the  width 
a y of  the  crosswind  spread.  With  this  assumption,  the  off  axis  functions 
for  a single  variance  source  (<j>  = 0)  are: 
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( D47 ) 


(D48) 


The  intermittency  at  any  off-axis  location  is  found  from  the  definition 
. 2 


pT  + 1 

yt  = T2 — r 

1T  + 1 

using  (D45)  and  (D46). 


(D49) 
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Closed  Form  Approximation  to  Crosswind  Averaged  Intermi ttency 

The  intermi  ttency  on  the  plume  centerline  y0-j-  for  a sampling  time  T$ 
is  related  to  the  short  term  centerline  i ntermi ttency  yQ  (for  sample  time 
Tsref)  by  (D19).  The  basic  assumption  embedded  in  this  equation  is  a 
Gaussian  probability  density  for  the  meandering  spread  a . This  meander 
spread  am  is  zero  at  the  short  sampling  time  of  Ts^e^  and  increases  with 
sampling  time  according  to  (D44). 

A closed  form  solution  to  the  integral  in  (D19)  could  not  be  found. 

This  is  unfortunate,  because  with  this  one  exception  all  the  other  equations 
used  for  conditional  and  time  dependent  variables  are  in  an  explicit  closed 
form.  The  lack  of  a closed  form  for  (D19)  is  one  of  the  major  obstacles 
to  developing  rapid  microcomputer  based  models  for  hazard  assessment.  In 
the  following  analysis  an  approximate  equation  for  the  integral  in  (D19) 
will  be  developed  to  avoid  the  time-consuming  numerical  integration. 

In  order  to  develop  an  approximate  form  we  must  first  look  at  the 
asymptotic  limits  of  the  integral  in  (D19)  for  large  and  small  values  of 
the  plume  centerline  concentration  fluctuation  intensity  i . For  very 
small  intensity  with  i « 1 we  find,  setting  i = 0 in  (D19)  that 

— = 1 .0  ( D50) 

Yo 

for  i =0 
o 

This  is  the  expected  result,  because,  when  the  plume  is  spreading  by 
diffusion  in  very  small  scale  turbulence  the  centerline  intensity  iQ  will  be 
small  and  there  will  be  little  meandering  so  that  y0  = yqj  = 1-0. 

At  the  other  extreme,  when  the  concentration  fluctuation  intensity  iQ  i 
large,  an  exact  solution  of  (D19)  is  possible  if 
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( D51 ) 


so  that  the  "1"  in  the  denominator  of  (D19)  can  be  neglected.  The  condition 
in  (D51)  is  quite  easy  to  meet,  because  in  the  fringes  of  the  plume  where 
the  dummy  variable  £ is  large,  the  exponential  term  in  (D51)  makes  the 
inequality  true  even  for  values  of  i greater  than  unity.  Neglecting  the  "1" 
in  the  denominator  of  (D19)  leads  to  the  exact  asymptotic  solution 


which,  surprisingly,  is  independent  of  the  concentration  fluctuation  intensity,  i . 

To  develop  an  approximating  equation  that  fits  smoothly  and  accurately 
between  the  asymptotic  limits  of  (D50)  and  (D52)  requires  a large  measure 
of  educated  guesswork.  Using  some  suggestions  for  J.W.  Rottman  (1984, 
private  communication)  the  integrand  in  (D19)  was  transformed  from  an  exponential 
to  a logarithmic  form,  and  then  approximated  by  a power  function  and  integrated. 
After  considerable  trial  and  error  the  best  approximate  form  was  found  to  be 
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for  "large"  i 
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The  error  in  using  (D53)  with  (D54)  and  (D55)  as  an  approximation  to  (D19) 
was  found  to  depend  on  the  value  of  the  centerline  intensity  iQ.  The 
approximate  value  of  the  absolute  error  (as  opposed  to  the  fractional  or 
percent  errors)  in  Y0-j-/y0  is 


error  < 


+ 0.02 


(056) 


For  values  of  iQ  greater  than  0.05.  At  very  small  values  of  iQ  - 0.01  the 
absolute  error  may  be  up  to  twice  as  large  as  predicted  by  (056). 


Relating  Meandering  of  the  Instantaneous  Plume  to  Additional 
Meandering  a 

The  actual  meandering  component  plume  spread  is  usually  applied  to 
the  ensemble  averaged  instantaneous  spread  which  occurs  for  zero 
sampling  time  when  Ts  = 0.  Gifford's  (1959)  meandering  plume  model  uses 
this  definition. 

The  analysis  in  this  report  takes  a different  approach.  Because  the 
concentration  fluctuation  statistics  are  only  known  for  the  total  (instantaneo 
plus  meandering)  plume  at  some  reference  time  Tsre^,  the  meandering 
component  am  is  the  additional  meandering  experienced  by  the  plume  for 
sampling  times  T$  that  are  longer  than  Tsre^.  While  it  would  be  better  to 
base  the  effects  of  meandering  on  the  instantaneous  plume,  our  ignorance 
of  its  concentration  fluctuation  statistics  prevents  this. 
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How  is  the  addi tional  plume  meandering  am  related  to  the  instantaneous 
plume  meandering  a^?  For  the  plume  at  sampling  time  T .p,  the  mean 
concentration  on  the  plume  axis  is 


Co  2TiUa^az 

For  sampling  time  Ts  the  axis  concentration  is 


( D57 ) 
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( D58) 


Vertical  meandering  is  generally  less  than  crosswind  meandering.  Because 
vertical  meandering  involves  turbulence  scales  much  smaller  than  horizontal 
scales,  Pasquill  (1974)  suggests  that  there  will  be  no  additional  vertical 
meandering  for  sample  times  longer  than  about  three  minutes.  By  requiring 
Tsref  > 3 min.  we  avoid  the  need  to  define  addi tional  vertical  meandering 
and  ozj  = gz-  Then,  taking  the  ratio  of  the  above  two  equations 


C -r  a 
oT  _ y 

c a _ 
o yT 

Equating  this  to  the  (D13)  leads  to 

2 2 x 2 
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( D59 ) 


( D60) 


This  sum  of  variances  is  the  expected  result  for  the  random  uncorrelated 
variables  of  plume  spread  a ^ at  ^sref>  and  The  addi tional  crosswind  meandering 
a . To  relate  the  additional  meandering  to  the  instantaneous  plume  meandering, 
a_jm,  we  note  that  the  total  spread  o^  at  ^sref  must  be 


2 2 

0=0-  + O.  n 

y ly  im,ref 

where  a.  .is  the  instantaneous  plume  meandering  at  sample  time  T 
im,ref  J 

Also,  by  definition,  at  sample  time  Ts , 


(D61) 
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Using  (D61)  in  (D60)  and  then  comparing  the  result  with  (D62)  leads  to 
the  simple  relationship 
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aim,ref 


( D63) 
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That  is,  the  total  meandering  spread  a-  is  found  by  adding  the  squares 

2 

of  the  instantaneous  plume  meanderinq  a.  r at  T and  the  additional 

r 3 l m , re  f sref  — ■ - — 

2 

meandering  am  that  occurs  as  the  sampling  time  is  increased  to  T . 

If  a model  for  the  concentration  fluctuation  statistics  can  be 
developed  for  the  instantaneous  plume,  it  may  be  more  useful  to  express 
the  effects  of  sampling  time  on  the  instantaneous  Ts  = 0 plume  than  on  the 
rather  arbitrary  reference  sampling  time  ^sref  However,  there  are 
difficulties  with  this  approach.  For  example,  the  mean  concentration  on 
the  centerline  of  the  instantaneous  plume  is 


C°i  = 2"Uaiyaiz 


( D64) 


Taking  the  ratio  of  (D58)  to  (D64)  yields 
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( D65) 


from  which  we  see  that  meandering  in  both  the  vertical  and  crosswind 

directions  must  be  considered  where  the  instantaneous  plume  is  used  as  a 

reference.  This  isn't  a serious  drawback  because  for  T$  > 3 min.  the 

ratio  of  a T/a • becomes  constant. 
zT  i z 


APPENDIX  E 


DOSAGE  FLUCTUATIONS 


For  a biological  receptor  the  integrated  concentration  exposure  is  often 
a more  important  factor  in  determining  damage  than  the  peak  concentration. 

For  a receptor  that  is  exposed  in  the  time  period  from  t = t-j  to  t = t-j  + T 
the  exposure  dose  is  defined  as 


D = 


VTe 

c dt 


(El) 


The  average  concentration  over  the  same  time  period  is 


tl  +Te 


c dt 


(E2) 


From  which  we  see  that  the  dose  is  related  to  the  exposure  time  by  the  simple 
equation 


(E3) 


For  risk  analysis,  receptors  which  appear  during  the  release  period  are  never 
considered,  and  all  receptors  are  assumed  to  be  exposed  to  the  entire  transient 
release.  Likewise,  in  a continuous  plume  the  time  origin  is  taken  at  the  instant 
that  the  receptor  is  first  exposed.  Thus,  in  both  cases  t-j  e 0 and  the  dose 
equations  simplify  for  fluctuating  concentration  c = c + c'  to 
T 
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D = 


(c  + c' ) dt 


(E4) 
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and 
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If  we  take  an  ensemble  average  of  many  exposures  to  find  the  mean  dosage  D 
caused  by  a fluctuating  concentration,  it  is  easy  to  show  that  the  fluctuating 
component  c'  cancels  in  the  ensemble  average  to  give 


(E6) 


o 

where  the  overbar  denotes  an  ensemble  rather  than  a time  average.  In  the  same 
way, 


c T 


and  using  (E3) 


(E7) 


c T Te 


(E8) 


where  c-p  is  the  ensemble  average  concentration  over  a sample  time  equal  to 
the  exposure  time  T . Comparing  (E3)  and  (E8),  D/cT  = D/Cj 

Dosage  Fluctuations  Caused  by  Concentration  Variance 

Defining  a mean  and  fluctuating  dosage  as  D and  D‘ , the  dosage  in  a 
single  realization  of  an  ensemble  is 


U = D + D'  (E9) 

Taking  the  mean  square  and  carrying  out  an  ensemble  average  (denoted  by  an 
overbar) 


(E10) 
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Using  the  mean  dose  D from  (E6),  the  dosage  fluctuation  can  be  found  from 

~2 

( El 0)  by  calculating  the  mean  square  D . Beginning  with  the  square  of  the  dose, 
and  using  two  time  variables  t and  t as  a reminder  that  one  may  vary  while 
the  other  is  held  constant. 


D 


2 


(Ell) 


Writing  c = c + c'  in  each  of  the  above  integrals,  and  expanding  by  multiplication 
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Take  an  ensemble  average,  noting  that  the  last  two  terms  will  vanish  because  the 
ensemble  average  of  c'  is  zero,  and  the  integrals  of  c and  c'  are  uncorrel ated. 
Recognizing  the  first  term  as  the  square  of  the  mean  dosage  reduces  the  average  to 


(t)dt 


(El  3) 
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We  may  relate  the  two  time  variables  t and  t by  defining  a new  variable  t' 
which  is  t = t + t'.  By  direct  comparison  with  (E10),  and  using  t1,  the  dosage 
variance  is,  for  t held  constant  while  t'  varies 


Te  Te 
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c' (t)c' (t  + t' ) dt  dt ' 
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This  is  related  to  the  fluctuation  variance  c with  an  averaging  time  T 

a a 
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equal  to  the  dosage  exposure  time  Te.  By  direct  comparison  with  equation 

*2 

(A36)  in  Appendix  A we  see  from  (El 4),  defining  cg  when  = Te 


(El  5) 


'2. 


where,  from  Appendix  A,  (A30)  relates  the  time  averaged  variance  to  the 
~rZ 

true  variance  cT  for  a sampling  time  duration  T and  averaging  time  T equal  to 
I s a 

the  dose  exposure  time  T 
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where^^is  the  integral  time  scale  of  concentration  fluctuations.  For  exposure 

times  T that  are  long  compared  to  the  concentration  time  scale  (A31)  shows 
e c r 

that  the  variance  ratio  reduces  to 
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Values  of  T a47c-j>  5 may  be  taken  as  "large".  Combining  (El 5)  and  (El  7)  the  dosa 
variance  is 
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(El  8) 


Taking  the  ratio  of  the  square  root  of  (E18)  to  (E8)  gives  the  dosage  fluctuatioi 
intensity  i‘D 

(E19) 


-E5- 


* 2 ^ 

The  variance  and  scale  j may  be  related  to  the  values  in  the  short  term 
plume  by  Appendix  D,  (D42) 


^cT  <7co 


(E20) 


where  i and<^7  are  the  concentration  fluctuation  intensities  for  the  short 
o Vco 

sampling  time  Tcre^  appropriate  for  the  wind  tunnel  simulation.  Using  (E8) 
and  (El 8)  in  (E19)  , 


(E21) 


Then,  using  (E20) 
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The  ratio  i-j-  /iQ  can  be  computed  from  Appendix  D,  (D32)  while  the  short  term 
time  scale  is  evaluated  in  Appendix  C for  Tsre^  = 3 min. 


APPENDIX  F 


PROBABILITY  OF  THRESHOLD  LEVEL 
CROSSING  DURING  EXPOSURE 

For  pollutants  whose  effects  are  cumulative,  a receptor  responds  to 
the  total  integrated  dose  to  which  the  receptor  is  exposed  during  a time 
interval  T . At  the  other  extreme,  exposure  to  highly  toxic  gases  such  as 
hydrogen  sulphide  may  be  best  described  by  a threshold  concentration  level,  c* 
below  which  there  is  no  response  and  above  which  any  short  exposure  will  be 
fatal.  Here,  estimating  the  chance  of  not  once  exceeding  c*  during  the 
time  T£  is  complicated  by  the  need  to  determine  the  probability  of  getting 
through  the  entire  exposure  without  ever  crossing  the  threshold  concentration. 

The  statistical  information  that  is  currently  available  only  allows  us  to  compute 
the  fraction  of  time,  V,  that  the  concentration  will  be  more  than  the  threshold 
for  a large  ensemble  average  of  exposures,  and  not  the  probability  of  a first 
threshold  crossing  occurring  during  an  ensemble  of  exposures.  Some  assumptions 
will  be  necessary  to  allow  "fraction  of  time"  probabilities  to  be  used  to  predict 
time  dependent  "first  crossing"  probabilities. 

Fraction  of  Time  Over  a Threshold 

Before  dealing  with  the  "first  crossing"  problem  it  is  useful  to  review  the 
statistical  analysis  in  current  use.  The  only  way  in  which  concentration  time 
series  are  now  analyzed  is  to  sort  the  time  series  measurements  into  concentration 
bins  to  specify  probability  density  p(c)dc  as  the  fraction  of  time  the  concentra- 
tion falls  in  the  bin  from  c to  c + dc.  Periods  of  zero  concentration  are 
usually  dealt  with  separately  through  the  i ntermi ttency  factory  and  a condition- 
al "plume"  probability  density  p (c)  which  describes  all  concentrations 
greater  than  zero, 

p(c)  = (1  - y)  6 (c)  + Y Pp(c) 


(FI) 


-F2- 


where  (1  - y)  is  the  fraction  of  time  with  zero  concentration  and  6(c)  is 
the  delta  function,  which  is  a spike  at  c = 0 with  the  properties  that 
6(c)  = 0 except  at  c = 0 and  an  area  under  6(o)  of  unity.  The  fraction  of 
time  that  the  concentration  is  less  than  some  threshold  c*,  is  given  by  the 
cumulative  distribution  function  ft(c*) . 


(F2) 


c* 


« (c*)  = (1  -y)  + Y 

0 


Pp(c)  dc 


(F3) 


The  fraction  of  time  during  which  the  threshold  c*  is  exceeded  is  V(c*),  and 
is  simply  the  complementary  cumulative  distribution  of  ft(c*) 

v(c*)  = l - a (c*)  (F4) 

Threshold  Crossing  Probability  for  a Series  of  Short  Exposures 

The  probability  V(c*)  may  be  thought  of  as  the  probability  of  exceeding 
the  threshold  c*  at  any  instant  during  an  exposure  time  JQ.  To  estimate  the 
probability  of  crossing  this  threshold  for  the  fi rst  time  we  must  specify 
the  hazard  rate  r.  This  is  the  probability  per  unit  time  of  exceeding  the 
threshold  c*  for  the  first  time,  and  this  "first  time  of  passage"  condition 
requires  that  we  have  always  had  c less  than  c*.  The  probability  of 
exceeding  the  threshold  for  a short  exposure  of  length  At  for  the  first 
time  is  just  rAt,  and  the  probability  of  not  crossing  the  threshold  is 
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Si  = (1  - rAt)  ( 

To  calculate  the  probability  of  not  crossing  for  N successive  exposures  each  of 
equal  duration  At,  two  assumptions  are  required. 

1)  The  hazard  rate  r for  for  a unit  time  of  exposure  is  equal 
to  the  probability  V of  exceeding  the  threshold  for  an 
instant . 

2)  The  probability  that  the  short  exposure  exceeds  the  threshold 
is  independent  of  the  probabilities  of  exceeding  the  threshold 

in  any  previous  interval. 

The  first  assumption  implies  that  the  time  interval  At  of  each  exposure  is  short 
enough  that  the  probability  of  exceeding  c*  remains  constant  over  the  time 
from  t to  t + At. 

The  second  assumption;  that  successive  exposures  are  statistically  indepen- 
dent, requires  more  justification.  We  know  that  the  concentration  time  series 
is  strongly  autocorrel ated,  with  a correlation  coefficient  c 1 (t)  c ' ( t + t) 
having  an  integral  time  scale <7c  of  30  to  100  seconds.  It  would  appear  from 
this  highly  correlated  time  series  that  the  probability  of  observing  a 
concentration  greater  than  c*  in  the  j^  exposure  should  be  strongly  dependent 
on  the  value  of  c in  the  previous  exposure  period  (j  -1).  In  fact  this  is 
true,  and  if  we  knew  the  value  of  the  concentration  during  the  (j  - 1)  exposure 
interval,  the  assumption  of  i ndependence woul d be  incorrect.  However,  all  we 
know  from  the  previous  exposure  interval  is  that  c lies  somewhere  in  range 
0 £C  < c*,  not  the  specific  value  of  c.  As  long  as  the  threshold  concentra- 
tion is  large  enough  for  a wide  range  of  possible  c values  to  occur  in  the 
previous  interval  then  assumption  #2  is  reasonable.  As  a rough  estimate, 
successive  exposure  intervals  should  be  statistically  independent  as  long  as 
c*  is  larger  than  the  mean  concentration  c during  an  interval.  For  most  risk 
analysis  this  condition  is  easily  met. 
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Using  these  assumptions,  the  hazard  rate  for  any  exposure  is  just 


(F6) 


and  for  N statistically  independent  exposures  the  probability  of  exposure 
without  once  crossing  the  threshold  is  just  the  product  of  the  series 
of  probabilities  of  not  crossing  during  each  At 


N 

s = n (i  -v  At) 

j=i  J 


(F7) 


where  n indicates  the  product  of  N terms.  Often  the  time  interval  is  expressed 
in  terms  of  the  respiration  frequency  of  an  exposed  person.  This  breathing 
rate  B expressed  in  "breaths  per  second"  defines  the  exposure  time  interval  as 


(F8) 


For  a series  of  N breaths,  combining  (F7)  and  (F8) 


N V. 

: n (1-  i) 

J=1 


(F9) 


Steady  Release  With  Long  Exposure  Time 

For  a continuous  steady  source,  the  plume  mean  concentration  c , 
fluctuation  intensity  i and  intermittency  y all  remain  constant  with  time. 

For  a given  threshold  level  c*,  let  the  exceedance  probability  V-  be  a constant, 
V,  for  each  exposure  time  interval.  For  N breaths  the  probability  (F9)  is 


(F10) 
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In  terms  of  the  total  exposure  time  Te,  the  number  of  breaths  N depends  on 
the  breathing  rate  B 


N = B T 

e 

so  that  in  terms  of  exposure  time  we  can  write  (F10)  as 


S(T 


N 


( FI  1 ) 


( FI  2) 


For  long  exposure  times,  the  number  of  breaths  becomes  large  and  we  use  the 
fundamental  definition  of  an  exponential  function 


lim  (1  -77-)  = exp  ( -x) 

N ->oo  IN 

to  write  the  probability  S of  not  once  exceeding  c*  as 


S(Te)  = exp(-V  Te) 
for  large  N 


( FI  3) 


( FI  4) 


Note  that  the  survival  probability  S is  independent  of  the  breathing  rate  B 
for  sufficiently  long  exposure  times.  Assuming  successive  time  intervals 
are  statistically  independent  for  exposure  probability,  the  exposure  time  does 
not  need  to  be  continuous  but  may  be  the  sum  of  a series  of  short  exposure 
times  (such  as  might  occur  if  a subject  was  entering  and  leaving  a contaminated 
area  several  times). 


Generalized  Model  for  Continuous  Exposure 

Conceptually  it  is  easier  to  visualize  the  exposure  as  a series  of  dis- 
crete events  rather  than  a continuous  process.  The  advantage  of  regarding  the 
exposure  as  a continuous  process  is  that  a time  varying  hazard  rate  that 
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occurs  for  a transient  release  can  be  easily  accounted  for.  If  we  define  the 
following  continuous  probability  functions: 

f(t)dt  = the  probability  that  the  time  to  the  first  passage  over 

the  threshold  will  occur  in  the  interval  from  t to  t + dt 

f(t)  = the  probability  density  function  for  the  first  passage 

★ 

time  over  the  threshold  c 

F(t)  = the  probability  that  the  threshold  c*  has  been  exceeded 
at  least  once  before  time  t. 

S(t)  = the  probability  of  being  exposed  to  c from  t = 0 to  t 
without  exceeding  the  threshold  c*. 

r(t)dt  = the  probability  of  exceeding  the  threshold  at  time  t given 
that  the  first  passage  over  the  threshold  c*  has  not  yet 
occurred. 

These  probabilities  can  be  related  using  the  hazard  rate  r,  which  we  now  allow 
to  be  a time  varying  function  r(t).  We  know  that  the  probability  f(t)dt  of 
fi rst  crossing  the  threshold  is  simply  the  product  of  the  probability  S(t)  of 
being  exposed  from  t = 0 to  t wi thout  crossing  the  threshold  and  time 
probability  r(t)dt  of  exceeding  the  threshold  in  the  current  time  step  from 
t to  t + dt.  That  is 

f(t)  dt  = S ( t ) r( t ) dt  (F15) 

By  definition  S(t)  is  just  the  complementary  probability  of  F(t) 


S(t)  = 1 - F(t) 


(F16) 
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and  F(t)  must  be  the  time  cumulative  probability  distribution  of  the  time 
probability  density  f(t).  That  is. 


t 


F(t) 


f (t 1 ) d t ' 


( FI  7) 


o 

where  t'  is  a dummy  time  variable  of  integration.  By  differentiating  (FI 7)  we 
find  the  usual  relation  between  the  cdf  and  pdf 


dF(t) 

dt 


- f(t) 


(FI  8) 


Then,  combine  (15)  and  (16)  to  obtain 


f(t) 


F(t)j  r(t) 


(FI  9) 


Equate  (18)  and  (19)  and  transpose  terms  to  obtain  a differential  equation  for 
the  cdf 


dF(t 


1 - Hi) 


r(t)dt 


(F20) 


We  then  use  our  first  assumption  that  the  hazard  rate  r(t)  is  equal  to  the 
instantaneous  probability  of  exceedance  V(t)  and  integrate  (F20)  over  an 
exposure  time  from  t = 0 to  t = T 


dF 


- F 


V ( t ) dt 


(F21) 


From  which 


T T 
e e 


-£n  (1  - F(t) ) 


V ( t ) d t 


( F22) 
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Because  F(0)  = 0 this  becomes 

T 

e 

£n(l  - F(Te))  = - V ( t ) d t ( F23) 

o 

Using  the  definition  of  the  survival  probability  in  (FI 6),  the  final  result  is 

(F24) 


' 

~ 

e 

S(Te)  = exp 

- 

V ( t ) dt 

L c 

) J 

which  shows  that  all  that  is  required  to  estimate  the  survival  probability  for 
an  exposure  time  Te  is  the  integrated  probability  of  exceedance  over  the  time 
interval.  For  the  case  where  V(t)  is  constant  in  a steady  release,  (F24) 
reduces  immediately  to  (FI 4).  If  we  define  an  average  probability  of  exceedance 
over  an  exposure  time  Te  as 


_ 1 


avg  T 


(F25) 


V ( t ) dt 

The  survival  probability  S is  then 

(F26) 

which  simplifies  the  task  of  estimating  risk  of  exposure  to  time  varying 


S(Te)  = exp 


- V T 
avg  e 


releases. 


APPENDIX  G 


STOCHASTIC  MODEL  OF  GAS  CONCENTRATIONS 

I . Introduction 

A mathematical  model  will  be  proposed  which  will  predict  concentrations 
of  a target  gas  released  as  a result  of  either  a pipeline  rupture  or  a 
continuous  plume.  The  model,  when  provided  with  estimates  of  several  basic 
parameters,  will  be  used  to  specify; 

a)  the  steady  state  probability  distribution  of  gas  concentration  p(c), 
and 

b)  the  probability  F(Te)  that  the  passage  to  a concentration  higher 
than  a specified  threshold  will  occur  at  least  once  during  a given 
exposure  time  Te. 

The  perspective  taken  is  that  the  initial  gas  concentration  on  exit 
from  the  source  is  cQ,  the  gas  exits  in  a continuous  process  but  is  quickly 
transformed  into  a series  of  discrete  eddies,  and  each  eddy  undergoes  a 
dilution  process  as  a result  of  air  entrainment  and  diffusion.  During  the 
dilution  process  and  as  function  of  the  time  t since  release,  the  gas 
concentration  of  each  eddy  is  a random  variable  and  the  average 
concentration  decreases. 

The  analysis  will  be  carried  out  as  follows: 

a)  In  section  II  the  dilution  process  will  be  modeled  as  a continuous 
random  walk  Wiener  process.  It  will  be  shown  that  if  the  gas 
concentrations  follow  this  type  of  process  then  the  probability 
density  function  of  gas  concentration  in  the  shcrt  run  will  be 
normal  (Gaussian)  and  in  the  long  run  an  exponential.  The  Gamma 
probability  density  function  is  suggested  as  an  approximation  to  the 
steady  state  probability  distribution  of  gas  concentration. 
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b)  In  section  III  it  is  noted  that  the  plume  over  a sample  site  is 


not  continuous  and  because  of  meandering  and  holes  in  the  plume 


there  are  intervals  of  time  that  the  plume  is  not  at  the  sample 
site.  The  meandering  process  is  modeled  as  a two  stage  Markov 
process  and  the  characteristics  of  the  model  are  highlighted. 

c)  In  section  IV  the  time  to  first  passage  past  a concentration 
threshold  will  be  discussed  and  an  expression  developed  for  the 
probability  that  a concentration  threshold  will  be  exceeded,  as  a 
function  of  time. 

d)  In  section  V numerical  routines  will  be  given  which  are  required  for 
calculation  of  the  relevant  variables  of  the  model. 

II.  The  Dilution  Process 

The  Wiener  process  has  been  used  as  a model  for  many  physical 
processes,  including  the  diffusion  of  one  gas  into  another.  The  Wiener 
process  may  be  regarded  as  a limiting  form  of  the  random  walk  in  one 
dimension.  Consequently,  no  discrete  jumps  occur  and  in  a small  interval  of 
time  the  process  is  restricted  to  local  transitions  and  has  independent 
increments.  The  displacement  per  interval  A t is  of  the  order  of  magnitude 
Atl/2f  and  in  any  single  realization  of  this  continuous  process  one  would 
expect  to  find  an  infinite  number  of  small  spikes  in  any  finite  time 
interval . 

The  equations  describing  the  gas  concentration  at  time  t from  release 
from  the  source  Cp(t),  for  a Wiener  process  are  also  a solution  of  the 
Fokker-Planck  equation  for  diffusion. 


G1 


= constant,  and  cQp  = variance  defined  by  the  process* 


where 
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This  one-dimension  diffusion  equation  with  an  external  field  of  force,  such 
as  gravity,  is  normally  used  to  calculate  the  density  of  particles  at  a 
point  x at  time  t given  the  initial  position  xQ,  p(x,t|xQ).  For  the 
dilution  process  the  point  is  fixed  to  the  frame  of  reference,  an  eddy  in 
the  plume,  and  the  distribution  for  the  concentration  within  the  eddy  is 
required.  Given  the  initial  condition  Cp(0)=co,  and  unrestricted  values  for 
Cp(t);  — 00  <Cp(  t ) < oo  , the  probability  density  function  of  gas  concentration 
will  be  Gaussian  with: 


c 

P 


c 

0 


dt 


t » where 


dCp(t) 

dt 


is  sometimes  called  the  drift 


G2 


and 


For  notational  convenience,  define  P 

dt 


dcP(o 

dt 


and 


G3 


For  the  dilution  process,  however,  the  gas  concentration  in  the  eddy  is 
restricted  to  be  greater  than  zero;  Cp(t)>0.  The  appropriate  boundary 
equation  for  this  case  is 


1^7 

2 cop 


d p (y* 1 co) 

dc 

P 


= 0 


= 0 
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It  is  known  that  the  solution  to  the  resulting  diffusion  equation  is  (e.g., 
Cox  and  Miller,  [1965],  p.224); 
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p(cp’t  l c0) 
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dc 
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op 
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where  $ (z)  is  the  cumulative  density  function  for  the  standardized  normal 

variate  Z.  The  equilibrium  or  steady  state  distribution  exists  only  for 
dc. 


[P  < 0 , and 

at 


P(C 


|co} 


dc 

dc  I 

_R 

( 2 

_J l\ 

dt 

...» 

dt  1 

) c > 0 

p(  ; p 


G6 


‘2 


op  op 

which  is  the  exponential  probability  density  function. 

If  examined  closely  it  becomes  apparent  that  p(cp,t|cG)  is  a 

combination  of  an  exponential,  a normal  and  another  distribution  which  is 

normal-like.  For  the  dilution  process,  the  drift  is  expected  to  be 

negative  since,  due  to  air  entrainment  and  diffusion,  the  tendency  is 

for  the  average  gas  concentration  to  decrease  as  time  since  emission  from 

the  source  t,  increases. 

dc 

When  the  product  P t is  small  the  distribution  of  gas 

d t 

dc 

concentration  is  predicted  to  be  approximately  Gaussian,  When  — t is 
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large  the  distribution  of  gas  concentration  approaches  an  exponential 

probability  density  function.  An  example  of  the  changing  shape  of  the  gas 

concentration  distribution  as  a function  of  time  as  predicted  by  the  model 

dc 

is  shown  in  Figure  G.l.  Also  note  that  E is  assumed  to  be  a constant 

dt 

and  is  independent  of  the  initial  condition  c0,  but  it  is  dependent  on  the 
process  being  analyzed.  Cp(0)=co  will  be  higher  than  the  final  value  of 


cd(  00  ) > since 


cp(“)  = 2 — 

Cop 
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— E may  be  either  positive,  zero,  or  negative,  but  for  the 

dt  

dc 

dilution  process  studied  here  cQ  is  expected  to  be  relatively  high,  E 

dt 

to  be  negative  and  the  average  concentration  will  decrease  to  a final  value 
of  cp(  oo  ). 

The  principal  predictions  of  the  model  are: 

a)  Immediately  after  release,  and  fcr  some  time  period  from  release, 
the  distribution  of  gas  concentrations  is  approximately  Gaussian. 

b)  When  time  since  exit  from  the  source  is  large  the  distribution  of 
gas  concentrations  will  follow  the  exponential  probability  density 
function . 

c)  If  gas  concentrations  are  examined  from  a sample  site  which  is 
positioned  such  that  the  gas  has  been  subject  to  the  dilution 
process  for  time  t,  then  the  steady  state  probability  distribution 
of  gas  concentrations  for  non-zero  concentrations  pe(c)=p(cp,t |cQ) . 
This  is  also  known  as  the  ensemble  average  in  the  literature. 

The  probability  density  function  for  gas  concentration  which  has  been 
derived  from  the  Wiener  process  is  too  complex  for  use  in  measurement 
experiments  and  it  will  be  approximated  by  the  Gamma  family  of  probability 
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density  functions.  This  family  was  chosen  because  it  has  as  special  cases 
the  exponential  and  deterministic  distributions  and  in  some  cases  (when  the 
probability  of  values  less  than  zero  is  small)  it  may  be  used  to 
approximate  the  normal  distribution.  The  Gamma  distribution  will  be 
referred  to  as  G(a,k)  and  its  general  form  is; 


k k-1  -qc 
a c e 


r(k) 


c > 0 
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For  this  application  the  distribution  becomes  G(a  ,k)  where 

ce  = expected  concentration  for  non-zero  concentrations  , and 
k = shape  parameter  of  the  Gamma  distribution. 


C ' 


k-1  - k 

c e 


Pe(c)  = 


c > 0 


r(k) 
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See  Figure  G.2  for  examples  of  the  Gamma  distribution  for  different  values 
of  the  shape  parameter  k.  Note  that  for  k=l,  the  resulting  distribution  is 
the  exponential  probability  density  function. 


: 


FIGURE  G.2  The  Gamma  Probability  Density  Function 
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III.  The  Meander  Process  Model 

The  plume  may  meander,  with  the  result  that  a sample  site  in  a 
stationary  geographical  location  will  not  experience  a continuous  plume, 
but  rather  will  experience  a series  of  intervals  when  the  plume  is  present 
and  then  not  present.  A two  stage  markov  process  is  used  to  model  the 
meandering  process.  Define 

Pij(t)  = The  transition  probability;  the  probability  that  the 

process  is  in  state  j at  time  t,  given  that  the  process  was  in 
state  i at  t=G. 

r-[j  = The  rate  of  transition  from  state  i to  i;  r^j  At  is 

the  probability  that  the  process  is  in  state  j after  some 
small  increment  of  time  At,  given  that  it  is  currently  in 
state  i. 

Let  the  state  i=0  represent  the  event  that  the  plume  is  not  present  and  the 
state  i=l  represent  the  event  that  the  plume  is  present,  then  it  can  be 
shewn  that 


10  r01  „~(r01  + r10)  Te 
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UU  e r01  r10  r01  r10 
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rCl  ~(r01  + r10*  Te 
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In  steady  state  as  t 

Po  = _2i_ 


and  PI 


r01 


Gil 


rCl  + no 


r01  + r10 


Consequently,  the  steady  state  probability  distribution  p(c)  can  be  written 


P(c)  = P0  + (l-Po)Pp(c),  G12 

where  Pp(c)  = the  steady  state  probability  that  the  gas  concentration  is  c 
given  that  the  plume  is  present. 

When  the  plume  is  present,  the  gas  is  not  a completely  homogeneous 
mixture,  but  consists  of  a number  of  discontinuous  eddies,  and  the  dilution 
process  within  each  eddy  is  assumed  to  follow  a Wiener  process.  In  steady 
state  the  probability  distribution  of  gas  concentrations  in  the  eddy  pe(c), 
will  be  distributed  according  to  the  Gamma  distribution  as  developed  in 
section  II  and  expressed  by  Equation  G9.  The  steady  state  probability 
distribution  of  gas  concentrations  can  now  be  written  as, 

P ( c ) = 1- y + YPe(c)  and 
where 

kc 


c* 


fi(c*)  = jp(c,t)dt  = I-Y+Y 

0 


Pe(c) 


pjc)  = 


KC'ck-l  g ce 


r (k) 


k =-L. 


and  1- y is  the  probability  that  an  eddy  is  not  present. 


Since  ft  (c*)  = 
ev 


'0 


pjc)dc  is  not  available  directly  a numerical  routine 


shown  in  section  V will  be  required  for  calculations. 
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IV  First  Passage  Times 

Define  {N(t);t>0}  as  a counting  process  for  the  number  of  events 
occurring  by  time  t.  If  the  arrival  rate  of  events  is  r;  r>0  and  if; 

a)  N(0)=C, 

b)  the  process  has  stationary  and  independent  increments, 

c)  P{N(h)=l }=rh+0(h)  for  h small,  and 

d)  P{N(h)>2=0(h). 

then  N(t)  is  a Poisson  process.  A function  is  of  order  0(h)  if 

lim  f(h)/h  = 0.  A process  has  independent  increments  if  the  number  of 
h ->0 

events  occurring  in  disjoint  time  intervals  are  independent.  The  increments 
are  stationary  if  the  number  of  events  occurring  in  the  interval  depends 
only  on  the  length  of  the  interval. 

It  is  argued  that  the  arrival  process  of  the  eddies  when  the  plume  is 
present  satisfy  the  conditions  specified  above  and  eddies  arrive  at  the 
sample  site  according  to  a Poisson  process  if  the  plume  is  present.  An 
important  characteristic  of  the  sampling  process  is  the  probability  that 
some  threshold  value  of  gas  concentration  c*,  has  been  exceeded  by  time  Te. 
In  the  study  of  stochastic  processes  this  is  called  the  first  passage  time 
distribution.  The  probability  of  threshold  exceedance  F(Te)  will  depend  on 
whether  the  plume  is  present  or  not  at  Te=0.  The  following  derivation  will 
assume  the  plume  is  present  at  time  Te=0  although  the  cases  of  plume  not 
present  and  no  information  can  also  be  dealt  with  using  this  type  of 
analysis. 

Consider  that  a threshold  exceedance  can  only  occur  when  the  plume  is 
present.  Using  an  absorbing  state  for  the  event  of  a threshold  exceedance  a 
three  state  Markov  process  can  be  defined.  It  is  important  to  note  that  if 
eddies  arrive  according  to  a Poisson  process  with  rate  r,  and  the 
probability  is  l~^e  (c*) that  the  gas  concentration  of  the  eddy  is  greater 
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than  c*,  then  the  process  of  threshold  exceedance  is  also  Poisson  with  rate 

r(l-ne(c*)). 

Defining  state  2 as  the  state  of  a threshold  exceedance  and  rj.2  as  the 
rate  of  threshold  exceedance  when  the  plume  is  present,  then  the  one-step 
transition  equation  can  be  written  as, 


w 


{Pq( t+  At)  ,P]_(  t+  At) , ?2( t+  At)  }=  { Pq( t ) ,Pi(t)  ,P2(t) } (Pij) 
here  (Pij}  is  the  one-step  transition  matrix  with  elements 


Pij 


G15 

defined  by 


1-roi  ^ rCl  At  Q 


(Pij)  = 


r10  At  l-(rio+ri2)  At  r\2  At 

0 0 1 


G16 


The  resulting  set  of  simultaneous  differential  equations  are; 

P0*(t)  = roiPi(t)  - riQPo(t) 

pl'(t)  = -(ri2+roi)Pl(t)  + riQPo(t)  G17 

P2’(t)  = ri2Pl(t) 

Examining  the  first  two  differential  equations  it  is  noted  that  p2.(t)  does 
not  appear  and  the  solution  to  this  smaller  set  of  simultaneous  first  order 
differential  equations  can  be  reduced  to  a single,  second  order  equation. 
Isolating  Pi(t)  in  the  resulting  second  order  differential  equation  has  the 
general  solution; 

pl(t)  = B]_  exp( b jt ) + B2  exp( b2t) 

where  b}  = [ -(rio+r01+r12)_s(lr ( (r01+r10+r12)^-4r01r  12)  ]/2 * and  G18 

b2  = [-(rio+r01+r12)+s<lr((r01+r10+r12)2-4r01r12)  ]/2- 
Bj  and  B2  are  found  from  the  boundary  conditions  for  P]_(0)  and  F2(0)-  Por 
the  case  where  the  exposure  interval  begins  with  the  plume  overhead  P}(0)=1 
and  Pq(0)=0.  Using  these  boundary  conditions,  a solution  was  found  for 
Equation  G18.  Noting  that  P2(Te)=F(Te)  is  the  probability  that  the 
threshold  has  been  exceeded  at  least  cnce  in  exposure  time  Te  and  is  equal 
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to  l-Po(Te)-Pi(Te) , then 

F(Te)  = l-(Bib2/roi)exp(-biTe)-(B2bi/roi)exp(-b2Te)  G19 

where  Bj  = l/2+{rio-roi+r 12)/ ( 2sclr( (r01+r10+r12)2~4r01r12) ) and  B2  = 1-Bp. 

A mathematical  model  of  the  discontinuous  dilution  process  has  been 
proposed.  Distributions  of  gas  concentration  in  an  eddy  as  a function  of 
time  since  emission  from  the  source  were  developed  from  the  Wiener  process 
model.  From  the  mathematical  model  the  steady  state  distribution  of  gas 
concentrations  at  a stationary  location  was  specified.  The  Gamma  family 
of  probability  density  functions  were  recommended  as  a an  approximation  to 
the  derived  gas  concentration  distribution.  An  expression  for  the 
probability  of  threshold  exceedance  for  some  exposure  time  Te  was  derived 
assuming  that  the  meandering  process  is  a two  stage  Markov  process  and  the 
arrival  process  for  eddies  when  the  plume  is  present  is  Poisson.  Numerical 
techniques  for  calculating  the  steady  state  probability  of  exceeding  a 
threshold  concentration  at  any  instant  for  the  Gamma  probability  density 
function  are  shown  in  section  V. 


V.  Numerical  Routines 


The  general  form  of  the  Gamma  used  in  the  report,  G f — , kl  where 


7 ls 


kc 

k 1 k-1  c 
c e p 


pp(c)  dc  = 


w 


dc 


(i  ‘I 


G20 


The  cumulative  distribution  is 


VV  =f 


kc 

k 1 k-1  c 
c e p 


r(k) 


dc 


G21 
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Normalizing  by  defining  a new  variable  3-  then 

cp  ’ 


C = c r-  and 

CP 


dc  = Cp  d 


c >1 


such  that 


3 


Pp  1-1  d 


, (0 


kc 

• k k-1 

k ^ r—  c ^ c — 

c_  - 


e "p  cpd 


Pi 


r(k) 


P; 


and 


k-1 


-kf- 

e cp 


r(kT 


k-1 


-k 


CP k 

r c \ c 

fr  e p 

' C ' 

= f 

‘-P  d rc  1 

J 

0 

^Tkl  - 

1 Pj 

When  provided  with  Cp  and  ip,  PD(C)  anc^  can 


Subroutine  For  r(k)  : 

For  k > 10  Stir!  ingsformul a will  be  used  : 


*n[r(k)]=  (k  - I)  zn  k - k + 1 tn  (&r)  + ~ ~ ~ 3 + 


G22 


G23 


G24 


be  found. 


J 

1260k5 


Then 
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r(k)  = exp  [nn  r(k)] 

and  for 

k < 10, 

r(k)  = r (a  + e)  . 

Where  a is  an  integer  and  0 < e < 1 . 

Therefore, 

r(a  + e)  = (a  + e " 1 ) (a  + e - 2)  x x (e  + l)(e)  r(e). 


where 


oo 


r(e)  = l c en 
n=l  n 


The  computer  listing  is  shown  in  Table  G1 . 


The  subroutine  for 


P 


P 


c 


P- 


G25 


where  x = — # is  found  from 
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dn  [p, 


'-§-]]  = in  k + (k  - 1)  an  (x)  - x - dn  (r(k))  , 

•cpJ 


and 


exp  [dn  [pp  ' ]]  . 


G26 


The  subroutine  for. 


k 


°p  r v 

p 0 


(k)k  f^-1 

e cp  d 

IcpJ 

G27 


is  proposed  by  BURGIN  [1975],  and 

ck  kc^~  ‘ c . ck 
e pa  — 


■ 


Ri 


r(k) 


kc 

c 


P xk-l  0-x 


e dx 


r(k) 


G28 


“p  f#-] 

P 


kc 

. k-1 

kc^  c 

e p 


r(k) 


l 

n=0 


n+1 


if  (k  + r) 
r=0 


G 29 


Computer  listings  for  Pp(c)  anc*  ftp(c)  are  s^own  ln  tables  G2  and  G3. 


-616- 

Table  Gl  Subroutine  for  r ( k ) 


10  PRINT  "THIS  PROGRAM  CALCULATES  VALUES  FOR  THE  GAMMA" 

20  PRINT  "FUNCTION  WHEN  SUPPLIED  WITH  THE  ARGUEMENT  K" 

30  PRINT  "WHAT  IS  THE  VALUE  OF  K" ; 

40  INPUT  K 
50  DIM  C < 1 2 ) 

60  FOR  L= 1 TO  12 
70  READ  C(L) 

80  NEXT  L 

90  IF  K <=  10  GOTO  120 

100  LGAMMA=  ( K — . 5 ) *LOG  <K>'~K+.  5*L0G  (2*3.  1 4 1 6 ) ( 1 / < 1 2*K ) ) - ( 1 / (360*10  ) 

110  GOTO  270 

120  LGAMMA=0 

130  FOR  L=1  TO  10 

140  K1=K-L 

150  IF  K 1 =0  GOTO  190 
1 60  LGAMNA=LGAMMA+LOG ( K 1 ) 

170  IF  KK1  GOTO  190 
180  NEXT  L 
190  K2=0 
200  K3= 1 

210  IF  K1=0  GOTO  260 
220  FOR  L= 1 TO  12 
230  K2=C (L) * <K1 ‘ L) +K2 
240  NEXT  L 
250  K3= 1 / K 2 

260  LGAMNA=LGAMNA+LOG (K3) 

270  GAMMA=EXP ( LGAMMA ) 

280  PRINT  "THE  VALUE  OF  GAMMA (K)  IS"; GAMMA 
300  END 

310  DATA  1.0,. 57721 , 65587, 042, . 16653, -.04219. -.00962 

320  DATA  7. 2169E— 03, 1. 1652E-03.2. 1524E-04, 1 . 2305E-04 , 2. 0 1 35E-05 


Table  G2 
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Subroutine  for  the  Gamma  Distribution  Pp(c) 


PRINT 

"WHAT 

IS 

INPUT 

IP 

PRINT 

"WHAT 

IS 

INPUT 

R 1 

PRINT 

"WHAT 

IS 

INPUT 

R2 

PR  I NT 

"WHAT 

IS 

> INPUT 

ITE 

> K=1/(IPA2) 

IS  THE  STEP  SIZE  FOR  THE  NORMALIZED  CONCENTRAT I ONS " ; 


10  PRINT  "THIS  PROGRAM  CALCULATES  VALUES  FOR  THE  NORMALIZED" 

20  PRINT  "GAMMA  PROBABILITY  DENSITY  FUNCTION  WHEN  GIVEN  THE  " 

30  PRINT  "CONCENTRATION  FLUCTUATION  PARAMETER  i AND  THE  LENGTH1 
40  PRINT  "AND  INTERVAL  OF  THE  NORMALIZED  CONCENTRATION  VALUES" 

0 PRINT  "WHAT  IS  THE  VALUE  OF  i 
(60 

THE  LOWER  VALUE  FOR  THE  NORMALIZED  CONCENTRAT I ON " 
BO 

THE  UPPER  VALUE  FOR  THE  NORMALIZED  CONCENTRAT I ON " 

'86 
I0 
100 
110 

la  20  GOSUB  1040 
130  FOR  L=R 1 TO  R2  STEP  ITE 

140  LNG= <K— 1 ) * (LOG (L) -LOG ( IP'2> ) -L*K-LGAMMA 

145  LNC=LNG-LOG ( IP'2) 

146  CONC=EXP (LNC) 

0 LPRINT  L , CGNC 

160  NEXT  L 
170  END 

1040  DIM  C ( 1 2) 

050  FOR  L= 1 TO  12 
060  READ  C (L) 

1 070  NEXT  L 

1080  IF  K <=  10  GOTO  1110 

090  LGAMMA= (K-. 5) *LOG  <K) -K+. 5* LOG (2*3.  1416)  + ( 1/ ( 12*K) ) - ( 1 / (360*K) ) 

100  GOTO  1260 

110  LGAMMA=0 

120  FOR  L= 1 TO  10 

130  K 1 =K— L 

140  IF  K 1 < =0  GOTO  1180 
1 50  LGAMMA=LGAMMA+LOG  < K 1 ) 


1 160  IF  K 1 < 1 GOTO  1180 
170  NEXT  L 
1 180  K2=0 

1185  IF  K1<0  THEN  K1=K 
190  K3= 1 

|200  IF  K1=0  GOTO  1250 
210  FOR  L= 1 TO  12 
220  K2=C(L)*(C1-  L)+K2 
1230  NEXT  L 
240  K3=l/K2 

i 250  LGAMMA=LGAMMA+LOG (K3) 

260  RETURN 
i 280  STOP 
i 1290  END 

1300  DATA  1.0,.  57721  , -.  65587 042 , . 16653,  -.  04219,  -.  00962 

|j310  DATA  7 . 21 89E-03, 1. 1652E-03,2. 1524E-04, 1 . 2805E- 04 , 2 . O 1 35E-05 
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Table  G3  Subroutine  for  fl(c)  and  v(c) 


10  PRINT  "THIS  PROGRAM  CALCULATES  THE  PDF  AND  CDF  OF  THE  NORMALIZED" 

20  PRINT  "GAMMA  PROBABILITY  DISTRIBUTION" 

30  PRINT  "WHAT  IS  THE  VALUE  OF  THE  FLUCTUATION  INTENSITY  i"? 

40  INPUT  IP 

50  PRINT  "WHAT  IS  THE  RANGE  OF  VALUES  OF  THE  NORMALIZED  PARAMETER?  MIN. 
60  I NPUT  RR 1 . RR2 

70  PRINT  "WHAT  IS  THE  STEP  SIZE  ACROSS  THE  INTERVAL"? 

80  INPUT  IITE 

90  PRINT  "WHAT  IS  THE  ACCURACY  REQUIRED  FOR  THE  NUMERICAL  ROUTINE"? 

100  INPUT  EPS 

110  PRINT  "WHAT  IS  THE  VALUE  OF  THE  I NTERMI TTENCY " ? 

120  INPUT  MIT 

130  LF'RINT  "THE  FLUCTUATION  INTENSITY  IS"?  IP 
135  LPRINT  "THE  I NTERM I TTENCY  IS"? MIT 
140  LPRINT  "CONC". "CDF". "VCDF" 

150  R 1 =RR 1 *M I T 
160  R2=RR2*MIT 
170  I TE= I I TE*M I T 
180  K=1/(IP''2> 

190  GOSUB  450 
200  SUM 1=0 

210  FOR  L=R 1 TO  R2  STEP  ITE 
220  SUM=0 
230  N=0 

240  LSUM= (N+l ) *LQG(L*K> 


250  FOR  R=0  TO  N 
260  LSUM=LSUM-LOG (K+R) 

270  NEXT  R 

280  SUM=EXP (LSUM) +SUN 

290  IF  (EXP (LSUM) /SUM )< EPS  GOTO  320 
300  N=N+ 1 


310  GOTO  240 
320  GOSUB  410 
330  GAM=EXP (LNG) *SUM 
340  SUM l=SUMl+( CONC* ITE) 

350  CUM= 1 — M I T +MIT*GAM 

360  VCUM= 1 —CUM 

370  L 1 =L/MI T 

380  LPRINT  Ll.CUM.VCUM 

390  NEXT  L 

400  END 

410  LNG=  (K-l  ) * (l.OG  (L)  -LOG  ( IP  '2)  ) — L*K— LGAMMA 
420  LNC=LNG-L06 ( IP  '2> 

430 . CONC=EXP (LNC) 


440  RETURN 
450  DIM  C ( 1 2 ) 

460  FOR  L= 1 TO  12 
470  READ  C ( L ) 

480  NEXT  L 

490  IF  K<  = 10  GOTO  520 

500  LGAMMA=  (K— . 5)  *LOG  (K)  -K+.  5+-L0G  (2*3.  1416)  -M  1/  <12*K)  ) - ( 1/  \360*K)  ) 

510  GOTO  680 

520  LGAMMA=0 

530  FOR  L= 1 TO  10 

540  K 1 =K— L 

550  IF  K1<=0  GOTO  590 


560  LGAMMA=LGAMMA+LOG ( K 1 ) 


570  IF  K 1(1  GOTO  590 
580  NEXT  L 


590  K2=0 

600  IF  K 1 < 0 THEN  K1=K 


MAX" 


620  IF  K 1 =0  GOTO  670 
630  FOR  L= 1 TO  12 
640  K2=C(L)MC1  D+K2 
650  NEXT  L 
660  K3=l/K2 

670  LGAMMA=LGAMMA-t-LOG  (R3) 

680  RETURN 

69  O DAT  A 1 . u , . 5772  1 . — . £>5507  , — . 042 . . 1 6653 . —.04219.  — . 00962 

700  DATA  7.21 B9E-03 , 1. 1652E-03.2. 1524E-04, 1 . 2805E-04 . 2 . 01 35E-05 
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